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XIV. BAND ZWEITES HEFT 1956 


Die pulsierende Quelle unter der freien Oberflaiche 
eines Stromes endlicher Tiefe 


Von E. Becker 


1. Ubersicht. Die folgende Arbeit behandelt das zweidimensionale Problem einer Quelle, die 
nter der freien Oberflache eines Stromes endlicher Tiefe pulsiert. Im allgemeinen laufen in einiger 
utfernung von der Quelle vier voneinander verschiedene Wellenziige von dieser weg, davon einer 
1 Gegenstromrichtung, die drei andern in Mitstromrichtung. Fir diese Wellen wird die Wellen- 
inge durch verhaltnismafig einfache, an das Dispersionsgesetz der Wellen ankniipfende Uber- 
gungen bestimmt. Weiterhin werden Formeln fiir die Wellenamplituden angegeben. Wenn man 
ut den Werten von Pulsationsfrequenz und Stromgeschwindigkeit gewisse Grenzen iiberscbreitet, 
ann verschwindet die stromaufwarts wandernde Welle und ebenso eine der stromabwarts wan- 
ernden Wellen. Die Formeln fiir Wellenlange und Amplitude werden besonders einfach fiir 
men unendlich tiefen Strom.! Als Grenzfalle werden hier die in ruhendem Wasser pulsierende 
Juelle 2 und die zeitlich stationare Quelle * erwahnt. 


2. Aufgabenstellung und Grundgleichungen. Im Abstandeh unter der Oberflache einer mit 
er Geschwindigkeit c iiber horizontalem Boden gleichférmig strémenden, schweren Filiissigkeit 
efinde sich eine Quelle, deren Starke Q sich pe- 
odisch mit der Zeit andert: 

Q=Qen! (V 
Abb. 1). Die im ungestérten Zustande ebene Ober- 
ache des Stromes wird durch die Wirkung der 
uelle verformt. Durch die nachfolgenden Rech- 
ungen sollen diese Verformungen bestimmt werden 
nd zwar insbesondere fiir groBe Entfernungen von 
er Quelle, wo sich harmonische, von der Quelle 
yrtlaufende Wellen ergeben. 

Unter der Annahme wirbelfreier Str6mung und 
ei Beschrankung auf kleine Oberflachenverformung 


uf das Strémungspotential ® folgenden Bedingungen geniigen (zur Koordinatenbezeichnung 
gl. Abb. 1): 


Abb. 1. Bezeichnungen, Koordinatensystem. 


PD PD 
ape gee: (2) 
ir z= 0 mu gelten ? 
a a\? Sa 
Rg] OF 5G —°- () 
‘ir z = — H (Boden) muB gelten 
OD 
oz Y; (4) 
nd in Umgebung der Quelle muf sein 
@ = — 2 Inr + regulare Funktion. (5) 


1 Zusatz bei der Korrektur: Wahrend der Drucklegung dieser Arbeit wurde dem Verfasser eine 
issische Veréffentlichung bekannt, in der auf etwas anderem Wege die Lésung fiir diesen Sonderfall 
benfalls gefunden wird: M.D. Khaskind, Uber Wellenbewegungen schwerer Flissigkeiten (in Russisch). 
rikl. Math. Mekh. 18 (1954), S. 15. 

2 H. Holstein, Z. angew. Math. Mech., 17, (1937), S. 38. 

3 H. Lamb, Hydrodynamics, 6th ed., Cambridge, 1953; Th. Havelock, Proc. Roy. Soc. A, 115, (1927), 
. 274. 

4 Vgl. Lamb, a. a. O. S. 364; die dort angegebene Bedingung ist als Grenzfall c = 0 in (3) enthalten. 
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AuBerdem miissen gewisse Ausstrahlungsbedingungen im Unendlichen hbefriedigt sein, die erst: 


spater formuliert werden sollen. ea | 
Wenn man eine den Bedingungen (1) bis (5) geniigende Lésung gefunden hat, erhalt man die» 


Oberflachenverformung ¢ (Auslenkung aus der horizontalen Ruhelage) als 1 


i /(G 0 | 
es 6 
3 g (i eae : ( y 
Da die Quelle harmonisch mit der Frequenz w pulsiert, machen wir den Ansatz 
Diez @(X, 2) Or 6 as (7)) 


hierbei ist c x das Potential der ungestérten Parallelstrémung. Setzt man den Ansatz (7) in die > 
Randbedingung (3) ein, so wird > 
—wrg—2iwed, +O Gxx + SY: = 9 (8) 
Tees. 0h | 

3. Harmonische Wellen als Lésungen der Grundgleichungen. Bevor das durch die Gleichungen | 


(1) bis (8) gegebene Randwertproblem gelést wird, sollen harmonische Wellen als partikulare } 
Lésungen der Gleichungen (2), (3), baw. (8), und (4) betrachtet werden. Hierzu machen wir den) 


Ansatz 


a (9)) 
P, = et iO {k > 0). (10) 


Der Ansatz (9) entspricht Wellen in Gegenstromrichtung, (10) in Mitstromrichtung. Beide An-} 
sitze befriedigen die Gleichungen (2) und (4); die Konstante 1/€o) k H wurde zur Vereinfachung ; 
spaterer Rechnungen eingefiihrt. Die Bedingung (3) 

“VA/rA haw. (8) ergibt | 


| 


fiir gy: fi(k) =(o + ekP—gkIgkH =O, (11) 
fir y,: £,(k) = (@—ckh)??—gkigkH=0. (129) 


WF -——"_ ya) ss Wir:« diskutieren zunachst Gleichung (11). Unter | 
der Voraussetzung k + 0 erhalt man 
U, 

WA, (g+e)-f%eH=VH. Cay 

27 (4<c) (4>0) 
y V(k) ist hierbei bekanntlich die zur Wellenzahl k 
~ gehérige Phasengeschwindigkeit bei Ausbreitung der 
as Wellen auf ruhendem Wasser 7. Wenn wir der An- 


zi schaulichkeit halber von der Kreisfrequenz w und 
3 der Wellenzahl k zu der Frequenz vy = /2 7 und der 
\ Wellenlange 1 = 2 7/k tibergehen, erhalten wir aus | 
Abb. 2. Dispersionskurve V/(A) und Ubersicht iiber die Lésungen (13) nach Wurzelziehen und einfacher Umordnung 


der Gleichungen (15), (18), (19). in Va) Fi (14) 
—VA=—C. 


Das untere Vorzeichen kann hier zu keiner Lésung fiihren, da definitionsgema8 alle in (14) vor- 
kommenden Gréfen positiv sind. Man erhalt daher 


Vy on (15) 


Die Lésungsméglichkeiten von (15) lassen sich leicht an Abb. 2 iibersehen. AuBer der ,»Disper- 
sionskurve’ V(A) ist hier V(A) —vA und V(A) + vA (vgl. weiter unten) aufgetragen. Man erkennt, 
daB es entweder zwei einfache Lésungen (mit J und 2 bezeichnet) gibt, oder eine Doppellésung | 
oder gar keine. Die beiden letzten Falle sollen zunachst auBer Betracht bleiben und der Fall zweier 
reeller Lésungen soll zuerst betrachtet werden. 
Zur Unterstiitzung der Anschauung schreiben wir (15) in der Form 
Viyece | 


1 Vel. Lamb, a. a. O., S. 364. 
2 Vel. Lamb, a. a. O., S. 366. 
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In dieser Schreibweise erkennt man leicht, daB die Lésungen von (15) Wellen entsprechen, deren 
Phase relativ zum ruhenden Wasser mit der Geschwindigkeit V stromaufwarts lauft. Da in beiden 
Fallen V>c ist, laufen die Wellen mit V —c an einem ruhenden Beobachter stromaufwarts 
vorbei. Gleichung (16) ist dann die Bedingung dafiir, daB die vom ruhenden Beobachter fest- 
gestellte Frequenz (Zahl der Wellen, die in der Zeiteinheit an ihm vorbeiziehen) den vorgegebenen 
Wert v hat. 

Man kann weiterhin leicht zeigen, daB die Lésung I einer Welle entspricht, deren Gruppen- 
geschwindigkeit U> c ist, wahrend bei Lésung 2 U <c ist (siehe Abb. 2). Nur Lésung J ent- 
spricht daher Wellen, die imstande sind, Energie stromaufwarts zu transportieren, denn U ist die 
Geschwindigkeit, mit der Energie relativ zum ruhenden Wasser transportiert wird. 

Die Funktion f,(k) wird im folgenden in der Form 


Silk) = (k — hy) (& — hes) fio) | 
mit ‘ (17) 
fie(k) + 0 a Oe J 

geschrieben werden. 
Gleichung (12) fiihrt in ganz analoger Weise auf folgende beiden Gleichungen: 


Vac (18) 
V+tyvAj=ote. (19) 


Der Anschauung entnimmt man sofort (Abb. 2), da® (18) und (19) fiir alle Wertepaare c, » jeweils 
eine einfache Lésung haben. Die Lésung 3 von (18) gibt Wellen, die sich relativ zum Strom mit 
der Geschwindigkeit V, absolut also mit V + c stromab bewegen; Gleichung (19) gibt Wellen 4, 
deren Phase sich mit V <c relativ stromauf, absolut also mit c — V stromab bewegt. Die 
Funktion f,(k) 1aBt sich immer in folgender Form schreiben: 


Silk) = (k —e ks) (k ar ky) fza(k) 
mit (20) 
faa(k) + 0 in Oi Cog 


4, Bestimmung des Potentials aus den Grundgleichungen. Nach dieser Vorbereitung 1aBt sich 
die in Ziff. 2 gestellte Aufgabe folgendermafen lésen: Das Potential p wird aus drei Bestandteilen 
zusammengesetzt : er 

P= 4 G4 G. (21) 
Das Glied y) ist ein geeignet gewahltes Potential freier Wellen der oben be- 
sprochenen Art. Das Glied py’ ist das Potential der in Abb. 3  skizzierten 
Quell-Senkenanordnung; es geniigt den Bedingungen (2), (4) und (5). Fir otf 
eo 0 ist yk 
p? = =o =0, (22) 
Das Glied p ”) tragt somit nichts zur Oberflachenverformung [nach (6)] bei. Da a 
p der Oberflachenbedingung (8) geniigen mu, erhalt man fiir ~ unter Benutzung l sp 


von (22) [py geniigt (8) identisch | 7, 
— oO + 2iweget Oba ter =—B py (23) 
mors == 0). 
Der Ausdruck "|, 1aBt sich leicht bestin.men. Es ist bekanntlich orf) 
‘ Qo Sur enH+h) Si)" @nH—h) y! 
a9. 5%| > x? + (2nH +h)? Pa 2+ (2nH—h) |° Ge 


Die einzelnen Summanden in (24) sind jeweils die von einem symmetrisch zur 


Oberflache gelegenen Quell-Senkenpaar an der Oberflache erzeugten Vertikal- 4yy.3. Brzeugung des 
Quellpotentials gp) 


geschwindigkeiten (Abb. 3). ToL QUT dopoh oka ge. 
Fir Pp wird folgender Ansatz gemacht: polegplice Otel 


Sait t: | Sea Golk Gee HY 
B= | Vo tylko) eg ae (25) 
0 
Damit sind (2) und (4) befriedigt. Einsetzen in die Bedingung (23) fiihrt mit (24) auf die Glet- 
6* 
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hung 
J Epa(h) fa(0e) e+ al) ful) e—* FY dk (26) 
pase Qn Hh) Gly rear |-| 


~ 8 | a t+ ent hy oy + nH —hp 
Hierin sind f,; und f, durch (11) und (12) definiert. Gleichung (26) wird nun in folgender Form ge- 
schrieben: 


oC 


[ po(k) fo(k) ek * dk i) ‘ps (— k) f,(— k) et *? dh 


| Qo iN CI" @nA+h) NS (— 1)" (nH —h) 
<6) Spee pu x + (2n H—h)? | 


) 


Zur Bestimmung von y, und y, aus (27) benutzen wir ein bekanntes Ergebnis aus der Theorie | 


der Fourierintegrale: Wenn 


; tk x _ 8% ee ee 
| Pwe dk, =O 


ist, dann gilt 
F(k) = 220 e—elel , 
Xt 


Wendet man dies auf jeden einzelnen Summanden in (27) an so wird 


mh _ _ 8% 5) (1)? eth pe 2s Ras 1)" ehh g—2 nk mes Ise 28 
valk) = 907 | 2,(—) P= 1 ). 28) 
Die Summen in (28) sind jetzt einfache geometrische Reihen und kénnen leicht aufsummiert 
werden mit dem Ergebnis 
_— gQ Gof k(H—h) 
Yl) a7) Cok 


Cel ee? (29) 
Damit wird 
r ikx rt neg © 
pe | ee 
290 ; So(k) Sik) | Sof?k H 


Benutzen wir nun fiir f; und f, die Darstellungen (17) bzw. (20) und nehmen gleichzeitig eine 
einfache Partialbruchzerlegung vor, so ergibt sich nach elementarer Rechnung 


foo) 


(30) 


shee Ostic| dew ee 1 1 F eee ( ie, aL Cojk(H—h)Gojk(H+-2) 7, 
ei F fash) (kg — hy) \h—keg bh — key Sia(l) (ky, — he) \k—k, kk Gof? k H ai 
Die hier auftretenden Integrale sind von der Form 
Paved | Retiha we o in 0<k<a@ 
| B(k) <———dk mit B(h) | ae 
0 
Es 1aBt sich leicht zeigen, daB folgendes gilt 1. 
2 ' } 
pad = +71 Blk) e'** 1 R, fir x> 0, 
| OS ae a} —mni Blk.) eb** 4 R, fir x <0, | 
mis \ (32) 
: ee [= — i B(k,) e—i*s* + R, ivigeee ssl (Ns 
| BU) eee CPA TIA ko ee eae tn 
0 J 


' Vel. z. B. Lamb, a. a. O. S. 412. Die dort gegebene Umformung des Integrals lat sich hier allerdings 
nicht in gleicher Form anwenden, da der Integrand in (31) auf der imaginaren k-Achse singular wird. Fiir den 
Beweis der folgenden Formeln geniigt es aber, den Integrationsweg z. B. auf einer Parallelen zur reellen Achse 
zu schheBen, die fiir kleine |k| in beliebiger Weise in den Nullpunkt gefiihrt wird, aber derart, daB vom In- 
tegrationsweg keine Singularitaten umschlossen werden. 


(27) | 
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vobei lim “RK, = 0 es bag as Ba 
|*| + 
Da wir uns, wie in der Einleitung angekiindigt, nur fiir die wellenartige Verformung in einiger 
Entfernung von der Quelle interessieren, wird hier auf eine Diskussion der R,, verzichtet. Eine 
jolche Diskussion kann nach den iiblichen Routinemethoden der Analysis vorgenommen werden !. 
Wir fiihren nun die folgenden Abkiirzungen ein: 
5 he j 
Bys(k; 2) — SOUR HA) Goh e+ 2) 
Fi(k) Co? k 1 
ind B,,(k;) ganz entsprechend. Dann erhalten wir unter Berwcksichtigung von (32) fiir x > 0 
OG t 1 re ee: 
We oe ee Fees [Bs,(kgs 2) e'*s * — Ba,(kys 2) ef be *] 
= 1 
k, —k, 


(33) 


, (34) 
[Byo(kys 2) e's * — Byo (kgs 2) € va ++ R | 
J 


mit lim gp = 0. Die Lésung fiir x <0 erhalt man aus (34) durch Vorzeichenvertauschung, 


ndem man also k, —k, durch k, —k, und entsprechend k,—k, durch k, —k, ersetzt. Freilich 
indert sich fiir x < 0 das Glied Pr in der Lésung, doch soll hier nicht darauf eingegangen werden. 

Die Wellen mit der Wellenzahl k, sind nach Ziff. 2 die einzigen Wellen, die Energie strom- 
aufwarts transportieren. Die in Ziff. 3 angekiindigte Ausstrahlungsbedingung befriedigen wir 
laher, indem wir iiber @ ein Wellenpotential gy" derart ttberlagern, da® fiir x 0 nur die Wellen 
mit k —k, ibrigbleiben, fir x <0 nur diejenigen mit k = k,, kg, ky. Nach (21) ist dann noch das 
Quellenpotential py zu addieren, das jedoch nichts zur Oberflachenverformung beitragt. Wir er- 
2alten dann fiir x > 0 


PEs ors Titi ea a 1 ( 

p= — 4 Opi ED oie 4 gh + gt (35) 
Geay Ss 
und fir x < 0 

Ber 2s 5 Le Oe ey 7 a | Ae ee 

P= rey ga(kg; 2) e' *s* — Bay(ky; 2) ek *] 
Qy i , a 

Pigs By (Kg, 2) e—'* + pR + 9% . | 


wobei nae P= =0. 

Die Gleichungen (35) und (36) gelten fiir den Fall, daB (15) zwei voneinander verschiedene, reelle 
Lésungen hat. Wahlt man w und ¢ derart, da (15) gar keine reelle Lésung mehr hat, da8 also 
i(k) += 0 in 0 <k <oo, so verschwinden in der Darstellung (35) und (36) die zu k, und k, gehérigen 
Wellen. Stromauf werden in diesem Falle gar keine Wellen mehr laufen und stromab nur noch die 
nuk, und k, gehérigen. Das Zusammenfallen von k, mit k, (also der Fall, daB in Abb. 2 die hori- 
zontale Gerade V = c die Kurve V — yA in ihrem Maximum beriihrt) bedeutet U = c, da U>c 
st fiir k, und U <c fir k, bei k, + k,. Andererseits erkennt man dann an (35) und (36), dab 
unsere linearisierte Theorie iiber alle Grenzen wachsende Wellenamplituden ergibt. Dieses Ver- 
halten ist typisch fiir U = c. Es tritt bekanntlich auch dann ein, wenn eine in ihrer Wirkung zeit- 


mabhangige Oberflachenstérung mit der Geschwindigkeit c = | g H (Grundwellengeschwindig- 
seit) iiberruhendes Wasser bewegt wird”; denn dann ist V = U = c = J/g H fiir die hinter der Stérung 


sich ausbildenden stehenden Wellen. Dieser Fall kann auch bei unseren Lésungen fiir wm — 0, ¢ = VeH 
sintreten, worauf hier ausdriicklich, aber ohne weitere Diskussion hingewiesen sei. 


5. Spezialisierung auf die in ruhendem Wasser pulsierende Quelle und die langsam pulsierende 
Juelle. Wir betrachten jetzt den Grenzfall c =0, m+ 0, also die in ruhendem Wasser pul- 
jerende Quelle. Es ist dann k, = ks, k, =k, = 0. Man erhalt fiir x > 0 


Bio (ky; 2) 
ky 


p=—£i heb pt +o. (31) 
fiir x < 0 ist die Lésung symmetrisch hierzu [(35) und (36) unterscheiden sich dann zwar formal um 
ine Konstante, aber das andert nur den willkiirlichen Nullpunkt des Potentials]. Die Deutung von 


‘ Vel. Lamb, a. a. O. oder Holstein, a. a. O. 
2 Vel. Lamb, a.a.O. S. 409, 414. 
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(37) ist klar: nach beiden Seiten laufen von der Quelle harmonische Wellen mit der Wellenlange 
A = 27 /k, weg. , 

Der andere Grenzfall ist c + 0, w+ 0. Es gehen dann k,, k; und die zugehérigen Ober- 
flachenverformungen ebenfalls > 0. Wir betrachten im folgenden daher nur die zu k, und k, 
gehorigen, nicht verschwindenden Potentiale. Fiir w = 0 sei k, = ky = ky. 

Dann wird fiir w— 0 

ky = ky — Ak +--+ . 
fk epee (eee) 


wobei lim 4k = 0 und die Punkte Glieder héherer Ordnung andeuten. Das Potential @ wird dann 


o> 
nae AB (hac eAr ey ere | 
ons 5 mm, ea tO t ema NCO) ! ot(k, + Ak) x 
2 a 6 Q% sci ky + 4k — kg : \s (38) 
2 Bs. (he ks ee 6 VER 
imt~34\%0 “7 p— i(k, — Ak) x u A 
+ g Qt e ko Ak k, e 7 Pr | 


In ®z sind hier sowohl diejenigen Anteile des Potentials zusammengefaBbt, die mit wachsendem |x| 
gegen Null gehen, als auch solche, die mit Ak baw. k, und k, verschwinden. Spaltet man auch in den 
beiden ersten Summanden in (38) die mit Ak bzw. k, und k, klein werdenden Anteile ab, so erhalt 
man nach kurzer Rechnung 


D's =i2ig 0, Lele) *| aca ane aces (39) 
0 


hierin hat w/Ak einen wohldefinierten Grenzwert. Zu seiner Bestimmung kniipfen wir an Abb. 4 
an: 


2rd = 2 [Aaj 


Nach der Definition von k und 7 ist andrerseits 


Gay 
AV on Ak, 
also 
7) 
d, 
Benutzt man die Definition der Gruppengeschwindigkeit U 
dV 
Ue V —) a? (41) 
so erhalt man aus (40) 
w 
Fp ee 
oder,.day))—= cltun mi Ost. 
: w 
ai. Ak = ¢ — U . (42) 
Gleichung (39) geht dann iiber in 
b— ex = 250, uleiad felt =a) sin ky x + Dp (43) 


0 


fiir x <0. Dieses Ergebnis ist einfach zu deuten: Fiir w = 0 erhalt man stehende Wellen strom- 
abwarts von der Quelle mit der Wellenlange A) = 22/k), wobei 2) bzw. ky durch die Bedingung 
V(Ao) = ¢ festgelegt sind. Wenn nun » =+- 0, aber noch geniigend klein ist, so bleiben die stehenden 
Wellen im wesentlichen erhalten, ihre Amplitude wird aber mit der Kreisfrequenz w moduliert. Be- 
merkenswerterweise kommt diese Modulation in einer Entfernung |x| von der Stérquelle mit einer 
der Signalgeschwindigkeit ¢ — U entsprechenden Phasenverschiebung an. Dies ist einleuchtend, 
da c — U diejenige Geschwindigkeit ist, mit der die Energie in dem stehenden Wellenzug (relativ zu 
einem ruhenden Beobachter) stromabwarts transportiert wird. 


6. Spezialisierung auf einen unendlich tiefen Strom. Es ist wohl angebracht, hier darauf hinzu- 
weisen, da man bei numerisch vorgegebenen Werten von H, h, wm und c die Zahlenwerte von kg 
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9 ky, ky und f,.(k), f54(k) unschwer ermitteln kann; dann sind nach (33) auch die Funktionen B,, 
zw. Bz, bekannt. Am einfachsten diirfte fiir die numerische Bestimmung ein graphisches Verfahren 
ein, wie es durch Abb. 2 nahegelegt wird. 


Fur den Sonderfall H = oo kann man leicht analytische Ausdriicke fiir k,,, fix, Bi, finden; es 
vird dann nimlich 1 


fo tu = ©, (44) 


= By = {eh 9), (45) 


ieee = 2 w c h = yi j (Reo )] (46) 
Ce | 


“ 2 g—2wc 
ge+2wc Ei oe 
jet BOs ene ee ip 47 
3.4 9 2 1 = g aie 9we ( ) 
Jie Wurzeln k,,, sind reell und voneinander ver- Me 
chi ve 
chieden, wenn wee 
& 
CO) Cpa 48 
4 (48) 2500 — oat 
st. Diese Bedingung 1aBt sich iibrigens leicht an 
edem nicht zu schnell flieBenden und nicht zu yy 
eichten Bach nachpriifen, dessen Oberflache man 
veriodisch, z. B. durch periodisches Eintauchen ““ | 
ines Stockes stért. Man stellt dann leicht fest, £50 | 
la nur dann Wellen direkt stromaufwarts laufen, 
venn Bedingung (48) erfiillt ist. ie 
A 
Y(A) +vA 50 sie As 
at 103) 
V-LA a5 ica 
ES joe 
o-V 4 i coe af A; 0 
& Wee sab s iS i} 
Ay 4 Ay—~ 0 
6, | 
wo 
Ap g a2 a 06 08 10 
Abb. 4. Zuammenriicken der Lésungen 2 und 4 Abb. 5. Lange der Wellen I bis 4 (in ‘dimensionsloser Darstellung) 
ee BA eRoTE Pliny Sel) E in Abhingigkeit von @ e/g. 


Durch (44) bis (47) sind alle in (35) und (36) vorkommenden, von k,, abhangigen GréBen durch 
» und c ausgedriickt. In Abb. 5 sind fiir die vier von der Quelle weglaufenden Wellenziige die 
Vellenlangen und in Abb. 6 die Amplituden in Abhangigkeit von w c/g dargestellt. Fiir den Grenz- 
all c = 0 erhalt man nach einer einfachen Rechnung 


a 


iene tl[mt+ —«x 5 49 
P—cx=Qie ee eden 2) 
ir x > 0. Hieraus ergibt sich als Oberflachenverformung (6) 
hw? > 
C= — Qo F cos(me + er) te (50) 
& 


nit limég = 0. Fir x < 0 ist symmetrisch hierzu. Gleichung (50) stimmt mit dem Ergebnis von 
Tolstein uberein. 14 : 
Fiir w = 0, d. h. fiir die stationare Quelle in einem unendlich tiefen Strom erhalt man 
2 aS 
6=—-2Q,e “cosh x +Cr (51) 


1 Vol. FuBnote 1 von S. 69. 
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fiir x <0; lim Cr=0. Fiir x> 0 fallt die wellenartige Verformung aus. Aus (51) 1aBt sich leicht}} 


x—— oo 


die Lésung fiir einen Dipol herleiten, die von Bedeutung ist fiir die Umstromung eines eingetauchten 
Zylinders !. 


7. Ausblick. AbschlieBend sei bemerkt, da® man selbstverstandlich durch Integration itber| 
linien- oder flachenhafte Quellverteilungen die betrachteten Lésungen verallgemeinern kann. Ab-) 
gesehen von méglichen Interferenzeffekten, die zwar die Wellenamplituden, nicht aber die durch ‘ie 


14500 


Welle 1: ey ar a 


Welle. 3: Xf NS 


hs | WOURPS NS Nema 


Welle 4: = [P\PPoo-- | 


Pfeile deuten die Richtung an,inder sich 
ale Wellenphase roriptianzt 


0500}- a | | 7 
-0 P 
N 
@c 

| | ee V-vA 

C OY } 08 10 ee . “ 
“ : ae Abb. 7. Ubersicht iiber die zu wellenférmiger Oberflachen- 
Abb, 6. Amplituden der Wellen 1 bis 4 (in dimensionsloser Darstellung) verformung fiihrenden Lésungen bei Beriicksichtigung der 

in Abhangigkeit von w c/g . Oberflachenspannung, 


Uberlegungen in Ziff. 3 gegebenen Wellenlangen beeinflussen kénnen, ergibt sich jedoch fiir die 
Wellenausbreitung, die in dieser Arbeit besonders interessiert hat, hierdurch nichts Neues. Es wird 


in diesem Zusammenhang darauf hingewiesen, das das Problem einer pulsierenden Quelle von Be- | 
deutung ist fiir die Behandlung der durch die Schwingungen eingetauchter Kérper erzeugten Ober- | 


flachenverformungen. Ein Anfang in dieser Richtung wurde in der Arbeit von Holstein gemacht; 
eine Zusammenfassung der seitherigen Bemithungen auf diesem Gebiet wird von Ursell? gegeben. 
Das Potential einer pulsierenden Quelle spielt dort die Rolle einer Greenschen Funktion. Die dem 
Verfasser bekanntgewordenen Arbeiten beschranken sich jedoch auf den Sonderfall c = 0.3 


SchlieBlich ist es einleuchtend, daB auch fiir allgemeinere Dispersionsgesetze (also z. B. bei Be- 


riicksichtigung der Oberflachenspannung) die Ausbreitung der Wellen, freilich nicht auch ihre Am- 
plituden, durch die Uberlegungen von Ziff. 3 bestimmt sein wird. In Abb. 7 ist angedeutet, daB 
man bei Beriicksichtigung der Oberflachenspannung nach Mafgabe der dann fiir Wasserwellen giil- 
tigen Dispersionskurve statt der vier Wellenziige wie seither u. U. deren sechs zu erwarten hat. 


(Eingegangen am 7. Juli 1955.) 


Anschrift des Verfassers: Dr. E. Becker, Gottingen, BunsenstraBe 10. 


Vgl. FuBnote 3 von S. 69. 


F. Ursell, Quart. Journ. Mech. and Applied Math., 8, (1954), S. 427. 
Vgl. FuBnote 1 von S. 69. 
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Uber gesteuerte Abkiihl- und Anwarmvorginge bei festen Kérpern 


Von A. Kneschke 


1. Einleitung. Der Abkihlvorgang eines bereits aus der Schmelze erstarrten Kérpers ist von 
echnischer Bedeutung. Er wird in seinem Verlauf von auBen dadurch gesteuert, daB die Tempera- 
ur y(t) der den Kérper umgebenden Luftgrenzschicht zeitlich monoton fallend geadndert wird. Der 
lamit verbundene Warmeentzug bedingt eine instationare Temperaturverteilung im Kérper, deren 
Kenntnis fiir die Durchfiihrung der Kiihlung wesentlich ist. Entsprechend ist der Verlauf des ge- 
steuerten Anwarmvorganges. 

In beiden Fallen wird von einem vollstandigen Temperaturausgleich zwischen dem AuBenraum 
ind dem Kérperinneren ausgegangen, so dafs fiir den Abkithlvorgang dem Korper die konstante 
Anfangstemperatur T,, dem AuSenraum die Temperatur ¢(0) = T,, fiir den Anwarmvorgang dem- 
sntsprechend die Anfangstemperaturen Null und g(0) = 0 zukommen. 

Im folgenden soll fiir beide Warmeleitvorginge die Temperaturverteilung in einem beliebigen 
Kérper, der zur Zeit t = 0 auf konstanter Temperatur ist, unter dem Einflu8 einer zeitlich verander- 
ichen AuBentemperatur g(t) ermittelt und aus der beiden Vorgingen zugrundeliegenden Randwert- 
ufgabe ein allgemein anwendbares Verfahren hergeleitet werden. Dabei wollen wir unsere Be- 
rachtungen an den gesteuerten Abkihlvorgang anschlieBen. Der Anwaérmvorgang wird, wie wir 
sehen werden, gleichzeitig mit erfaBt. 


2. Formulierung der Randwertaufgabe. Im Inneren des Kérpers K vom Volumen V gilt fiir 
lie Temperatur T(P,t) im Punkte P die Fouriersche Warmeleitgleichung 


D{T} = — AT =0, (1) 


vobei a? die Temperaturleitzahl des Kérpers bedeutet. Im Falle des Warmeiiberganges in die den 
Kérper umgebende Luftschicht erhalt die Randbedingung, wenn wir unter / die Warmeleitzahl 
ind unter « die Warmeiibergangszahl verstehen, fiir n als 4uBere Normale der Kérperoberflache 
F die Form 

A OT 
% On 


UT }r=( 4 ay, =p)... (a= 0). (2) 


3. Lésung der Randwertaufgabe. Wir setzen 
2 Pt) TAP) TP.) 
ind fiihren damit die gestellte Aufgabe auf die beiden Ausgleichsprobleme 
DSR iia) wal yee ag) =a | (a) 
Dad (Beye =) (t> 0) | 
ind 
DAT (Pp e0, AT (PO) = 0 | (b) 
LEP Mme olen 0) | 
purtick. 

Die Aufgabe (a) beschreibt einen natiirlichen (ungesteuerten) Abkiihlvorgang, bei dem der Kor- 
er mit der konstanten Anfangstemperatur T', seinen Warmeinhalt an die Umgebung der Tem- 
veratur Null abgibt. Sie wird durch den Ansatz 

TCP y her" ou P) 


uf die Integration von 
Au +0" u =0, (rat) () 
a On FP 


uriickgefiihrt. Bekanntlich besitzt dieses homogene Randwertproblem abzahlbar unendlich viele 
eelle Eigenwerte, die sich im Endlichen nicht haufen. Je zwei zu verschiedenen Eigenwerten 0; 
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gehérige Eigenfunktionen u,(P) sind zueinander orthogonal, und jede zweimal stetig differentiierbare 
Funktion, die den Randbedingungen geniigt, lat sich in eine gleichmaBig konvergente Reihe nach 
den Eigenfunktionen u,(P) entwickeln. Damit wird 


Ty(P,t) Tes) Au (Pye (3) 


V— 


so da sich fiir die normierten Eigenfunktionen die Fourierkoeffizienten 
A, = fu,(P) dv 
k 


ergeben.! : : 
Die Integration von Ausgleichsproblemen von der Art der Aufgabe (b) kann man mittels eines 
Satzes von Duhamel? durchfiihren. Danach wird die Lésung durch ein Integral der Form 
t 
air 
aa | 90) G(P, t— 1) dr 
0 
dargestellt. Fiir unsere Zwecke ist es vorteilhaft, den Duhamelschen Integralansatz durch ein Glied 
y(t) G(P, 0) zu erweitern, so daB wir 


t 


T,(P, t) =4 | 9 G(P,t — x) dr — g(t) G(P, 0) 
6 
setzen. Nach Ausfiihrung der Differentiation und nachfolgender Teilintegration wird dann 


T,(P, t) = (0) G(P, t) + f y' (x) G(P, t— 1) dr — g(t) G(P, 0) . (4) 


Wir erkennen daraus, daB die Funktion ¢(r) G(P, t— 1), physikalisch betrachtet, die Temperatur 
zur Zeit tin einem Punkte P des Kérpers wiedergibt, die unter dem EinfluB der AuBentemperatur 
p(t) im Laufe der Zeit t — 7 entstanden ist. G(P, t) hat also den Charakter einer EinfluBfunktion, 
deren physikalische Bedeutung sich sofort zeigen wird. Der Ansatz (4) befriedigt bereits die An- 
fangshedingung 7,(P,0) = 0. Aus 


D{T,(P, t)} = y(0) D {G(P, 9} +fo'@ D{G(P, t—1)} dr + p(t) a AG(P, 0) = 0 


folgt, wenn D {G(P, 1) = 0 gesetzt wird, die Bedingung AG(P, 0) = 0. Die Randbedingung ist 
entsprechend 


L{T(P, t)}r = 9(0) L{G(P, t)} r + p(t) L{G(P, t—1)}r dt — g(t) L{C(P, hr = 9) 


dann erfiillt, wenn fiir jedes t> 0 


L{G(P,t)}}r =0, G(P,0)=—1 
ist. Wir erhalten damit fiir die EinfluBfunktion G(P, t) die Randwertaufgabe 
DAC P ht == 0... E(P, 0) 


L{G(P, t)}r = 0 (i= Oye 
Sie ist von der Art der Aufgabe (a). Es gilt mithin 


TAP, ) en 
PR Pon og bere 


GP. 1) = 


Die Funktion — G(P, t) ist demnach die Temperaturverteilung, die durch einen natiirlichen Abkiih|- 
vorgang mit der Anfangstemperatur der EKinheit im Kérper gegen den AuSenraum mit der Tem- 
peratur Null hervorgerufen wird. 


* Derartige natirliche Warmeleitvorgange sind behandelt bei H. S. Carslaw, Introduction to the mathe- 
matical theory of the conduction of heat in solids, London 1921; Gréber-Erk-Grigull, Die Grundgesetze der 
Warmeiibertragung, 3.Aufl., Berlin, Gottingen, Heidelberg, 1955; F. Berger, Z. angew. Math. Mech. 11 (1931) 
S.45 und 8 (1928) S. 479. : 


* Vel. R. Courant u. D, Hilbert, Methoden d. mathematischen Physik, Bd. 2, S. 183, Berlin 1937. 
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Auf Grund von (4) lautet die Lésung der urspriinglichen Aufgabe 


t 
(Bast aoe AE ORD yp 
T(P,1) = TAP, 1) + ot) — 910) 8 — | p'(e) 9 ae. 
0 
un Falle der Abkihlung ist hierin (0) = T,, fiir den Anwirmvorgang, der durch die Aufgabe (b) 
viedergegeben wird, ist T,(P, t) = 0 und ¢(0) = 0 zusetzen. Damit wird in beiden Fallen die Tem- 
eraturverteilung im Inneren des Kérpers durch 


: / a P,t— 
T(P, t) = o(t)—| oe) = ar (5) 
0 a 
largestellt. Die Lésung des gesteuerten Abkiih]- bzw. Anwarmvorganges ist also auf die Loésung 
les zugehérigen ungesteuerten (natiirlichen) Abkihlvorganges zuriickgefiihrt. Darin liegt die 
hhysikalische Bedeutung des Ergebnisses (5). 
Setzen wir (3) in (5) ein, so lautet die Lésung 


oO t 
T(P, t) = o(t) —> A, u,(P)- [p'(r) e— 8  —) dr. (6) 
y=1 0 
Jurch Teilintegration erhalt sie die Form 
. Ay Ae 1 1 28 ny Ea: —T 
T(P, 1) = oft) — > AE) | (0) — gO) e-8! — fop"(a) ede |. (6a) 
yv=1 ¥ 0 


Jer WarmefluB ® durch die Kérperoberflache F' ergibt sich aus dem Warmeinhalt 
(oe) t 
Q(t) =ey[ T(P,t) dv =cy Volt) —ey SA | g(t) er —9 dr 
K y=l 9g 


urch Differentiation nach der Zeit t und Beriicksichtigung der Beziehung V = SY A;. Er betragt 


vl 

ic) t 2 

D(t) =A 3) A} 6; [p(t e- 8) dr. (7) 
el 0 
Jurch Teilintegration wird daraus 
foe) t ° i 

P(t) = cy SAF y'() —G'(0) e—8' — | p(x) ede | . (7a) 

y=1 0 


‘emperaturverteilung und WarmefluB in der Form (6a) und (7a) lassen in dem Zeitfaktor 
t 
Z,(t) = gy’ (t) — y’ (0) e— 81 | @"(z) e—276,(t—t) dr 
0 


hre Abhangigkeit von der Geschwindigkeit ~’ und Beschleunigung p” der Kiihlung bzw. Anwar- 
nung erkennen. 


4, Anwendungsméglichkeiten. a) Es sol] die Temperaturverteilung eines unendlich langen 
Creiszylinders vom Halbmesser R und der konstanten Anfangstemperatur T,, bei dem einfachen 
‘all der linearen Abkiihlung 

p(t) = Ta— ot 
nit der Kiihlgeschwindigkeit ~ ermittelt werden. 

Fiir den natiirlichen Abkiihlvorgang, erhalt man aus der Integrationsaufgabe unter Ver- 
vendung der Besselschen Funktionen hekanntlich' die Temperaturverteilung 


2 


: a \ J,(cr) cy aT, t 
EAE aed a 2: 6, (Tke:) + F(e,)) Jo (i r) e ; 
Jabei sind die c,-Werte die Wurzeln der transzendenten Gleichung 
Ie) _ Ae 
J,(c) aR - 


1 Vgl. Gréber-Erk-Grigull, a. a. O. S. 53. 
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Aus (5) ergibt sich sofort die gesuchte Temperaturverteilung zu 


T(r, ji T, espe ne> J,(cr) ae | ( ee 5t 


c3 (Jo ey + J Cy)) 


Die Ubertemperatur 4 T(r, t) = T(r, t) — q(t) einer beliebigen Stelle im Zylinder zur jeweiligen 
Temperatur g(t) des AuBenraumes strebt mit wachsender Zeit t dem Wert 


me p J,(Cv) a 
lim A T(r, t) 2S ce (Je(c y ) + J?(cr)) Jo ( | 


t— oo 


zu. Die hierbei auftretende Reihe 1aBt sich para tone es ist 


LAT (Re es r). 


9 
me 4a 


In der gleichen einfachen Weise lassen sich mittels der Integraldarstellung (5) fur andere Korper- 
formen auch im Falle nichtlinearer Kiihlung bzw. Anwaérmung die Temperaturverteilungen er- 
mitteln. Setzt man g(t) z. B. als periodische Funktion an, so ist man in der Lage, den fiir Re- 
generatoren grundlegenden Wert des Ausnutzungsgrades warmespeichernder Kérper zu bestimmen. 
Uber die Auswertung des hier angegebenen Verfahrens fiir die technischen Belange soll spater 
ausfihrlich berichtet werden. | 

b) Die Lésung in der Form (5) ist auch dann brauchbar, wenn der Verlauf von T,(P, t) nicht | 
analytisch festliegt, sondern fiir einen bestimmten Punkt des Kérpers empirisch durch einen na- | 
tiirlichen Abkiihlversuch ermittelt worden ist. Der Temperaturunterschied | 


t=t : t’=1 
ATP, ) =| pi) | ee 
j a e a 
ist dann durch die von 
' nr PP, ’) 
Y (t t ) os Tr - 


1 


und der t’-Achse begrenzte Flache fiir jeden festen Zeitpunkt t bestimmt. Das angegebene Verfah- 
ren ist mithin auch auf die gesteuerte Abkuhlung oder Anwarmung von Kérpern anwendbar, die 
einer analytischen Behandlung nicht zuganglich sind. 


(Eingegangen am 11. Juli 1955.) 
Anschrift des Verfassers: Professor Dr.-Ing. A. Kneschke, Freiberg/Sa., Richard-Wagner-Str. 13. 
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Hydromechanische Grundlagen fiir die Behandlung von stationaren 
und instationéren Grundwasserstrémungen 


II. Mitteilung 
Von G. Heinrich und K. Desoyer 


In einer vorangegangenen Arbeit! wurde die Theorie der Grundwasserstrémungen einerseits 
if die Behandlung von instationiren Fliissigkeitshewegungen, anderer-eits auf den Fall aniso- 
oper inhomogener Festkérper erweitert. Es sollen nun die in der angegebenen Arbeit durchgefihr- 
n Betrachtungen dahingehend verallgemeinert werden, daB die dort getroffene Voraussetzung 
inktférmiger Berithrung der Kérner des Festkérpers fallen gelassen wird und das Problem auf 
lgemeine Festkérperskelette erweitert wird. 

Wir betrachten fiir das folgende ein zeitlich unveranderliches ruhendes Festkérperskelett, von 
m. wir nur voraussetzen, dafs auf makroskopischen Flachenelementen jedes nicht speziell ge- 
hrten (statistischen) Schnittes das Verhaltnis der von Flissigkeit bedeckten Flache zur gesamten 
ache des Elementes als Funktion des Ortes angegeben werden kann. Dies setzt eine statistische 
erteilung der Festkérpersubstanz voraus. Es ist dabei gleichgiiltig, ob das Festkérperskelett aus 
nem einheitlichen Stoff besteht, oder ob es aus verschiedenen Komponenten zusammengesetzt 
| (z. B. bei Beton aus Fiillmasse und Zementpaste). 

Fiir ein makroskopisches Flachenelement der Hiullflache seij das Verhaltnis des von Fliissigkeit 
deckten Anteils der Flache zur gesamten Flache. Diese GréSe wird fiir verschiedene Arten von 
shnitten (statistische Schnitte oder speziell ausgesuchte Schnitte) verschiedene Werte annehmen. 
ir zeigen zunachst, daB auch fiir anisotrope Festkérper fiir statistische Schnitte J gleich sein muB 
m relativen Porenvolumen n. Dazu betrachten wir ein durch statistische Schnitte heraus- 
trenntes quaderférmiges makroskopisches Volumelement mit den Kantenlangen dx, dy, dz in 
nem beliebig orientierten Koordinatensystem (x, y, z). Die in diesem Quader enthaltene Fliissig - 
itsmenge n dx dy dz muf sich nun auf drei Arten ermitteln lassen: bezeichnen /,,/,,/, die/-Werte 
- die statistischen Schnittflachen senkrecht zu den Kantenrichtungen x, y, z, so muB gelten 


ndx dy dz =/j,dx dy dz =j,dy dx dz =A, dz dx dy. 
ieraus erkennt man, daB, auch fiir anisotrope Festkérper, fiir statistische Schnitte 
Ay =Ay =A, =A SN (1) 


lt (Verallgemeinerung des sog. Delesseschen Gesetzes?). 

Fiir spezielle (nichtstatistische) Schnitte gilt (1) i. a. nicht. Auch solche Schnitte sind jedoch 
aktisch von Bedeutung. So wird man z. B. bei der Beurteilung der Rutschgefahr bei kérnigen 
estkérpern von vornherein nur solche (nichtstatistische) Schnittflachen betrachten, die die Kérner 
-s Festkérpers nicht durchschneiden, sondern durch deren Berithrungsflachen hindurchgehen. 
ir wollen daher im folgenden, um auch diese Méglichkeit zu erfassen, mit einer zunachst offen- 
Jassenen Gréfe A rechnen, die wir auch bei anisotropen Festkérpern als richtungsunabhangig vor- 
issetzen, die aber ortsabhangig sein kann. 

Wir denken uns durch eine beliebige (statistische oder speziell ausgewahlte) geschlossene Schnitt- 
iche einen endlichen Teil des Gemisches (Festkérper und Fliissigkeit) herausgegriffen und wollen, 
nlich wie in der ersten Arbeit, die Krafte, die auf den im Inneren diese. Hiillflache befindlichen 
liissigkeits- bzw. Festkérper wirken, analysieren, um die Bewegungsgleichungen bzw. Gleich- 
wichtsbedingungen aufzustellen. Im Gegensatz zu der Behandlung in der ersten Arbeit erscheint 
fiir die vorliegende Verallgemeinerung erforderlich, die von Flissigkeit bedeckten Teile der aube - 
n Hiillflache mit samtlichen im Innern befindlichen Flissigkeitsgrenzflachen gegen das Festkérper- 
elett zusammenzufassen. Diese Flache umschlieBt dann den ganzen herausgegriffenen Fliissig- 
sitskérper. Bedeutet do, 7 ein von Fliissigkeit bedecktes vektorielles Oberflachenelement der 
iillflache (nach auBen orientiert) und do; ein Oberflachenelement der inneren Grenzflachen der 


1 G. Heinrich u. K. Desoyer, Ing.-Arch. 23 (1955), S. 73. Bei Verweis auf Gleichungen dieser Arbeit wird 
n betreffenden Gleichungsnummern der Hinweis I. vorangestellt. 
2 4. Delesse, Annales des mines (4) 13 (1848), S, 379. 
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Fliissigkeit (gegen den Festkérper orientiert), dann ist die Resultante aller auf diesen Flachen auf 
die Fliissigkeit einwirkenden Normaldriicke gegeben durch 


—Ip da Fi — Ip dp; ’ 

wobei p den Druck in der Fliissigkeit bedeutet. Da diese Grenzoberflache den Fliissigkeitskérper 
vollstandig einhiillt, kann man den Gaufschen Satz anwenden und erhalt 

|p doar + fp 40; = [pp Vn = [pp-ndvV, (2) 
wobei dV, ein Volumelement des Flissigkeitskérpers, dV ein Volumelement des Gemisches be- 
deutet. Die iibrigen Krafte werden ebenso eingeftihrt wie in der ersten Arbeit und zwar iy als die 
von der Hautreibung am Festkérperskelett herriihrende auf die Fliissigkeit wirkende Quasivolum- 
kraft, bezogen auf die Volumeinheit des Gemisches. Wie bereits in der ersten Arbeit gezeigt wurde, 
heben sich die auf den von Flissigkeit bedeckten Teil der Hillflache wirkenden Reibungskrafte bei 
der makroskopischen Mittelwertsbildung auf. Ferner ist VU die Schwerkraft, bezogen auf die 
Masseneinheit. Es gilt dann 


[(-" Vp + Rw — nVU)av =0, (3) 


wo yw das spezifische Gewicht der Fliissigkeit und g die Fallbeschleunigung bedeutet. 
Da das obige Integral fiir jedes beliebige Volumen verschwinden mu, so erhalt man durch 
Nullsetzen des Integranden 


ity = (Vp + VU). (4) 


Da (4) mit (I. 1) tbereinstimmt, spielt also die Voraussetzung der Punktberiihrung der Kérner fiir 
die Bewegung der Fliissigkeit titberhaupt keine Rolle. Es bleiben daher alle Gleichungen der ersten 
Arbeit, die sich auf die Bewegung der Fliissigkeit beziehen, auch fiir den verallgemeinerten Fall 
unverandert. \ 

Um die Gleichgewichtsbedingungen fiir den herausgeschnittenen Anteil des Festkérperskeletts 
aufzustellen, denken wir uns, ahnlich wie in der vorangegangenen Arbeit, die Schnittspannungen, 
die im Festkérperanteil der 4auBeren Hillflache wirken, in der makroskopischen Betrachtung iiber 
das ganze makroskopische Hiillflachenelement gemittelt. Derin einem solchen Element herrschende 
Spannungszustand kann durch Einfiihrung eines symmetrischen Spannungstensors © dargestellt 
werden. Bei kérniger Festkorperstruktur kénnen zur Beurteilung der Rutschgefahr spezielle (nicht- 
statistische) Schnitte so gelegt werden, daf keines der Korner geschnitten wird. In diesem Fall 
stellt 6 den Spannungstensor der Korn-zu-Korn-Krafte dar. Bei statistischen Schnitten (hier wire | 
das Festkérpermaterial durchschnitten) stellt © den statistisch gemittelten Spannungstensor im 
Festkérpermaterial dar. Dieser statistische Mittelwert des Spannungstensors ist iiberall dort maB- 
gebend, wo es auf das elastische oder plastische Verhalten des Festkérpers ankommt. Dies wird 
spater noch begriindet. 

Die Gleichgewichtsbedingung fiir den herausgeschnittenen Teil des Festkérperskeletts lautet 


[S- do. — [aw av + [pao— [A —m) PU dV 208 @ 


Hier bedeutet do, den nach aufen gerichteten Vektor eines Hiillflachenelementes. Das zweite und 
dritte Integral sind die Reaktionen der betreffenden auf den Fliissigkeitskérper wirkenden Krifte 
[vgl. (3)]. Das letzte Glied stellt die auf den Festkérperanteil wirkeude Massenkraft (Schwerkraft) 
dar. yx ist das spezifische Gewicht des Festkérpermaterials. 
Aus der Bedeutung der oben eingefiihrten GréBe J ergibt sich der Zusammenhang 

dq jap = A dd, : (6) 

Aus (2) folgt dann 7 
[p do, = [Vp -ndV — [pA do, 


und daraus mit Hilfe des GauBschen Satzes 
[p do; = [Vp —V (Ap) av. (7) 


Durch Einsetzen von (7) in (5) ergibt sich nach Anwendung des Gaufschen Satzes auf das erste | 
Integral von (5) . 


[|v -E— Rw — (L—n) “VU + (n—A)Vp— pVAlav =o. (8) | 
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Ja dieses Integral fiir jedes herausgegriffene Volumen verschwinden mu8 erhalt man durch Null- 
etzen des Integranden und Verwendung von (4) 


V-G=AVp4 *e— nee) 7 4 pV A. (9) 


’ + © stellt die auf die Raymeinheit bezogene Resultante der Schnittkrafte im Festkérperskelett 
lar. 

Bei kérniger Struktur des Festkérpermaterials ist fiir punktformige Beriihrung der Korner 
= 1 zu setzen, wenn Schnitte durch die Beriihrungspunkte der Kérner betrachtet werden. Dann 
reht (9) in die Gleichung (I. 3) tiber: 


VG =Vp + "eh py, (9a) 


fiir statistische Schnitte (A = n) nimmt (9) die Form an 


V-G=nDVp-4 Kane) CU Vr. (9b) 


Damit ist eine alte Streitfrage aufgeklart, die in der hydraulischen Literatur bis heute noch nicht 
sereinigt werden konnte, die Frage namlich, ob bei der Wirkung der strémenden Fliissigkeit auf den 
Festkérper der Druckgradient in voller Starke oder mit dem Faktor n versehen einzufiihren ist!, 
Es ist dies also lediglich davon abhangig, ob der betrachtete Schnitt ARSE gefiihrt oder durch 
lie Beriithrungszonen der Korner hindurch gelegt wird. 

Analog wie in der vorangegangenen Arbeit kann man den Druck p in (9) durch die Standrohr- 
spiegelhéhe h ausdriicken gemaB (I. 11) 


U P 
7 + ms i) (10) 
und erhalt 
A—n ene U 
mit (I. 12) 
Vh =—-1-y. (12) 


Aus (11) und (12) erkennt man, da die auf die Volumeinheit des Gemisches bezogene resultierende 
Schnittkraft im Festkérperskelett sich aus drei Anteilen zusammensetzt: der erste Anteil rihrt von 
der Strémung der Fliissigkeit her. Der zweite Anteil entspricht dem Gewicht des Festkérper- 
skeletts, das um einen i. a. reduzierten Auftrieb vermindert ist. Im Falle der Punktberiihrung ist 
fiir ausgewahlte Schnitte durch die Bertihrpunkte (A = 1) der volle Auftrieb wirksam. In diesem 
Fall stimmt (11) mit (I. 48), bei deren Ableitung diese vorausgesetzt wurde, tiberein. Fiir statisti- 
sche Schnitte (A = n) verschwindet der Auftrieb zur Ganze, unabhangig davon, ob Punktberiihrung 
vorliegt oder nicht. Fiir ausgewahlte Schnitte mit 7 < n wechselt der Auftrieb das Vorzeichen. 

In allen Fallen also, in denen es auf statistische Mittelwerte der Spannungen im Festkérper- 
skelett ankommt (z.B. wenn man Spannungs—Dehnungs-Beziehungen verwendet, um das elastische 
oder plastische Verhalten des Festkorpers zu untersuchen), darf, wie (11) zeigt, nicht mit einem 
Auftrieb gerechnet werden, da in diesem Fall der Festkérper statistisch zu schneiden ist (A = n). 

Letzteres ergibt sich aus der folgenden Uberlegung: In einer Kontinuumstheorie kénnen die 
Deformationen bzw. Deformationsgeschwindigkeiten nur einen Zusammenhang mit den Schnitt- 
spannungen besitzen, sie kénnen aber nicht direkt in Zusammenhang gebracht werden mit zusatz- 
lichen Kraften wie z. B. Druckkraften, die von der Flissigkeit auf das Skelett wirken. AuBerdem 
mu verlangt werden, dafs der Zusammerihang zwischen den Deformationen bzw. Deformations- 
geschwindigkeiten und den Schnittspannungen im trockenen und fliissigkeitserfiillten Skelett der- 
selbe bleibt, um die im trockenen Zustand ermittelten Materialkonstanten in beiden Fallen ver- 
wenden zu kénnen. Man mu daher die Schnitte so fihren, dab im fliissigkeitserfiillten Skelett 
keine zusatzlichen die Schnittspannungen andernden Auftriebskrafte auftreten. Diese Bedingung 
arfiillt nur der statistische Schnitt. 


1 Veil. K. Terzaghi, Erdbaumechanik, S. 128, Leipzig und Wien 1925, wo der Druckgradient in voller 
Starke eingefiihrt wird, und E. Mayer- Peter, H. Favre, R.Miiller, Schweiz. Bauz. 108 (1936) S.37, wo der Druck- 
gradient mit dem Faktor n versehen ist. 
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Es sind daher alle in der Literatur vorhandenen Ansatze, die in diesen Fallen ein Auftriebsglied 
enthalten, fehl am Platze.1 

Der dritte Anteil in (11) tritt nur bei ortsabhangigem / auf. Praktisch von Bedeutung ist vor 
allem der Fall, da® man fiir zwei langs einer Flache aneinander grenzende Festkérperskelette mit 
zwei verschiedenen A-Werten (A, und A,) langs dieser Flache zu rech- 

nen hat. : 
Es sei ¢,, der Einheitsvektor senkrecht zur Flache, vom Kérper | in 
” den Kérper 2 weisend. Wir ersetzen zunachst die Unstetigkeit von / 
durch einen stetigen Ubergang 2 =A(v) im Intervall y = —¢, bis 
y =e, und gehen dann zur Grenze ¢, > 0, €,—> 0 iiber. Wir errichten 
iiber einem Flachenelement dF, der Grenzflache einen kleinen Zylinder 
mit der Hohe ¢, + ¢, (s. Abb. 1). Bildet man nun das Volumintegral 


von Ay (r — al VA [drittes Glied von (11)] tber diesen Zylinder und 


fiihrt den angegebenen Grenzitibergang aus, so erhalt man bei Be- 
Abbal. achtung von (10) 


lim | Yn (h— Ad Fs Fhemeibe i pVA(v) dF,.dv =e,,4F, lim | po dv =eydFyyp(y—7a)- (13) 


6,70 
€.70 »y=—¢, 


Ist also die Resultante aller Schnittkrafte, die langs einer geschlossenen Hullflache auf das Fest- 
kérperskelett wirken, zu bestimmen, so gilt, wenn im Inneren Unstetigkeitsflachen obiger Art vor- — 
kommen, gemaB (11) und (13) 


i; 


S [0-8 aV,= 3 [yyw Ph + (yj — FM yw) = vlay, 
J A] 6 
jth 


ay 
5) 


i 
+3 [wp Visav, + & | oy Pde Say) aa (14) 
: V; jel Fi] 
Dabei sind die Volumintegrale iiber die durch die Unstetigkeitsflachen F; voneinander abgegrenz- 
ten Einzelvolumina V; zu nehmen. Das letzte Integral ist tiber alle Unstetigkeitsflachen zu bilden. 


Bei ortsabhangigem/ beeinfluBt also der absolute Druck p in der Fliissigkeit die inneren Spannungen 
im Festkérperskelett. 


(Eingegangen am 11. Juli 1955.) | 


Anschrift der Verfasser: Prof. Dr. G. Heinrich und Dr. K. Desoyer, Wien IV, Karlsplatz 13, 
Technische Hochschule. 


1 Siehe z. B. I’. Tolke, Ing.-Arch. 2 (1931), S. 295ff. Vgl. auch die Polemik zwischen P. Fillunger und 


O. Hoffmann in ,,Die Wasserwirtschaft (1929) Heft 18, 20, 21, 33 : : ae 
Riteueeaen. (1929) » 33, (1930) Heft 4, 25, die keine endgiiltige : 
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Die Knickformeln fiir den Stab mit n Teilstiicken 
konstanter Biegesteifigkeit 


Von S. Falk 


__ 1. Einleitung. Knickformeln fiir den abgesetzten Stab mit nur zwei Teilstiicken sind in der 
technischen Literatur seit langem bekannt; fiir drei und mehr Teilstiicke jedoch findet man héch- 
pstens noch stark symmetrische Stabe durchgerechnet!. In der vorliegenden Arbeit werden nun die 
Knickformeln fiir den Stab mit beliebig vielen Feldern hergeleitet und an zwei einfachen Beispielen 
ivorgefiihrt. 


| 2. Der gewohnliche Knickstab. Fiir den gewohnlichen Knickstab mit konstanter Lings- 
tsraft P und konstanter Biegesteifigkeit EJ (Abb. 1) gilt die Differentialgleichung 


ya) + y"(x) =0 mit PF (0); (1) 


{mit den vier linear unabhangigen Lésungen 1, x, cosy x und sin y x, Wahlt man als Integrations- 
ysonstanten Durchbiegung y,, Neigung yj, Biegemoment M, und Querkraft Q, am linken Balken- 
onde (x = x9), so lautet die Lésung von (1) mitsamt den ersten drei Ableitungen 


cos ¥ (x—x,)) — 1 ‘sin » (x — xX) —» (% — x) 


ya) = ¥9 + (# — ¥) 70+ 5 M, + - Qo 


EA at ne eau — = Sa) M, ae cos v(x ae Xo) — 1 On A 
eH 2 M(x) = cosy (x—%,) M,-+ ~ sin y (x — x) Qo, 
— EK y’ =QA(a) = —ysin y (x — x) M, + cos y (x — %) Qo 


+der, indem wir die Koeffizienten dieses linearen Gleichungssystems zur ,,Feldmatrix” 7}(x) und 
ie GréBen y, y’, M und Q zum ,,Vektor“ t) zusammen- 
yeassen, 
(x) = B(x) Yo- (3) 
ssetzen wir x = x, ein, so geht die Feldmatrix in die Einheits- 
aatrix © iiber, wie man in (2) leicht nachpriift. Es ist also 
(x) = ©) = Yo, wie es sein muB. Dagegen wird fiir x = x, Abb Gewahnlicher Knickstab mit konstante 
Dy = %(ac,) Do = Q 5 (4) iegefestigkeit EJ und konstanter Lingskraft P. 
er ist L == %(x,) die sogenannte ,,Leitmatrix’’, die nach (4) den linearen Zusammenhang von 


gpurchbiegung, Neigung, Biegemoment und Querkraft an beiden Enden des Stabes vermittelt; 
‘aan erhalt sie, indem man in dem Koeffizientenschema des Systems (2) x = x, einsetzt; und zwar 


LS 


( 


t 


a 


hird wegen x, — x, =/ 


if l (cosy 1—1)/P nia ane Wg 

0 1 —ysinvyl/P (cos y 1 — 1)/P 
: v= (%) =| 0 0 cos vl = sin vl |. (5) 
0 0 —ysiny 1 cos y 1 } 


| 
i. nach Art der Lagerung des Stabes treten zur Differentialgleichung (1) noch vier Randbedin- 
sungen; sind diese linear und homogen in den Integrationskonstanten, so hat (1) nur dann eine 
jon y(x) = 0 verschiedene Lésung, d. h. der Stab kann nur dann ausknicken, wenn die Deter- 
sainante D dieses Gleichungssystems verschwindet; somit ist D = 0 eine Bestimmungsgleichung 
\ ir die Kraft P, und die in der Praxis meist allein interessierende kleinste positive Wurzel dieser 
‘jranszendenten Gleichung heift die Knicklast P,,. 
‘ 1 Eine groBe Anzahl solcher Formeln findet man zusammengestellt bei A. Pfliiger, Stabilitatsprobleme 
Jer Elastostatik. Berlin, Gottingen, Heidelberg 1950, S. 249—267. 
fi 
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Besonders einfach sind nun solche Stiitzungen des Stabes, die an jedem Ende das Verschwinden 
zweier der GriBen y, y’, M, Q verlangen; in der Leitmatrix @ sind dann offenbar lediglich zwei 
Spalten und zwei Zeilen zu streichen, und die verbleibende Unterdeterminante zweiter Ordnung 
ergibt, gleich Null gesetzt, bereits die gesuchte Knickbedingung fiir P. Die Tabelle 1 zeigt sechs 
einfache solcher Lagerungen mit den zugehorigen Randbedingungen; die eingeklammerten Be- 
dingungen am rechten Ende in b) und e) sind bereits durch die bequemeren dritten Bedingungen 
links mit erfillt, somit iiberfliissig, wovon man sich leicht iiberzeugt. Die Nummern oberhalb der 
Stabenden geben an, welche Spalten und Zeilen in der Leitmatrix {2 (5) zu streichen sind, Die gleich 
Null gesetzten verbleibenden zwei- bzw. einreihigen Unterdeterminanten findet man in der zweiten 
Spalte der Tabelle zusammengestellt; die ersten vier von ihnen sind die bekannten Eulerschen 
Knickbedingungen. 


Tabelle 1. 
Art der Stiitzung mit Rand- " ‘ ick- 
bedingungen. Nummern der in @ Knickformeln fiir tis ee ai as oe 
zu streichenden k Cis non 
Spalten Zeilen ein Feld n Felder (alle 9; # 0) 
i ee | ee 
fab ff zZ Ds elf 
3 4 : l ~ 2 
, A 
y=0 y=0 ™n } 
M=0 M=0 Z 
bey 1 
a & 
A Hs — 
cos vl = 0 CaO) : 
y=0 M=0 Pe > P(e) 
‘ f a 
y=0 (Q=Py) 7 
Y=0 
Cae 3 ; (he Gna 
4 4 2—2cosvl=yvlsinvl (Cn 1), 1) 


BE a: h. Se S, =», US, 


y-0 y-0 Pass ec, | cdc wegen (13) 
y-0 ¥-0 2 yaar da) Beane Pk 
vA 
d i é 
Se too —— ol 20 =a | or ee 
y-0 y=0 n 
y=0 M=0 


COST —1() Co 


A, >P (zy 


§ 

3. Der Knickstab mit n Teilstiicken, Nun sei ein Stab mit abschnittsweise verschiedener, aber 
konstanter Biegesteifigkeit EJ; nach Abb. 2 gegeben. In jedem Teilstiick gilt natiirlich nach wie 
vor die Differentialgleichung (1), jedoch mit verschiedenen Parametern »;. Feld- und Leitmatrix 
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haben daher immer noch die Form (2) baw. (5), nur ist jetzt y durch »; und x — x, durch x — x; 
bzw. / durch]; zu ersetzen. Da nun an jeder Feldgrenze Durchbiegung, Neigung, Biegemoment 
und Querkraft stetig ineinander iibergehen, ist einfach 


Y, = 21 Yo> Yo = Vol » 


Dn = AS Yn—1 ? (6) 
also 
Yn ae ue eal see ae a Yo (7) 
er oe ee (8) y 


i - 5 ; Abb, 2. Mehrteiliger Knickstab; Lingskraft P im gesamten Stab 
Die Gleichung (8) aber hat dieselbe Form wie konstant, Biegesteifigkeit EJ; stiickweise raat 
(4), und das hei®t, daB in der Produktmatrix 

%, die gleichen Spalten und Zeilen zu streichen sind wie in der Leitmatrix { des gewdhnlichen 


Knickstabes, womit zumindest formal das Problem des n-teiligen Stabes auf diesen zuriick- 
gefiihrt ist. 


Wir multiplizieren jetzt nach (7) die ersten beiden Leitmatrizen miteinander und bekommen 
pmit den Abkiirzungen cos 7,1, = c,, siny, 1, = s,, cos v9 1, = cg und sin y, 1, = sy 


¢—1 Sah [ co— 1 Sy — Vel, 
; a l, 3 Pa | | A est 
eed Vat! qt eo) Pe! 
| a= | 0 : P P j = | 0 1 Pe Spe 
| e Q Cs $1; 0 0 Cy S9/Vo 
| 0 0 ——¥, 5; Ci ) 0 0 == Sy Cy } 
yp. Vv 
[ (« Bs S| Sp 1 ce, s, + — s,¢,) —», (ly + Ie) \ { (9) 
) eng Re "2 3 ee es att 
| : ip Py; | 
v yp 
—, (« Se a Sy | (« = 7 S\ s] —] 
Qe Que ] ee ee ! Ph oe les 
w2 1 P {2 ¥ 
Vs, 
| 0. 0 (c Cg —— 8} s] (c Sy + mai ea) | 
vy Vo 
(0 0 —~Vpy\C, Sp + — S$; Cy €, Co — — $4 Se ] 
2 Vo V; 


‘der gewohnlichen Funktionen sin und cos gewisse allgemeinere Funktionen auf, fiir die wir ge- 


q erkennt, daB die Produktmatrix , <’, die gleiche Form hat wie 0, und Q,, nur treten statt 
=: Bezeichnungen einfiihren wollen. Mit den Abkiirzungen 


(12) 


I : 7 © “A Vv; a {. . 
Cc; = COs p; l; Py S$; = 81Nn V; l; 9 l; == l; + one l, -|- l, “ ee; = ar . Ee; = 7 (10) 
vay U 
\ schreiben wir 
4 4 “” 
i C; = Cy—1 Cj; — @i—1 Si—1 Si > So = Oh Oy a Cal Se | 
2 4 > 2, a is (11) 
} Ce CC eee eS | Sins AS ra oes | oe = ST St Oe | 
Inshesondere sel 
| 65. =e y= Aas ee ee | 
f} dann ist 


by 10, = = C08 y, hy, S, = Ses siz, | 


| AuBerdem besteht ein leicht zu bestitigender Zusammenhang, der an den Satz des Pythagoras 


a . 
)- erimnert: 


| C; C; -{- S; S; = ted Cay + Si S;—1 pray OS Sa Cs Cy | $y Sy = C9 C9 +. So Sy = i ; (13) 


«Da nun, wie ein Blick auf (9) zeigt, auch bei allen weiteren Multiplikationen die Form der Produkt- 
« matrix erhalten bleibt, 1a8t sich mit den durch (11) eingefiihrten Funktionen direkt die in (8) be- 


| 7 
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notigte Matrix ‘f,, angeben: 


{ ¢, —1 Sn — nl 
Te ait 3 a 
a Ss Cd pcm! my 1 
m,—|0 1 “gt oo ( = 31). (14) 
i=1 
0 +0 Ch $n/Pn 
0 0 Un Sn i ) 


Sie stimmt duGerlich bis auf sie Zeichen ~ und ~ genau mit der Leitmatrix eines gewohnlichen 
Knickstabes iiberein, und damit sind auch schon die in der dritten Spalte der Tabelle 1 angegebenen 
verallgemeinerten Knickbedingungen gefunden. Die Faktoren », in e) und 1/», in a) sind dabei 


nicht etwa iiberfliissig, da ja der n-te Abschnitt des Balkens durchaus starr sein kann (y, = /P/EJ n=0 
wegen EJ, = 00), worauf man gegebenenfalls achten muf. 

Interessiert schlieBlich auch die Biegelinie des Stabes selbst, so setzt man die berechnete Knick- 
last in die Leitmatrizen @; ein, berechnet nach (6) die Vektoren {,, ),... 1), und hat damit die 
Biegelinie im Feldei als erste Komponente des Vektors 


Hi(x) = Gila) Yio - | (15) 
4. Die explizite Darstellung der Knickformeln. Um nun die in den Knickformeln auftretenden 
Ausdriicke (11) in eine handliche Form zu bringen, kann man zwei Wege einschlagen: 
a) Wir dividieren die Formeln (11) nach jedem Rekursionsschritt durch den Faktor c; = cos »; 1;, 
dann entsteht mit den Abkiirzungen 


~x Ci ~ C; a 3; ~ x 3; } 
Pi = G g++ + Cis cj =, ¢; =—, si. —">=" 5 sj; =, tel (16) 


das leichter zu ttbersehende Formelsystem 


C= C1 — Ct G Si Sj; = ¢j_1t; + @j_1 Si_i1, | 
Sak 2 Oe isd ~ * ~* ae ind ~* (17 
CR == OE en UR OF ig Se SS Cp a Ge SHE o J ) 


Die ersten drei Funktionen c? und s} sind zum Beispiel 


cf =1 Ce, | 

hs “- sn Sx “~ 

cg =1l—e,t,t,, 82 = 1, 4+ et, (17a) 
~* “an “~ n~ “~ vas SAS n 

c3 = 1— eé, tot, — eg t, (tg + e,4,), S3 = (1 — @, tg t,) ty + @ (to + @, ty) | 


usw., ¢} und s} genauso, wenn man e; durch e; ersetzt. Oder noch deutlicher mit den gegebenen 
Werten »; und den ausfiihrlich hingeschriebenen Produkten p; aus (16) 


“& 
Q x 
te Vy Sy 
pe ae i ——$. is Vy t, 5 
Cy cy 
c t; Vous 
a Wats ae 223 05 
142, +2 — yt, + yt, 
Cues Vo Cie, {shania (17b) 
Co t, t. Ve $ 
=U | sats 2 __(y i yt LEE CR ty 
ea Vo bs Se SSH Vo t > t. — = ss 
cee ¥, 14 ate) os mela 14 + 2 tg + Ys ty Mity Ya fs 
usw. 
toe ee il SH t 
Gi : aia eS 
Vy Cy UA 
Se — us Yo to S» aes by by 
ae, yy Pa ba ire iaea tea (17¢) 
26) Cy Vy Vo 
C5 by ty t s t 
are aie, rate—(; {2 Deis 3 Sets ty re tg ue fe ts 
1&2 C3 1 1 Ys Vg Cy Cy Cy Vy Vy V3 2? Vy, 
usw. 


b) Mit Hilfe der Beziehungen 


2 (cos w cos 6 — k sino sin B) = (1 + k) cos (x + B) + (1—k) cos (« — f) , | 
2 (cos « sin B + k sin x cos f) = (1 + k) sin (x + 8) + (l1—k) sin(x—f) | (184 
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gehen die Formeln (11) iiber in 


poe) ey +(L—@4) ea, 2s; = (1 + e_1) ee @— 41) $21, | 
¢ @ ~ ~ A 19 
Be (1 - e;_4) ee gie( Ee €55)) 24, 2s; = (1 + e;_)) go + (1 — e;_1) Sane z 
wobei die Symbole © und © bedeuten, daf simtliche in den Ausdriicken CN, Clits ap und: s: 
vorkommenden Argumente 6 + Si— coal 


im Falle des Zeichens zu ersetzen sind. 


{® durch »; 1; + 6 | 
|© durch y; 1; — 6 J 
(Fir die ersten vier c; zum Beispiel erhalten wir mit der voriibergehenden Abkiirzung 1; 1; = i 
iL c, E—=COS)(1'), 

2 ey = (1 + .¢,)-cos.(2 + 1) + (I—¢,) cos (2 — As 

Ac, = (1 + @,) [(1 + @,) cos (3 + 2 + 1) + (1— e,) cos (3 + 2— 1)] 

+ (1 — e,) [(1 + e,) cos (3 — 2 — 1) + (1 — @,) cos (3 — 2 + Ls 

e,=— (1 + es) {d - ey) [(1 + e;) cos (4 + 3 +2 -+-1) + (1 — @,) cos (4 + 3 + 2— 1)] 


(20) 
+ (1 — eg) [(1 + @,) cos (4 + 3 —2— 1) + (1—@,) cos (4 + 3—24 DJ}, 
( 


| + (1— @,) {(1 + es) [(1 + @) cos (4 — 3 —2— 1) + (1— @,) cos (4 —3—2 + 1)] 


; (1 — e,) [(1 + @,) cos (4—3 42+41)+(1 e,) cos (4— 3 +2—1])]}, 
vusw., S; bzw. ¢;, 8; genauso, wenn man die Zeichen cos durch sin, bzw. @; durch é; ersetzt. 

te) als reine Summen von cos- bzw. sin-Funktionen im allgemeinen praktischer sein, zumal sie 
‘sich durch ein umlaufendes Vektordiagramm leicht veranschaulichen lassen. Die Anzahl der Sum- 


alle @; <1 ae & <1 
alle 1-@; >0 : alle 1-8; >0 


| Da im (17) lauter Produkte von trigonometrischen Funktionen auftreten, werden die Ausdriicke 
i 


1 a G6 
| | fe El e= Oi | EG, <E<E-<Ehy 
! 


y By = Dy y U> yr > Wy 


Abb. 3. Einseitig verjiingte Stabe. 


jmanden betragt in (17) wie in (19) 2"—", allerdings geht sie bei symmetrischen (nicht unbedingt 
‘symmetrisch gelagerten) Staben wegen », 1, =, 1;. Yn—1ln—1 = V2 ly USW. auf die Halfte zuriick. 
‘Bei zu groBem n wird man mehrere benachbarte Abschnitte einmal durch den gréften und ein 
Jandermal durch den kleinsten Wert der vorkommenden Biegesteifigkeiten ersetzen und so mit 
‘Hilfe der exakten Lésungen fiir diesen auf wenige Teilstiicke reduzierten Stab die wirkliche Knick- 
‘last in Schranken einschlieBen. 


_ 5, Untere Schranken fiir die kleinste Knicklast. Insbesondere betrachten wir jetzt solche Stabe, 
idie sich nach rechts (bzw. links) verjiingen (Abb. 3). Da dann offenbar alle Zahlen | — e; (bzw. 
1 — é;) positiv ausfallen — die 1 + é; und 1 + @é; sind es ja ohnehin — werden sadmtliche Vorfak- 
ttoren in den Funktionen ¢;. $; (bzw. ¢;, §;) nach (19) ebenfalls positiv, und daraus folgt, dafs die 
‘kleinste positive Nullstelle der Summenfunktion weiter rechts liegen mu als diejenige des Summan- 
iden mit dem maximalen Argument 0 = 7,1, +++-¥2l, +74. Damit ergeben sich dann die 
‘in der letzten Spalte der Tabelle 1 angegebenen unteren Schranken fiir die gesuchte kleinste Knick- 
last P,. Aus Symmetriegriinden folgt noch, daB die Abschatzung in den Fallen a) und e) auch dann 
‘gilt, wenn der Stab nach rechts bzw. links verjiingt ist. 


6. Stabe mit starren Zwischenstiicken. Falls einige Abschnitte des Stabes starr, ihre Biege- 
‘steifigkeiten also unendlich groB sind, werden die zugehérigen »; = 0. Man schreibt dann in den 
(Forme!n (17) bzw. (19) wieder ¢; = v;/v41. @ = vi4i/vi und fiihrt mit allen Werten »;. die gleich 
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Null zu setzen sind, soweit sie im Nenner vorkommen, die Grenziibergange 
sin »; l; . ter; 1; 
le lm et 
Vj 


lim 


v.>0 Vi v0 
i 7 


(21) 


durch, was sich in den Formeln (17) besonders bequem gestaltet, wahrend man in (19) erst die 


Additionstheoreme heranziehen mu. Die zu»; = 0 gehérige Leitmatrix wird nach (5) 
1 l; 0) 0 
1 0 0 
Ne = } 0 il i fiir y= 0 5 (22) 


0 0 0 1 
und das heiBt, da der Stab in diesem Teilstiick unverbogen bleibt. 
7. Beispiele. a) Erstes Beispiel: Man berechne die Knicklast des Stabes der Abb. 4 mit 
EJ, =9EJ, EJ, =4EJ, EJ, =EJ,1, =1/4, ly =1/4, 1, =1/2. Nach Tabelle 1 liegt der Knickfall b) vor; 


die Bedingung heiBt somit cs =0. Mit der Abkiirzung y= V P/IEJ stellen wir zunachst samtliche 
benétigten GréBen zusammen: 


y vl y y 2 ox 5 in l 
WET 8 Vyl =T5> Oh ge ee Le epre 3 a aes 
vil a Vp 1 A 3 ne edl 
(Raia ge) Val==> bi ameter eral Fi De Oh = sys eed eta 
| vi 
1S arate Use eo a 


Da sich der Stab nach rechts verjiingt, ist nach Spalte4 der Tabellel mit @ = y, 1, + v2 1, + 51, 


= 17 yl/24 
p St anaes 
Poses, =P l7vl 
WAIN I oe 
=P (Soe 2) ae 


Abb. 4. Knickstab eingespannt-frei, drei Teilstiicke. 


Der genaue Wert folgt aus (20), dritte Zeile: 
a 1 i: Tvl 1 ral 


oN Gates On ce seo 


cos 


fs I7vl 1 
3 24 


| 
oy eae! 


oder einfach mit » 1/24 = z 
24¢, = 0 = 15 cos 172+ 3cos 1324+ 5cos7z-+ cosllz. 


Den Verlauf der vier Summanden und der Funktion 24 ¢, zeigt Abb.5; man findet als erste Null- 
stelle z = y 1/24 =0,11 und damit P, = 6,98 EJ/I?. 

Zweites Beispiel: Man berechne die Knicklast des Stabes der Abb.6 mit EJ, =co, EJz=EJ, | 
EJ; =oo, |, =l, =2 1/5, 1, =1/5; Knickfall a), Bedingung $,/v, =0. Hier ist nun ¥, =%3=0, also 7] 
gibt die dritte Zeile von (17c) nach (21) und wegen cos0 = 1 ; 


S3 -=h4+2+h—hytl,, 
2 


Pou leuCOSioll ace 
so OF cha Is C nee G ee ) tg M 
2 


Fir die angegebenen Werte findet man hieraus mit »,l, = z 


~ ada +h) ed 3% on 
StS nt = eo 


und entnimmt der Abb.7 die kleinste Wurzel der Gleichung tg z=f (z) als z=y,1,=1,51. Also ist 
5 21y\2 Peale . 
(1 1)? = 5" = se p7 = 151 =2,28; somit P, = 2,28 ( aud — 14,95 EJ/-. 


41? 


8. Knickstibe mit elastischen und starren Zwischenstiitzen und sprunghaft veranderlicher Langs- 
kraft. Wir nehmen zunichst an, der Stab sei an beliebig vielen Zwischenstellen durch Federn und _ 


7 
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Drehfedern elastisch gestiitat; auch mége sich an jeder Feldgrenze die Langskraft sprunghaft 
andern. Die Differentialgleichung (1) gilt dann immer noch, auch die Feld- und Leitmatrizen 


werden davon nicht beriihrt. Da jedoch in den Feldrindern nicht mehr alle vier GréBen y, y’, M 


und Q stetig ineinander iibergehen, hat man gewisse von den Federkonstanten und Langskraften 
abhangige Ubergangsmatrizen l|; einzufiihren, und die Produktmatrix 8, lautet jetzt allgemeiner 


Sn = eee Rn Us os Us oe ue . (23) 


Oe ccc; bree! : Ree ; 
‘ 2 er Matrix sind wiederum zwei Zeilen und zwei Spalten zu streichen; man braucht daher 
auch von vornherein nur je zwei Spalten gleicher Nummer der 2n +1 Matrizen der rechten Seite 


Abb. 5. Zeichnerische Lésung fiir den Stab von Abb. 4. Abb. 7. Zeichnerische Lésung fiir den Stab von Abb. 6. 


von (23) zu berechnen. Fiir Ul; =€ (i=1, 2,...,n-+1) geht (23) wieder in (7), (8) tiber. Bei 
starren Zwischenstiitzen ist die Methode etwas abzuandern. Naheres dariiber findet man bei 
W. Schnell! und in einer Arbeit des Verfassers?. 

Zwar lassen sich auch jetzt noch allgemeine Knickformeln angeben, doch werden diese schon 
bei drei Feldern so umfangreich, da} man in der Praxis besser anders vorgeht: ein fester Wert P 
wird in die Leit- und Ubergangsmatrizen eingesetzt, nach (23) die Produktmatrix ‘f,, und aus dieser 
der Wert einer gewissen durch die Randbedingungen des rechten Endes festgelegten zwei- 
reihigen Unterdeterminante berechnet und iiber P aufgetragen. Durch Interpolation bestimmt 
man dann den Nulldurchgang, ein Verfahren, das fiir Schwingungen von Kurbelwellen in der 
Ingenieurmechanik unter dem Namen von Holzer-Tolle seit langem bekannt ist. 


(Eingegangen am 14. Juli 1955.) 
Anschrift des Verfassers: Dr.-Ing. Sigurd Falk, Braunschweig, Celler Strabe 93a. 


1 W. Schnell, Z. angew. Math. Mech. 35 (1955) Nr. 6/7, S. 269. 
2 §. Falk, Abhandlungen der Braunschweigischen Wissenschaftlichen Gesellschaft 7 (1955), S. 60. 
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Theorie der Hochdruckréhrenfeder 
Von W. Wuest 


1. Grundgleichungen. Hochdruckréhrenfedern nach Art des Bourdonmanometers werden 
zum. Messen von Driicken im Bereich von etwa 100 bis 10000 kg/cm? benutzt. Das MeBelement 
besteht dabei aus einem kreisférmig gekritmmten Rohr von ovaler Querschnittsform (Abb. 1), dessen 
eines Ende verschlossen ist, wahrend das andere Ende fest mit der Druckzufiihrung verbunden 
ist. Bei Einleiten eines Druckes sucht sich der ovale Querschnitt der Kreisform anzunahern. Wegen 
der Kriimmung der Feder werden dabei ,, Ringspannungen“ hervorgerufen, welche den Kriim- 
mungshalbmesser zu vergroBern trachten. 

Es ist bekannt, daB diese Ringspannungen bei diinnwandigen Federn die elastischen Defor- 
mationen des Querschnitts wesentlich beeinflussen. Man kann also die Formanderungen von — 
Niederdruckfedern nicht berechnen, ohne yon vornherein die — 
Riickwirkung der Ringspannungen zu _ beriicksichtigen. Diese _ 
Riickwirkung ist vernachlassigbar klein, wenn! der Rohrenfeder- 
kennwert 

i a 
A= dR, Sel 
ist, wobei a die mittlere groBe Halbachse des Ovalquerschnittes, — 
d die Wandstarke und R, den Kriimmungshalbmesser bedeutet. 

Diese Bedingung ist bei Hochdruckréhrenfedern stets erfiillt. 
Man hat also die Vereinfachung, daB zuniachst die Formande- 
rungen des Querschnitts genau so berechnet werden kénnen, wie — 
wenn die Rohrachse gerade ware. Erst hinterher werden die — 
Ringspannungen und die Kriimmungsanderung der Feder er- | 
mittelt. 


Abb. 1. Schematische Darstellung einer : a é 7 
= Figs ee Ts Den Spannungszustand in einem geraden Rohr kénnen wir 


aber auf einen ebenen Formanderungszustand zuriickfihren. 
Denken wir uns eine ,,Scheibe* herausgeschnitten, so kénnen sich die parallelen Ebenen wegen 
ihres Zusammenhanges mit dem tibrigen Korper nicht ,,ausbeulen“*. Wenn wir die x-Achse mit 
der Rohrachse zusammenfallen lassen und eine Torsionsbeanspruchung des Rohres ausschlieSen, — 
lauten die Gleichgewichtsbedingungen? 


O0y 
oy 


Cie 
iy ah ee 2 es ee 
1 Oz ap Me (1) 


Diese Bedingungen lassen sich durch Einfiihrung der Airyschen Spannungsfunktion F befriedigen: | 


ee Cr er 


Cina s Da Bye? ee eRe ce (2) Z 


Die Spannungsfunktion mu dabei der folgenden biharmonischen Differentialgleichung geniigen: 


eer is ia Sas a4 Coe 
A AR a aad 5a) F= 0. (3) 


Fur die axialen Spannungen erhalt man 


Ox =" (Gy +0;) tay +Bzt+y, (4) 


wobei a, 8 und y Konstante sind und 4 ~0,3 die Querkontraktionszahl bezeichnet. Die drei letzten 
Glieder entsprechen einem Spannungszustand, der aus Zug und Biegung zusammengesetzt ist. Wir 
setzen zunachst ~—/ =y =0 und iiberlagern erst am SchluB® diese Zusatzspannungen. 


* Eine Begriindung dieser Bedingung wird an anderer Stelle gegeben werden; vgl. auch W. Wuest, Ing.- 
Arch. 20 (1952), S. 116. 


* Vel. C. B. Biezeno u. R. Grammel, Technische Dynamik, 2. Aufl., Bd. 1, S. 104. Berlin, Gottingen, 
Heidelberg 1953. 
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Hinzu treten noch folgende Bestimmungsgleichungen fiir die Verschiebungen: 


ov 
Oy = Er, b= Oy — - (o, Ie Oz). (5a) 
é 
E = = &,E =o,—p(ox +0,) (9b) 
ov Ow 
Cz 4 = Ty: (5c) 


Dabei ist der Schubmodul G = E/2 (1 + p). 


2. Der Spannungszustand im geraden Rohr von flachbogenférmigem Querschnitt. Bei der 
Herstellung von Réhrenfedern geht man von einem dickwandigen Rohr von kreisférmigem Quer- 
schnitt aus, dem man durch Flachschlagen oder Flachpressen in heiBem Zustand die gewiinschte 
ovale Querschnittsform gibt. Wir bezeichnen diese 
Form als Flachbogen-Querschnitt und ideali- 
sieren sie entsprechend Abb. 2. Wir zerlegen den 


ad 
0-74 $4 


Abb. 2. Zerlegung des Flachbogenquerschnittes in die Abb. 3. Krafte und Momente 
rechteckigen Mittelteile und die Randbégen. am Rechteckteil. 


Querschnitt, wie in Abb. 2 angedeutet, in zwei halbe Ringflachen und zwei Rechtecke. An den 
,,Nahtstellen“ miissen die Spannungen und Deformationen stetig bleiben. Durch Ansatz von 
Reihénentwicklungen kann man diese Forderung sicher erfiillen. Wir wollen aber einen anderen 
Weg gehen und suchen eine Naherungslésung dadurch zu gewinnen, daB wir fiir den Ringbereich 
und den Rechteckbereich einfache Ansatze fiir die Spannungsfunktion wahlen, welche die Rand- 
bedingungen am Auffen- und Innenrand erfiillen, wahrend an der Nahtstelle nur die Integral- 
werte der Spannungen, Momente und Deformationen (iiber die Wanddicke integriert) stetig sind. 
Wir machen uns in diesem Fall das Prinzip von De Saint Venant zunutze, nach dem Stérungen und 
Spannungsverteilung in einer gewissen Entfernung von der Lastangriffsstelle abklingen. Wenn 
man das Ergebnis anderer bekannter Falle tibertragt, diirften die Stérungen etwa in einer Ent- 
fernung von der GréBenordnung der Wanddicke abklingen. Man wird diese Naherung also unter 
Umstanden fiir sehr dickwandige Federn nicht mehr anwenden diirfen. (Vgl. hierzu Abschnitt 2e.) 

a) Rechteckteil. Im Rechteckteil (Abb. 3) wahlen wir fiir die Spannungsfunktion das 
Polynom 


F = a,2* + a, + ayz* + a52° + -y” (by + 6,2 + baz” + 6, 2°). (6) 
Durch Einsetzen dieses Ansatzes in (3) folgt dann 
| b= 34,). by = —54:, (6a) 
und man erhalt fiir die Spannungskomponenten 
Gy = 2a, + 6a,z + 124, (22 — 0,5 y7) + 204, (2 — 1,9 y*2) , (7a) 
6, =2b, + 2b,2— 6a,2? — 10a,2', (7b) 
Tyz =—2b,y + 12a, yz + 30a; y2?. (7c) 


Die Randbedingungen lauten 


ar a 


(8) 
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i 
Die vier Bedingungen ergeben i 
{ 
= 2S ; 
Gea Or by = 
as 


P 
4 
Le io ; 
h 


Nun ist aber die gesamte Langskraft offenbar gleich der von den Halbkreisenden ausgeiibten Kraft, 
also : 


oP 
~ 403” @ 


h | 
f{oydz =4ha,=pr. (10) | 


—h 


Es bleibt somit nur die Konstante a, statisch unbestimmt. Wir kénnen sie noch durch das Mo-— 
ment M, ausdriicken, das an der Nahtstelle ibertragen wird: 


h < . , i 
Mee tet. 3 M, Ein mbeke ee: 
[_ fovea =— My = 40,8 + (5 —z]p. 64 =a (it ao we ayy 


—h 
Das Moment M, ist dabei statisch unbestimmt. 
Insgesamt lauten jetzt die drei Spannungskomponenten 


ite) If Ne a EH Visoe eee es ensa, l | 
ey =|retalg) +2 oh OR |S eee: (i2— 


[te 1h fae 3 P—y?/3 zy Be 2 3 M 
poet DN eairrdica he ee 2m | 


ee) ie = ia (a) [PX +A 


Die bei der Integration auftretende willkiliche Funktion F(z) mu aber verschwinden, da aus | 
und (12) fiir die Verschiebungskomponente w durch Integration 
1 ab |e I /f eaXe 
Ew=(—)|—3 Bk v(t) | ie 
r L/c2Ne 3 2—y?* z a i 
rg recrgea ears ateerseriya pote ee tar): 
Die zunachst willkiirliche Funktion F,(y) bestimmt man aus (5c) zu: 


S458 pei 3 Mer | 
F,(y) = — |0,0925 gh Saye) oe) YBa lie) a Va en (14a) 


wobei cy eine spater zu bestimmende Integrationskonstante ist. 


_ b) Gekrimmter Abschnitt. Es empfiehlt sich hier, Polarkoordinaten zu wihlen (Abb.4). | 
Die Differentialgleichung (3) der Spannungsfunktion schreibt sich dann 


picay Sete ie bee Ay ¥ 
AA F=(q r Ot | Te ae) ai (159 
und die Spannungskomponenten lauten 
= a Or ior ~*~ CF Gaelok 
"72 gy? r or?) 99 = Be? cris (16) | 
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Fir die Dehnungen gelten die Bezichungen 


Ee, =6,—p (0, +0.) = ES, | 
r 
Ri atcwee ees eapnpAar 1 0Aq | 
¢ = % —[ (ox +4,) B( Pree, | (17) 
= 1 a4 aA 
CYro — Tr = G(; i | z zr) . 


Wir haben dabei die Verschiebungen eines Punktes (r,p) mit Ar, Ap bezeichnet. Es ist zweckmabig, 
von der allgemeineren Spannungsfunktion zunachst den achsensymmetrischen Belastungsfall abzu- 


trennen. Die Spannungsfunktion F, fiir diesen Sonderfall ist dabei als Lésung der Differential- 
sleichung 


f+ ent Hee AC 
44 =(at; iz) Fi =0 (15a) 
gegeben durch ; 
k= €& + ¢,Inr + ¢,r? + cr? Inr. (18) 
Die Spannungskomponenten werden dann 
1 dF. 
Or1 ee te + ¢,(1 42 Im r) +t : | 
dF. 
Col =F = —F+2q+e,(8+2Inz), (19) 
Tr gl — 0, 
and die Konstanten c,, ¢,, c, bestimmen sich daraus, daB 
Ted ot ea Cope ere 
) fUY, f= Ty 5 Or1=—Pp, 
and. 5 


2 


; 
[ee rdr M,+pP 2 (7; + 72) Bild 4. Krafte und Momente am 


zy Randbogen. 
" 


‘ein mul}. Hierbei ist M, das an den beiden Enden iibertragene Moment (Abb.4). Die drei Be- 
ingungen liefern 


1 60 : 5.9/9 P 
cy = [2rdln 2 (rep —2 Mo) — BA pl | 
Ly a Ss a 2 2 2 ‘ | 
| C= Ny 2 r— 1, +2 (r3lnr, — rj ln r) (pry %2—2M,) 
| : (20) 
+ (1 +2Inr,) r? rz In a ae ri) P , 
1 
lees ‘ ire | 
cy =| CE FD (Prete —2 My) — 2 rim 
N = (r3 — r?)? — 4 rir} In? 4 : (20a) 
fiir die Dehnung bzw. die Verschiebungen erhalt man durch Integration der Gleichung (17) 
| BAr, = (1 +1) 4 + (1 2p) 2eyr +e,(2 1—2p)rInr—r)], | (21) 
— EAg, = 4 (1—p*) ¢; (p — 9%) - | 


Wegen Ag, =0 fiir py =2/2 muB oy = m/2 sein. 
r2 

Nun ist zwar f % 9; dr =pr,, und die an der Nahtstelle (y =0) iibertragene Langskraft wird vom 
achsensymmetrischen Ansatz erfaBt. Wegen des durchgangigen Verschwindens der Schubspannung 
‘ehlt aber in dem Ansatz der EinfluB der Querkraft V (Abb.4). Wir miissen also noch einen An- 
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teil hinzufiigen, der die iibertragene Querkraft erfa®t. Ein solcher Ansatz ist gegeben durch? 


F, =fylr) sing. (22) 
wobei man durch Einsetzen in (15) fiir f, die gewohnliche Differentialgleichung 
Gos Al iG! 1\2 
(za ied =) fe =9 (15b) 
erhalt mit der Losung 
fa =qr +eo;rlnr+or+—. (23) 
Die Spannungskomponenten werden dann 
a 2 Ghee 
Or2 =(4) sin @ =(2ar+3— 3 sing, 
Woe 6 Cs 74 Cy . 2 i 
Gan. fs al CG peat stage sing, ( ss) 
d c 2¢ | 
yeep a7 cos g = — (2ear+ 3 — a cosy. | 


Fir g=0 und x verschwinden o,2 und o,2, so daB wir diese Stellen als freie Enden betrachten — 
konnen, Gleichzeitig haben die Schubspannungen ihren GréBtwert. Dies entspricht also einem _ 
Belastungsfall, bei dem am freien Ende eine Querkraft iibertragen wird. Da die Radialspannunga,. | 
sowohl fiir r =r, als auch r =r, verschwinden soll, kénnen zwei Konstante eliminiert werden: 


r 


2 2 22 
ry + ro Tier aNe « 
Org ne Cy (* = 73 sing, 
04 2 ree 
n+re LET OLN 7G 
Oy2 a0, (3,458 — Ail sing, 


itn an 
Trp2 =—2a(r <= y + 1s} e08 


Die Konstante c, bestimmt man daraus, da an der Stelle ~ = 0 


[tn dr = V = pl 
sein muf} (vgl. Abb. 4), woraus folgt 
oe pl ae 25) 
“4 ~ 21@t Fre) arin —(— A)" Ce 
Fiir die Verschiebung erhalt man aus den (17) 


2 
r 


Late . db 
Ban =| 0p} + yan @ to 3 
(26) | 


E Ap, =— {7 (1) +21) 45 cog — P+ RO. | 


Die bei der Integration auftretende Funktion R(r) muB aber verschwinden wegen Ag, =0 fir 
y= 2/2, und die Funktion ® bestimmt man mit Hilfe der dritten Gleichung (17) zu 

@ = U, cosy + Uy sing — 2c, (1 — 2) sing. 
Da Ag, = 0 fir p =x/2 ist, muB U, = c, (1 —y?) sein. 

Das Glied mit U, entspricht lediglich einer Parallelverschiebung in der Richtung der langen 
Achse, die nicht interessiert, so daB U,=0 gesetzt werden kann. 

c) Ermittlung des statisch unbestimmten Momentes. Das an der Stelle py =0 iiber- 
tragene Moment M) ist bisher unbestimmt geblieben. Wir ermitteln es aus der Forderung, daf der 
Mittelwert der Verdrehung zu beiden Seiten der Nahtstelle gleich ist: 

h 


aff Lesele Bele: 
milf) all or “) (29) . 
r, g=0 | 


—h 


‘ A. und L. Féppl, Drang und Zwang 1, S. 307. Miinchen 1941. 
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ist. Aus (13) findet man aber mit u = 0,3 


h 
: a Lee l “ M, ‘ 
xi (| se ae] =+ (0.455 je se 0.013) p= 1,365 ae (28a) 
—h yy=l 
Andererseits liefern (21) und (26) 
jie (ak eA) i . 
eer er Tey [(r Ap), — (dp), ] 
e p=0 (28h) 
5 = — 1,82 4 C3 Saeed [6.24 (r5 a r;) -+- 1,04 ee In al A | 
ty 
2| be 
i a 74 /4=05, 
. 40 1a 
Bey 
-4 ee] 
40. 
4 
Zor 
-3 


Abb. 5. Das an der ,,Nahtstelle** zwischen Rechteckteil und Randbégen iibertragene Biegemoment in Abhangigkeit 
von den Abmessungsverhiltnissen. 


.us der Gleichsetzung von (28a) und (28b) und durch Einsetzen der Konstanten c, und c, aus (20) 
nd (25) erhalt man fiir das Moment M, bei einer Querkontraktionszahl , =0,3 
3 
3,64 (—) he 2 eG 
= i =p (29) 
10,92 7) +H 


M, 
d2 


@ == hea ty Ty 5 


6,24 (rg — rf) + 1,04 (r} + 73) In (“2 
, 1 


F = 0,026 + 


v) 


2|(% + 14) In 2 — ( — ri) 
dig 


1,82 a 2.2 l9 
Ge lc — ri) r4yrg—2 re In ("| 4 
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wobei N durch (20a) gegeben ist. Dieses Ergebnis ist in Abb.5 dargestellt. 


d) Gré8te Beanspruchung. Nach Anschauungen, die insbesondere von A. Nadai* ver- 
treten werden, ist die Beanspruchung bei dickwandigen Rohren unter Innendruck nicht durch 
die gréBte Hauptspannung oder Hauptdehnung gegeben, sondern durch den Gréfitwert der spezi- 
fischen Formanderungsarbeit. Man erhalt danach fiir die FlieBgrenze 


Ge /+ [(oy a4)? + (62—o3)? + (¢3—4;)"] - (30) 


Im vorliegenden Fall ist die Beanspruchung an der Innenseite der Randbégen, also fiir r =r, und | 
p =n/2 am gréBten, weil dort Gréfitwerte von Biege- und Zugspannungen zusammentreffen. Es | 
ist dort 


Oy = 0, =— Pp, (31a) | 
am ——fbta ta +2inn)—4ieg—z) | 
Taupe (31b) 
= —— p =| 2(e 3) ae (rz — ri), | | 
ce Oy = Ox = fb (Or +49) +> (31e) | 
wopel | 
FB x 
y=(1—2 yu) FP, (32) | 
1 | 


Fy, =r?a +A4lr, die lichte Innenflache, 
F, = (2? — ri)a +41 (r, —1,) die Querschnittsflache 
ist. Formel (30) kann noch mit Hilfe von (3lc) umgeformt werden zu 
6, =V(l—p + pw) (0, — 0)? +(1— 4 + 4p?) 0,0, +7? —y (L—2 4p) (0; +0). (33) 


Zahlenergebnisse fiir das Verhaltnis der gréBten Beanspruchung 6, ma, zum Innendruck p sind 
in Tabelle 1 in Abhangigkeit von den Parametern r,/d und I/d zusammengestellt. 


Tabelle 1. Gréfte Beanspruchung Coma 

e Le l 
S 0,25 0,5 1,0 1,5 2,0 2,5 

Id a. | 
0 1,7918 1,9486 2,3102 2.7005 3,1579 3,9350 
0,1 2,3382 2,4172 2,8165 3,2912 3,7878 4,2944 
0,25 SRB) | SAD 3,6801 4,2637 4,8889 5,5298 
0,50 5,2106 4,9058 Dp | 6,1419 6,9725 71,8347 
0,75 7,5681 6.9396 | 71,4172 8,3277 9.3559 10,4352 
1,00 10,3526 | 9.3303 9,7615 | 10,8126 12,0303 13,3236 
leo 132502500 12,0742 12,4142 13,5909 14,9924 | 16,4944 
1,50 17,1968 15,1669 15,3710 16,6580 18,2353 19,9430 
2,00 25,7360 | 22,3914 22,1846 23,6453 BSH! | 27,6607 
3,00 47,8906 | 40,9651 39,3613 40,9736 42,4531 | 46,3117 
4,00 76,8014 65,1729 | 61,2163 62,7051 65,6285 69,1663 


In Abb.6 sind fiir eine angenommene Streckgrenze von 100 kg/mm? die maximal erreichbaren 
Driicke aufgetragen, wenn an keiner Stelle die Streckgrenze tiberschritten werden soll. Fiir héhere 
Streckgrenzen sind die Druckwerte mit einem entsprechenden Faktor zu multiplizieren. Vielfach 
werden jedoch noch Beanspruchungen zugelassen, bei denen 6rtlich die theoretische FlieBgrenze 
iiberschritten wird. Beim erstmaligen Unterdrucksetzen der Feder tritt dann in einem gewissen — 


1 4, Nadai, Der bildsame Zustand der Werkstoffe, Berlin 1927. 
2 Vel. zB. W. Wuest, Z. angew. Math. Mech. 32 (1952), S. 90. 
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Bereich (Abb.7) plastisches FlieBen ein, so daB bleibende Verformungen auftreten. Beim Ent- 
lasten der Feder bleiben innere Spannungen zuriick, die den Betriebsspannungen entgegenwirken 
(,,Autofrettage der Feder). Weiteres FlieBen tritt jetzt nur dann ein, wenn der beim ersten Mal 
erreichte Druck iiberschritten wird. Bei der Betrachtung von Abb.6 fallt auf, daB die Kurve 
r,/d=0,5 fiir groBe Werte von I/d die anderen Kurven schneidet. Dies ist offenbar auf eine Span- 
nungserhéhung durch ,,Kerbwir- 
kung“ zuriickzufihren. 


j theor Beanspr 


HA AN ®E 


Bild 7. Autofrettage der Réhrenfeder. Durch 

Ufa Einleiten eines hohen Druckes wird in der am 
G7 015 «G2 03 G4 G8 O6 G708097 45 iy iB X héchsten beanspruchten Zone die FlieBgrenze 
uberschritten. Die Spannungen wachsen dann 

Abb. 6. Verhaltnis von Héchstdruck zu Héchstbeanspruchung in Abhangigkeit von nicht mehr weiter und der Werkstoff wird pla- 


70 


den Abmessungsverhaltnissen. Fiir eine héchstzulassige Beanspruchung von 100 kg/mm? stisch verformt. Beim Entlasten ergeben sich 
gibt die Ordinatenskala unmittelbar den Druck in kg/em? an. Die Kurven gelten fiir bleibende Formanderungen und innere 
rein elastische Verhaltnisse, also ohne Beriicksichtigung einer Autofrettage. Spannungen. 


Ofp 
-7 — b. =I 
a ” 7 2 3 rp ae a Ba “7 0 1 2 
Bild 8a. Tangentialspannung zu beiden Seiten der Nahtstelle nach Bild 8b. Tangentialspannung zu beiden Seiten der Nahtstelle 
den Naherungsansitzen (r,/d = 1; l/d = 1). nach den Naherungsansatzen (r,/d = 0,5; I/d = 1). 


e) Spannungsverteilung an der Nahtstelle. Wegen unserer Naherungsannahmen hat 
die Spannungs-. und Deformationsverteilung an der ,,Nahtstelle“ zwischen dem geraden und dem 
gebogenen Abschnitt eine Unstetigkeit. Um uns ein Bild ttber die GréBe dieser Unstetigkeit zu 
machen, ist fiir r,/d = 1, l/d =1 sowie r,/d =0,5 und I/d = 1 die Verteilung der Tangentialspannung o, 
bzw. o, zu beiden Seiten der Nahtstelle berechnet worden. Wie Abb.8a und 8b zeigen, ist der 
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Spannungssprung der Naherungslésung nur gering, vor allem wenn man ihn mit der maximalen 
Tangentialspannung an der Stelle my =z/2 vergleicht, die in den beiden Fallen 10,12 p bzw, 9,61 p 
betragt. Man kann also vermuten, da die Naherungslésung auch noch in extremen Fallen ein 
einigermaBen zutreffendes Bild der Spannungsverteilung liefern wird. 


3. Kriimmungsiinderung der Feder. Wegen der Kriimmung des Federrohres rufen die Quer- | 


schnittsdeformationen, wie in Abb.9 schematisch angedeutet, ,,Ringspannungen“ hervor. Wenn 
die Querschnittsabmessungen klein gegen den Kriimmungshalb- 
messer sind, kann man naherungsweise mit einem mittleren 
Kriimmungshalbmesser Ry, rechnen und die Ringspannungen 
sind dann gegeben durch o; = E w/Ry, wenn w die Deformation 
in der radialen z-Richtung ist. Die Ringspannungen rufen 


Kriimmungshalbmessers nachgibt. Bezeichnen wir die Kriim- 
mungsdnderung mit ~ =A R/ Rp, so ist das resultierende Moment 


M=||Exza—|[zo% ar=o 
le Ry i Tae 3 


J 


wenn die Feder nicht durch adufere Krafte belastet ist. Die 
Integration ist dabei itber die gesamte Querschnittsflache F, 
zu erstrecken, und z bezeichnet den Abstand von der Sym- 
metrielinie, also z = z + r, + d/2 im geraden Abschnitt und 
= —recos g im Randbogen. Daraus folgt aber fiir die Kriimmungsanderung 
w2dF 
0 =e 


Bild 9, Ringspannungen o,, infolge der Kriim- 
mung der Réhrenfeder. 


(34) 


In Abb. 9 ist angedeutet, da die Ringspannungen durch radial gerichtete Spannungen O; 
im Gleichgewicht gehalten werden miissen. Damit wirken also die Deformationen vermittels der 
Ringspannungen auf den urspringlichen Spannungszustand zuriick. Bei dickwandigen Federn ist 
der Einflu8 dieser Riickwirkung nur gering und kann im allgemeinen vernachlassigt werden. Bei 
diinnwandigen Federn muf dagegen die Riickwirkung beriicksichtigt werden, wie eingangs erwahnt 
wurde. 

a) Berechnung der Kriimmungsanderung. Bei der Berechnung des Ausdrucks (34) zer- 
legen wir den Integrationsbereich in den geraden und den gekriimmten Abschnitt. Im geraden Teil 
ist w durch (14) und (14a) gegeben. Die Integrationskonstante c, ermitteln wir aus der Bedingung 
daf an der Nahtstelle y = 1 bzw. y = 0 

h To 
ftoids =aaf Ardr 


—h Ti 
sein muB, wobei Ar = Ar, + Ar, mit Ar, aus (21a) und Ar, aus (26). Man findet so 


4 
yee 2.275 (7) + 1,384 Be 0,05043 d 


1\2 M, . | 
— |5:46(7) — 05] a5 0 | ae SE 04 ree. (35) 


L 


rlnr,—r? Inn, 


ai: os 
Vs a ry 
Hierbei ist wieder 4 = 0,3 gesetzt. 
Im gekriimmten Abschnitt ist 5 
w = Ar cosg—r Ag sing 
mit 
Ary aw nraey lie lp = Ag, + Ag, . 


Setzt man noch fiir das Tragheitsmoment der Querschnittsflache um die groBe Achse 


Niaes 4 ; g 
[pH gE = Delany ee Gao Oi) (36) 


1 


ein Moment hervor, dem die Feder durch VergréBerung des | 
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so erhalt man schlieBlich fiir den Ausdruck (34) 


E 1/406 dS d\s 
Oo = sla Bs 0,819 — 0.5266 (7) + 0,0561 (7) — 0.0325 () a 


; d\'d M 
0.0379 (7) | +4 Mons 


; fe 
ene 0.5 (7) | + 4e bdl 


IRENA Dy fe eee 
es E 2D 6, + 0,200 (r$—r?) cg ey) 


+ (0,200 (73 In rg — rj In7,) + 0,0500 (r4 — ri)) c 


‘3 


— (0,59982 (r3 — r3) + 1,7542 r? 73 (r2— 1,)) “ | ay 

j Pp 
Dabei ist b =r, + d/2 die mittlere kleine Halbachse des Querschnittes. Der Summand jy/p am 
Schluss von Formel (37) ist dabei in folgender Weise entstanden: In den Formeln fiir die Defor- 
-mationen sind bisher die linearen Glieder in der axialen Spannung Ao, = f z + y vernachlassigt 
worden, wobei j = (1 — 2) (Fy/F,) p und 8 = «w E/R, ist. Diese Anteile liefern fiir die Deforma- 


tionen v* und w* 


av* (Quay Ow* 2) 
— - | ara es Zz ua 
oy uaF ar az lg? + &) 
also 
| <a Hata aca Ty ieee OS 
a u (F zy z=) t+ V(z), ey (i ze i) + W(y). 
Dabei ist aber V(z) = 0, weil v* = 0 fiir y = 0 wird. Ferner findet man aus 
dv* Ow* 
a renaye Te 


die Beziehung 
oD) 


Wy) = = R y? + const , 
so daB schhieBlich 


F, Fy 
wird. Durch Auswerten der Formel (37) ist die Tabelle 2 entstanden. 


Tabelle 2. Kriimmungsdénderung  E/p. 


i 7 7 
r,/d | | | 
0.25 0,50 | 1,00 | 1,50 2,00 2,50 
| | 

0 0,1210 0,3425 0,8267 1,3005 1,7610 2,2087 
0,1 0,2329 0,5312 1,2115 1,9014 | 2,63607 | 3,4183 
0,25 0,5551 0,9943 1,9554 3,0309 | A,2817 | 5,9675 
ya55 Ibias, | 2,1305 Sloan =| 5,0394 | 71,8004 | 10,136 
0,75 2,7650 4,1099 6,4126 9,2333 12,450 16,093 
1,0 4.9370 6,6148 LOSERS Hae 14,075 18,561 23,622 
11715, 8,1309 10,389 15,084 | 20,339 26,903 32,911 
ies 12,613 15,522 7 2933 28,238 agi 44,158 
2 26,642 30,915 39,849 | 49,801 60,974 | 73,347 
3 853431 91,074 105572 One eS LLO 425.94: | 176,398 
4 205,878 210,281 227,162 Psy WAH 5) 280,040 353,199 


| der Querschnitt symmetrisch in bezug auf die Liniez = Oist. Dies trifft aber nur in erster Naherung 
gu. Bei der Herstellung von Hochdruckfedern wird namlich das urspriinglich gerade Rohr in 


H 
- rotgliihendem Zustand (also bei niedriger FlieBgrenze) gebogen. Dabei werden die dem Kriim- 


b) Exzentrizitatskorrektur. Bei der bisherigen Berechnung wurde angenommen, daf 


i 


mungsmittelpunkt zugewandten Fasern gestaucht, wahrend im auferen Bereich die Wandstarke 
8 
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verringert wird. Der Querschnitt wird also unsymmetrisch. Wie in einer friiheren Veréffent- 
lichung! nachgewiesen wurde, tritt aber bei Rohren mit unsymmetrischer Querschnittsform unter 
der Wirkung eines Innendrucks ein zusatzlicher Kriimmungseffekt auf, der im vorliegenden Fall 
dem urspriinglichen entgegenwirkt. Zur Berechnung dieses Effektes wollen wir die vereinfachende | 
Annahme machen, da® die Verformung der Feder ideal-plastisch vor sich gehen soll. Dann gelten | 
die Beziehungen 
: | 
Tea BAR RF ae iB SA 
Nun ist aber 
8 ee BE 
"oT Rgt are ay oe ay ane 
also 
1 yz = ih Ze 
=— > =———4+W 
v FRE e w 4 ue (v) 
Wegen v = 0 fir y = 0 muB V(z) = 0 sein. Ferner gilt 
ov Ow 
te By Gie 
so dah 
ly De es 
SScen: + Wy) = 9 
und demnach : 
NB Sipe ieee 
peta nS Be He Xi R, (38) 


wird. Diese Formeln gelten natiirlich nur, solange die Verschiebungen klein gegeniiber den ur- 
spriinglichen Abmessungen sind. 


Hs Anwendung auf den Kreiszylinder. Hier ist 
Ar=veosy+wsing, y=rcosgp, 2z=rsing. 


Damit ergibt sich 
2 


r . 
Near ee (39) 
Die zusatzliche Kriimmungsanderung ist aber gegeben durch 
FAs p 
wo, = (1— 2p) tee (40) 


wobei Fy der lichte Innenquerschnitt, J das Tragheitsmoment des Wandquerschnittes in bezug 
auf den Schwerpunkt, As der Abstand zwischen dem Schwerpunkt des Innenquerschnittes und des 
Wandquerschnittes ist. Die Verschiebung des Innenquerschnittes ist aber nach (38) 

a as 


ee ee 
1 aR, 
und die des Wandquerschnittes 
sees r3—?r} 5 re + 
A Ry (r3 — r?) ATR 
so daB man 
2 
As =s,—s,= ik 


hat. Damit wird die zusatzliche Kriimmungsinderung (bei Beschrankung auf Glieder erster Ord- 
nung) 


el) tae (41) 


ri—ri E 


2. Anwendung auf den Flachbogenquerschnitt. Wir teilen hier den Querschnitt wieder 


in den geraden Innenteil und die Randbégen ein. Die Verschiebu ngen im Innenteil sind durch die 
Beziehungen (38) gegeben. 


' W. Wuest, Z. angew. Math. Mech. 32 (1952), S. 90. 


7 
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In den Randbégen setzen wir fiir die innere Begrenzung (Abb. 10) 
y=l+r, sng, Z= 1, cosg, 
so dafi die Beziehungen (38) dort die Form 


bie 1 r, cos p(l + 7, sin p) 


2 Ry ‘ 

1 (J +r, sin y)? — rj cos p 

(i) = ~ = 
4 Ry 


annehmen. Die radiale Verschiebung in den Randbégen ist dann 


Ar =vsing + weosq, 


so dal} 
l : 
Ar = aR, (I? a ri) cos YP 


wird. Das Flachenmoment der deformierten flachbogenférmigen Querschnittsflache in bezug auf 
die alte Symmetrieachse ist dann 


M(r, 1) =+ Eras ! 


2 (j2__ 72 
R, | 3 r (1 NG 


7 
A 
wobei r = r, fiir die innere Begrenzung, r = r, fiir die 


auBere Begrenzung gilt. Fiibrt man noch die innere 
Querschnittsflache 


Po=4lr+ra 


und die Wandquerschnittsflache 


=S) 2 2 
F, = 41(r, —1) + (2—14) 2 “ 
Bild 10. Plastische Formanderung (iibertrieben) beim 
: 2 : ‘ : Herstell i der Ré : 
ein, so betrigt die durch die Deformation des Quer- ne apc give Ota Naa aN ions Bk 
schnittes beim Herstellungvorgang verursachte zusatzliche Kriimmungsadnderung beim Unter- 


drucksetzen der Feder 


F 
[M(ra, 2) — M(ry D1 Ms) 
eo = — (l— 2p) a ee: R, =. (42) 
In Tabelle 3 sind Zahlenwerte fiir die GréBe — w, E/p zusammengestellt, die durch Auswerten 
von Formel (42) gewonnen worden sind. 


Tabelle 3. Exzentrizitdtskorrektur w, E/p. 
a es = See ee ee 


r,/d | | 
0 0,0160 0,0450 0,1067 | 0,1654 0,2215 0,2759 
0,1 0,0216 0,0518 0,1098 0,1722 0,2282 0,2925 
0,25 0,0279 0,0601 0,1228 0,1814 | 0,2374 0,2918 
0,50 0,0344 0,0697 0,1349 0,1945 | 0,2511 0,3056 
0,75 0,0377 0,0754. 0,1435 0,2050 0.2624 | 0,3178 
1,00 0,0387 0.0771 0,1495 | 0,2144 0,2718 0,328] 
1,25 0,0382 0,0792 ~. 0,1530 | 0.2188 0,2793 0,3368 
1,50 0,0367 0,0784 Osaal 0.2230 =| 0,2852 0,3441 
2,00 90,0313 0,0735 . 0,1547 | 0,2270 0,2930 0,3547 
3,00 ~ 0,0153 0,0545 0.1418 | 0,2222 0,2952 0,3634 
4,00 —0,0043 0,0304 OUL90 Ny 0,2094 0.2852 0.3593 


4, Winkeldrehung und Federweg. Technisch kann man sowohl die Winkeldrehung als auch die 
Verschiebung des Federendes ausnutzen. Um beides zu berechnen nehmen wir an, daB die Feder, 
die im Ausgangszustand einen Kreis mit dem Halbmesser R, und dem Windungswinkel y, bildet, 
unter der Wirkung des Innendruckes in einen gréferen Kreis mit dem Halbmesser R und dem 


8* 
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Windungswinkel y iibergeht. Wenn die Dehnung in der Achsenrichtung des Rohres mit ¢ bezeichnet 
wird, mu offenbar fiir die Federlange gelten 
Ry — Ry po = € Ro Y- 
Ferner gilt 
R = ee ql ae Wk) ° 
wobei w, = w + w, die wegen des Exzentrizitatseinflusses korrigierte Krimmungsanderung ist. 
Demnach gilt 


1 ; 
W = Wo i ~ Wo (1 + ¢— oz), 
und die Winkeldrehun g ist : 
Ay = Wo — YP = Yo (we — €) - (43) 
Dabei ist die Langungskorrektur [vgl. (31 ¢)] 
F 
UE 6 Oe wee Sy ie) eae) ee) ae (44) 


Durch Auswerten dieser Formel ist Tabelle 4 entstanden, in der Zahlenergebnisse fiir ¢ E/p in 
Abhangigkeit von r,/d und I/d zusammengestellt sind. 


Tabelle 4. Ldngungskorrektur ¢ E/p = y/p. 


r,/d 
0,25 | 0,5 1,0 the} | 2,0 2,5 al 

I/d Sh | 

— 
0 0,0167 0,0500 0,1333 0,2250 0,3200 0,4167 
0,1 0,0232 0,0590 0,1442 0,2366 -0,3319 0,4288 
0,25 0,0313 0,0706 0,1589 0,2526 0,3487 0,4461 
0,50 0,0415 0,0862 0,1800 0,2765 0,3742 0,4726 
0,75 0,0491 0,0985 0,1977 0,2973 0,3970 0,4968 
1,00 0,0549 0,1084 0,2128 0,3155 | 0,4174 0,5188 
129 0,0596 0,1165 0,2258 053.3 7aneul 0,4359 0,5390 
1,50 0,0633 0,1233 0,2371 0,3462 0A527 0,5575 
2,00 0,0691 0,1340 0,2558 0,3709 0,4820 0,5905 
3,00 0,0766 0,1485 0,2827 0,4082 0,5279 0,6436 
4,00 0,0810 0,1577 0,3012 0,4350 | 0,5622 0,6845 


Die Koordinaten des Federendes sind im Ausgangszustand 
XEy = — Ry siny, 
Y¥Eq = Ry (1 — cos yp) 


und im deformierten Zustand 


Xp = — Ry (1 + x) sin (vo ee) P 


+/+ 
= 


Die Verschiebung des Federendes 


$= V (= — x59)" + (ye — yn)" 
wird daher, wenn man nach den kleinen GréGen ¢ und a, entwickelt, 


s = Roy ax |(2 + yo — 2 Yo sin yy — 2 cos yy) — 2 (¥) — Yo Sin Yo) + (,) 8] (45) 
Die Kriimmungsinderung w E/p ist in Abb. 11 in Abhangigkeit von der GréBe I/d fiir verschiedene 
Werte von r,/d aufgetragen. In das gleiche Kurvenblatt sind strichpunktiert die Kurven Ay El, p 
eingetragen, die sich von den anderen nur durch die Exzentrizitits- und Langungskorrektur unter- 
scheiden, Fir kleine Werte von I/d sind die Unterschiede merklich, dagegen fiir 1/d > 1 nur sehr 
gering. Abb. 12 stellt eine Kombination von Abb. 11 mit Abb. 6 dar, in dem hier Ay E/tpo oy tiber 
104 p/o, aufgetragen ist. Fiir eine Grenzbeanspruchung o, = 100 kg/mm? gibt demnach die 
Abszissenskala unmittelbar den Nenndruck in kg/cm? an. Die berechneten Kurven gestatten nur 
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ol/p 
493 (2V DVR 


— op 
———Ayl/ypp 


10 


Ula 


10" 
0 7 2 3 7 


Abb. 11. Kriimmungsinderung w E/p (ausgezogene Kurven) und Winkeldrehung (Ay/y,) (E/p) (gestrichelte Kurven). 
AYE Yo % 


ee 10714, 
“80 100 2 3 YB 6: 6 71000 E 2 3 4 6000 


und Elastizitatsmodul E und héchst- 


Abb. 12. Winkeldrehung Ay in Abhangigkeit yom Nenndruck bei gegebenem Windungswinkel y, 
unmittelbar den Nenndruck in 


zulassiger Beanspruchung o,,. Die Abszissenskala gibt fiir eine zulassige Beanspruchung yon 100 kg/mm? 
kg/cm? an. Die Einhiillende kennzeichnet die bei rein elastischer Beanspruchung erreichbaren Bestwerte, 
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im rechten unteren Teil eine Einhiillende als Optimalkurve einzuzeichnen. Immerhin zeigt diese 
Darstellung, da® der Federquerschnitt mit wachsendem Nenndruck immer kreisahnlicher werden. 
mu. Die Winkeldrehung nimmt mit wachsendem Nenndruck immer mehr ab. 

Es ist allerdings zu bemerken, daf die Kurven nur fiir rein elastische Spannungszustande gelten. 
Wie schon erwahnt, wendet man bei Hochdruckréhrenfedern stets eine Autofrettage an. Hierfur 
ist aber eine groBere Wandstarke giinstiger, und es kann daher vorteilhafter sein, unterhalb der 
elastischen Optimalkurve zu bleiben. ; 

Aus Abb. 13 kann die GréBe ¢/w, entnommen werden und in Abb. 14 ist die durch (45) be- 
stimmte Funktion s/R, w, aufgetragen. Fiir die praktische Anwendung ist es jedoch zweckmafiger, 
den Federweg nicht aus Formel (45) zu berech- 
nen, weil die Bewegung gewohnlich doch nicht 


ee, 
8 Ro Ox 


2 
0 Yo 
0 WG OY 06 08 40 780° 210° 240° 210° 500° 530° 560" 
Abb. 13. Abhingigkeit der HilfsgréBe ¢/m ,, von r,/d Abb. 14. Abhangigkeit der Bewegung des Federendes s vom 
und I/d. Windungswinkel y,. 


am theoretischen Federende abgenommen wird, auf das sich diese Formel bezieht (vgl. Abb. 18). 
Richtiger ist es, zunachst den ,,Polarpunkt’* der Feder zu bestimmen. Dies ist der Punkt, um 
Pry den sich eine am Federende befestigte 
Ebene bei kleinen Druckanderungen 
dreht. Wenn man den Polarpunkt kennt, 
braucht man lediglich den Polstrahl k 
bzw. k’ (Abb. 15) abzugreifen, d.h. die 
Entfernung vom Polarpunkt zu der Stelle, 
an der die Bewegung abgenommen wird. 
Der Federweg ist dann 
Sk Aa, (46) 
und die Richtung der Verschiebung steht 
senkrecht auf dem Polstrahl. Fiir den 
Fall c/w, = 0, der bei dimnwandigen 
Bourdonfedern stets vorliegt, ist die Po- 
larkurve bereits friiher! berechnet wor- 
den. Im folgenden soll der Polarpunkt 
fir beliebige Werte von ¢/w, berechnet 
werden. Wir gehen dabei von den friiher 
abgeleiteten Gleichungen fiir die Rast- 
-polbahn (Herpolhodie) des Federendes 
aus die auch hier giiltig sind: 


BUT 43 a2 O70 


dyp 

Xp = XE st “dy > 
Abb. 15. Lage des Polarpunktes in Abhingigkeit vom Windungswinkel y, Vp = VE— dig 
und von der HilfsgréBe ¢/«,. Ate s dy 


Als Polarpunkt bezeichnen wir denjenigen Punkt der Rastpolbahn, fiir den ¢>0, w,—>0 dh. 
den Momentanpol fiir unendlich kleine Bewegungen). 


i asgya Wuest, Z. Instrumentenkde. 63 (1943), S, 416. 


\ 
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Wenn wir beriicksichtigen, daB 


OL, 


dy =—w(1 rey da , 


so erhalten wir als Koordinaten des Polarpunktes 


Cosa Oy sin y 1 
=— R : Pe ese envied d Peer dS 

Xp OF g, Ope: Yoo = Ro (1 Prag mers 5 , (45) 
Fur 180° = Yo < 360° und 0 <e/a, <0,5 ist diese Beziehung in Abb. 15 aufgetragen. Im 
gleichen Bild ist fiir y) = 270° und ¢/w, = 0 bzw. 0,5 angedeutet, welche Lage der Polstrahl k zum 
theoretischen Federende und der Polstrahl k’ zur tatsaich ichen Bewegungsabnahme hat. Damit 
sind wir in der Lage, die Bewegung des freien Endes der Réhrenfeder vollstandig zu beschreiben. 


5. Vergleich mit der Mitteldruckfeder. Wenn in den hier abgeleiteten Formeln r,/d sehr grof 
wird bzw. d/r, gegen 0 geht, gelangen wir zum Grenzfall der ,,Mitteldruckfeder. Diese Bezeichnung 
soll andeuten, dafs dieser Grenzfall naherungsweise bei 


Rohrenfedern mittlerer Druckbereiche (Nenndruck etwa Tabelle 5. 


zwischen 50 und 200 kg/cm?) verwirklicht ist. Streng ge- a F, F, 
nommen ist namlich der Grenziibergang d/r, —> 0 in unsern 

Formeln nicht erlaubt, weil man dann in den Bereich der 0 1,7776 02427 
diimnwandigen Niederdruckfedern kommt. Zur Inter- 0,1 1,6605 0,3089 
polation ist es jedoch niitzlich, den Fall r,/d—>oo als ne ey eee 
mathematische Grenzkurve einzutragen, der man sich 0.4 1.1757 0.3897 
allerdings nur soweit nahern darf, als die Bedingung 0.5 0.9875 0,3736 
1 0,6 0,7926 0,3363 
a?/dRy< 1 (wobei a =r, + es d+ l 0,7 0,5941 0.2787 
0,8 0,394.4. 0,2024 
die mittlere groBe Halbachse des Querschnittes bezeichnet) te hee ey 


noch erfillt ist. Wenn wir diese Grenzkurve mit den 

fiir endliches r,/d erhaltenen Werten vergleichen wollen, ist es zweckmaBig, als Parameter das 
mittlere Achsenverhaltnis « = b/a = (r, + d/2)/(r, + d/2 +1) einzufiihren und fiir die groBte Be- 
anspruchung o, und die Winkeldrehung Ay folgende Ansatze zu machen: 


= (Gq) Fier) 


A Sai 
=e as (3) 5 Polo r,/d, ps) ; 


| f; (c,7,/a) 


O05 
Wo 
f (a,7/d) 
4 lea 
7/2 =0,05 


a 0 
— 
0 02 O4 06 o8 10 0 2 4 6 (oe 10 
Abb. 16. Hilfsfunktion zur Berechnung der gréften - Abb. 17. Hilfsfunktion zur Berechnung der Winkeldrehung 
Beanspruchung. des Federendes. 


Die Hilfsfunktionen F, und F, haben im Grenzfall r,/d co und fiir u = 0,3 die Werte! von Tabelle 5. 
Wir erginzen die bisher fiir verschiedene r,/d berechneten Werte dadurch, da® wir F,; und F, 


1 Wine Herleitung der Funktionen F’, und F,, wird an anderer Stelle gegeben werden, da die Mitteldruck- 
feder gleichzeitig auch als Grenzfall der dinnwandigen Niederdruckfedern betrachtet werden kann. 
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auch fiir « — 0 ermitteln. Dabei wird F, aus Formel (33) errechnet, wahrend 


4 (ald)? + (ald) + 1/4 
FE; (0, r,/d, 1) eats (1 pe") (r,/d)? + (4/d) + 1/3 


ist. 

Durch Auftragung dieser und der aus den Tabellen 1 bis 4 umgerechneten Zahlenwerte sind 
die Kurventafeln Abb. 16 und 17 entstanden. Man erkennt, daf fiir die Berechnung der Winkel- 
drehung bzw. des Federweges der Grenzfall der Mitteldruckfeder eine recht brauchbare Naherung 
darstellt, solange das mittlere Achsenverhaltnis « < 0,5 ist. Dagegen liegen die Kurven fiir die 
Berechnung der Beanspruchung erheblich auseinander. 


6. Vergleich mit Messungen. Zum Vergleich mit der vorstehend entwickelten Theorie standen 
zwei MeBreihen an industriell gefertigten Réhrenfedern hoher Druckbereiche zur Verfiigung. Die 
erste MeBreihe bezieht sich auf U-férmig gebogene Stahlrohr- 
federn, bei denen das freie Ende entsprechend Abb. 18 gerade 
verlangert ist. In Tabelle6 sind Zahlenangaben fir diese 
Federn zusammengestellt, und der errechnete ist mit dem ge- 
messenen Federweg verglichen. Mit einer solchen Feder wurde 
auch der ,,Polarpunkt“* P in der Weise gemessen, das am Feder- 
ende B ein feines Drahtnetz befestigt wurde. Bei periodischem 
Unterdrucksetzen der Feder bewegt sich auch das Drahtnetz, 
und mit einem auf einem Kreuztisch befestigten Mikroskop kann 
man diejenige Stelle abtasten, die bei der Bewegung in Ruhe 
bleibt. In Abb. 18 ist der so gemessene Wert neben dem ge- 
rechneten eingetragen. Die Ubereinstimmung ist offenbar be- 
friedigend. Vergleichsweise wurde dieser Versuch in der Weise 
wiederholt, daB das Drahtnetz an der Stelle A befestigt wurde. 
Tatsachlich wurde auch, wie es theoretisch zu erwarten ist, fiir 
P nahezu die gleiche Stelle gefunden (Abweichung etwa 1 mm). 


ie 


Pp : 
‘\gerechnet | 
gemessen | 


88 


Abb. 18. Industriell hergestellte U-férmige 


Stahlrohrfedern. 


Tabelle 6. Mefwerte I. U-formige Stahlrohrfedern. 
R, = 44mm; k’=75mm; p=a; E=2,1-10%kg/cm®. 


| | | ] 
P | d | ca | i | qs" Ay ei Ay - 108 | Sger. Sgem. 
kg/em? | mm | mm | mm | kg/mm? ae ¥ | mm mm 
(Ge Gee tee das me ial: 180 | 686 | 162 | 1,58 
1000 | 2,29 | 185) 2,65 113 11,8 5,62 | 1532 1,65 
1600 |, 2,5" | 2,00 | 2,45 152 8,4 6.40 | 1,48 1,58 
2000 Bast NaS, 10 eo 95. lh aealgs 5,8 5,52 | 1,30 1,44 


1 nach Kurvenbl. Abb. 6 berechnet. 2 nach Kurvenbl, Abb. 11 berechnet. 


Die zweite MeBreihe bezieht sich auf Héchstdruckfedern gréferer Abmessungen und einem wirk- 
samen Windungswinkel y = 250°. Bei diesen Federn liegen die elastischen Beanspruchungswerte 
durchweg tiber der FlieBgrenze des Materials. Es sind also ausgepragte FlieBzonen anzunehmen, die 
eine Autofrettage der Feder bewirken. Derartige Federn miissen stets mit einem Sicherheitsschutz 
(z.B. Schutzrohr oder Schutzgehduse) versehen sein, weil nach einer gewissen Anzahl Lastwechsel 
mit einem Bruch gerechnet werden mu und bei dem hohen Fliissigkeitsdruck die Bruchstiicke mit 
groBer Gewalt weggeschleudert werden. In Tabelle 7 sind Zahlenangaben iiber diese Federn 
zusammengestellt. 
Tabelle 7. Mefwerte II. Héchstdruckréhrenfedern. 

Ry = 92mm; k’ = 120mm; y = 250° = 4,36; E = 2,1- 10° kg/cm?. 


P | d ry l ) 


0 Ay Bt | Ay 103 | Sger. Sgem. 
kg/em?~ | mm mm mom kg/mm? ie ee mm mm 
6000 | 5,2 0,95 2530 300 0,82 2,34 1,22 1,14 
6000 | 7,0 1,65 2,39 278 0,88 2,02 1532 0,98 
6000 | 7,0 1,65 2,39 278 0,88 Zeon 132 1,05 
8000 | 4,85 0,78 Phe 354 0,57 2,17 ies} 0,84 
8000 | one 0,65 2,00 386 0,50 1,90 1,00 1,05 
8000 6,4 0,85 2,49 383 0,54 2,05 1,07 1,00 


» Die Werte fiir dy E/y p muBten extrapoliert werden, da die vorkommenden r,/d noch unter den gerechneten liegen. 
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Bei einer Eintragung dieser Werte in Abb. 12 erkennt man, daf offenbar die Abmessungen keines- 
wegs optimal gewahlt worden sind. Bei einer weiteren Héchstdruckréhrenfeder, die bis 8500 kg/cm? 
belastbar ist, liegen besonders sorgfiltige Messungen vor. Die Abmessungen waren hier wesentlich 
ginstiger gewahlt, so daB sich ein o, = 250 kg/mm? ergab und der gemessene Federweg lag 12°, 
unter dem gerechneten. 

Fiir einen Vergleich zwischen Theorie und Messung ist es iibrigens wesentlich, die Stahlzu- 
sammensetzung bzw. die elastischen Materialeigenschaften zu kennen, weil bei hochlegierten 
‘Stahlen, wie sie im Hochdruckgeratebau verwandt werden miissen, die Legierungsanteile einen 
wesentlichen Einflu8 auf den Elastizitatsmodul und die Querkontraktionszahl haben. 


7. Sonderfall der Kreisrohrfeder. Bei diinnwandigen Réhrenfedern ergibt ein kreisférmiger 
Querschnitt keine Bewegung des Federendes. Bei sehr dickwandigen Federn kann dagegen auch 
ein Kreisrohr eine Bewegung des Federendes liefern, die mit den flachgeschlagenen Federn ver- 
gleichbar ist. Die Theorie der Kreisrohrfeder ist in den vorstehenden Berechnungen mit enthalten. 
Da hier jedoch eine geschlossene Lésung méglich ist, wollen wir diesen Fall noch einmal gesondert 
betrachten. 

Fir die Tangential- und Radialspannungen gelten hier die bekannten Laméschen Formeln 


On — a us =e (1 fa ) P> | 
2 1 
2 P (46) 


Ox = ph (o, se O;) a te 
Der Wert von y wird daraus berechnet, dah 


fa, dia? p 


sein mu, woraus 


Nee rere He egs P 


Fur die Dehnungen und Verschiebungen gilt 
Ar 
Ee, =E— =0,—p (G;, ai Oz) 


Die Radialverschiebung wird damit 


Ar =aty a2 rt tp) a od ) (47) 
Die Krimmungsanderung bestimmt man wie friiher aus 
Jf Arr? sin? pdpdr ff Ar r?dr 
OO ffrsint pdpdr — frédr_ 
mit dem Ergebnis 
wo = ea 3 72 + (1—2y) rl te (48) 


Nach Beriicksichtigung der Exzentrizitatskorrektur (41) erhalt man schlieBlich fir die Kriimmungs- 
anderung 


2 
aT 


WkKorr. = oy ae [(2 a 2) r ie (1 pea /A) ry (48 a) 
2 1 
und fiir die Winkeldrehung 
Aimer. ae Perec tls tog 2 49 
|p poem if ee c AWA AE ( ) 


Die durch (48), (48a) und (49) bestimmten Funktionen sind in Abb. 19 aufgetragen. Im gleichen 
Bild ist die Winkeldrehung fiir gegebene Beanspruchungsgrenzen eingetragen und die gestrichelte 


110 Wuest: Theorie der Hochdruckroéhrenfeder Ingenieur-Archiv | 


Kurve laBt ferner erkennen, daf die korrigierten Werte fiir die Winkeldrehung etwa das 0,7-fache 
der nichtkorrigierten betragen. 

Ich habe 1943 einen Versuch mit einer gewendelten Kapillarfeder von kreisformigem Querschnitt 
bei Driicken bis 2000 kg/cm? durchgefiihrt. Der AuSendurchmesser des Querschnittes betrug dabei 
2,5 mm. Die Winkeldrehung wurde durch einen aufgesetzten Spiegel gemessen. Kin weiterer Ver- 
such wurde im Juni 1945 mit einer fiir 8000 kg/cm? bestimmten dickwandigen Kreisrohrfeder groBer 
Abmessungen unternommen. Die Versuchsunterlagen sind verloren gegangen. Doch betrugen 
nach meiner Erinnerung die gemessenen Winkeldrehungen etwa das 0,55-fache der nichtkorri- 
gierten gerechneten (die Korrekturen waren damals nicht bekannt). Da die hier angegebenen Kor- 
rekturen den Ausgangswert nur auf etwa 
das 0,7-fache erniedrigen, bleibt noch eine 
kleine Differenz zwischen Theorie. und 
Messung. , 

Im Jahre 1945 habe ich (in einer Pa- 
tentanmeldung) den Vorschlag gemacht, 
den EinfluB der exzentrischen Verlagerung 
beim Herstellungsvorgang der Feder da- 


gangsmaterial ein Rohr verwendet, das 
im entgegengesetzten Sinne ausmittig ge- 
bohrt ist. Noch ginstiger ist die Wirkung, 
wenn man die Gegenexzentrizitat gréBer 
wahlt als zum Ausgleich erforderlich ist. 
Man verbindet so die Wirkung des Hoch- 
druck-Bourdonrohres mit der des aus- 
mittig gebobrten MeBrohres, fiir das ich 
in einer friheren Arbeit! Berechnungs- 
grundlagen entwickelt habe. 


8. Zusammenfassung. Zur Messung 
hoher Driicke 1m Bereich von 100 bis 
10000 at werden in der Technik vielfach 
dickwandige Roéhrenfedern nach Art des 

Abb. 19. Kriimmungsanderung der Kreisrohrfeder ohne und mit der IIR A peu ee Ba 
pes: EPS RET ares are dern werden durch Flachschlagen eines 
Kreisrohres und anschlieBendes Kriimmen 
hergestellt. Die Theorie derartiger Hochdruckfedern wird auf einen ebenen Formanderungszustand 
zuriickgefihrt und die Spannungen und Deformationen werden berechnet, wobei die Querschnitts- 
form durch ein rechteckiges Mittelstiick und zwei Randbégen (Flachbogenquerschnitt) idealisiert 
wird. Als Sonderfall wird auch der kreisformige Querschnitt behandelt. Bei kreisférmiger oder 
kreisahnlicher Querschnittsform werden verschiedene Korrekturgréfen von wesentlichem EinfluB, 
die ebenfalls abgeleitet werden. AbschlieBend wird auch das Bewegungsgesetz solcher MeBfedern 
untersucht. 

Die hier dargelegte Theorie der Hochdruckfeder habe ich mit weitergehenden Vereinfachungen 
und ohne die Exzentrizitats- und Langungskorrektur bereits 1943 entwickelt. Ergebnisse wurden 
auszugsweise in einem Aufsatz in ,,Die Technik“ 3 (1948) S. 23—30 wiedergegeben. In der vor- 
hegenden Neufassung sind alle Tabellen neu berechnet worden, wobei ich von Frau Otto in Gottingen, 
wesentlich unterstiitzt wurde. 


0 Of 02 03 o¢ O05 06 “O7 


(Eingegangen am 3. August 1955.) 
Anschrift des Verfassers: Dr. W. Wuest, Gottingen (Hann.), Béttingerstr. 6/8. 


1 W. Wuest, Ing.-Arch. 19 (1951), S. 12. 


durch zu beseitigen, daB man als Aus- — 
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Kritische Drehzahlen gewisser Rotorformen unter Beriicksichtigung 
der Kreiselwirkung bd 


Von L. Fehrle 


1. Einleitung. Die Berechnung der kritischen Drehzahlen von endseitig zweifach kugelig ge- 
lagerten Wellen mit konstanter Biegesteifigkeit, die ganz oder teilweise kontinuierlich mit Scheiben 
von. linear veranderlichem Halbmesser besetzt sind, fiihrt auf ein technisches Eigenwertproblem, 
dessen Auflésung im allgemeinen unter Verwendung von Naherungsmethoden erfolgt, da die strenge 
Lésung, soweit ihre Existenz mathematisch gesichert ist, meist nur durch umfangreiche Rechen- 
arbeit zu gewinnen ist. 

Ausgehend von dem durch eine Differentialgleichung einschlieBlich Rand- und Ubergangs- 
bedingungen gegebenen Eigenwertproblem — der Eigenwert ist proportional dem Quadrat der 
kritischen Drehschnelle — hat K. Karas’ die angenaherte Auflésung mit der Methode von W. Ritz? 
durchgefiihrt. Diese Methode liefert, was sich unter bestimmten Voraussetzungen exakt beweisen 
la8t, obere Schranken fiir die positiven, untere Schranken fiir die negativen Eigenwerte. Von 
praktischer Bedeutung sind lediglich die positiven Eigenwerte; das Auftreten negativer Eigenwerte 
bei beriicksichtigter Kreiselwirkung im Falle synchroner Prazession im Gleichlauf besagt, da die 
Kreiselwirkung der Scheiben fiir w>@ (@ entspricht dem gréBten positiven Eigenwert) ein kri- 
tisches Verhalten der Welle verhindert. 

In der vorliegenden Arbeit erfolgt die angenaherte Berechnung der kritischen Drehzahlen unter 
Beriicksichtigung der Kreiselwirkung mit einer Methode, die, ankniipfend an ein von R. Grammel?® 
entwickeltes Verfahren, die Aquivalenz zwischen dem durch die Differentialgleichung samt Rand- 
und Ubergangsbedingungen und dem durch eine Integrodifferentialgleichung gegebenen Eigen- 
wertproblem benutzt. Das Verfahren von R. Grammel, mit welchem sich die kritischen Drehzahlen 
bei unberiicksichtigter Kreiselwirkung berechnen lassen, liefert ebenfalls obere Schranken; unter 
gewissen Voraussetzungen ist streng beweisbar, daB diese genauer sind als die mit dem gleichen 
Naherungsansatz (s. u.) errechneten entsprechenden Ritzschen Lésungen. Diese charakteristische 
Eigenschaft der Grammelschen Methode tritt in ihrer hier durchgefiihrten Verallgemeinerung auch 
bei beriicksichtigter Kreiselwirkung noch im Falle synchroner Prazession im Gegenlauf auf. Bei 
synchroner Prazession im Gleichlauf allerdings gelangt man zu Naherungen, die z.T. oberhalb, 
z.T. unterhalb der strengen Lésungen liegen kénnen. 

Ausgangspunkt fiir die Berechnung der kritischen Drehzahlen, deren Ursache die nie ganz zu 
vermeidende exzentrische Lage der Scheibenschwerpunkte ist, ist das nachstenende Aquivalenz- 
prinzip?: 

,,Man berechnet die kritischen Drehzahlen als diejenigen, mit denen die Welle in ausgelenktem 
Zustand stationar um die Ruhelage umlaufen kénnte, wenn die Scheibenschwerpunkte genau mit 
der Wellenachse zusammenfielen.” 

Die mathematische Erfassung des durch diese Vorschrift charakterisierten Problems kann auf 
verschiedene Weise erfolgen. Die Rotationssymmetrie des stationaren Umlaufs der ausgelenkten 
Welle (die Durchsenkung infolge des Eigengewichtes des Rotors bleibt unberiicksichtigt) legt es 
nahe, sich der Methoden der Statik zu bedienen und die Wellenauslenkung im kritischen Zustand 
als eine Folge gewisser duBerer Krafte anzusehen. Diese Uberlegungen fiihren auf eine Integral- 
bzw. Integrodifferentialgleichung. Eine andere Ansatzméglichkeit ist durch das Hamiltonsche 
Prinzip gegeben; bendtigt werden hierzu potentielle und kinetische Energie des aus Welle und 
Scheiben bestehenden Systems beim Roti¢ren mit der kritischen Drehschnelle. Dieser Weg fiihrt 
auf das durch eine Differentialgleichung samt Rand- und Ubergangsbedingungen definierte Eigen- 


wertproblem. 


* Auszug aus der Frankfurter Dissertation. 1. Referent Prof. Dr.-Ing. K. Karas, 2. Referent Prof. Dr. 
R. Moufang. oi 

1 K. Karas, Die kritischen Drehzahlen wichtiger Rotorformen. Wien 1935. 

2 W. Ritz, Crelles Journ. 135 (1909), S. 1. 

3 R. Grammel, Ing.-Arch. 10 (1939), S. 35. J 

4 R. Grammel, Neuere Untersuchungen iiber kritische Zustinde rasch umlaufender Wellen. Ergebn. d. 
exakten Naturwissensch., Bd. 1, S. 92. Berlin 1922. 
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Als auGBere Krafte treten in erster Linie Fliehkrafte und Kreiselmomente in Erscheinung, tiber 
die kurz folgendes gesagt sei: 

a) Fliehkraft: Bezeichnet man mit y(x) die wegen der Rotationssymmetrie des Problems als — 
ebene Kurve dargestellte Wellenauslenkung im kritischen Zustand, mit w die kritische Dreh- | 
schnelle, so tritt an einer im Punkt x, befindlichen Scheibe der Masse m eine Fliehkraft F=m (Xo) w 
auf, 

b) Kreiselmoment: Im ausgelenkten Zustand hat im allgemeinen jede Scheibe die Ebene, in 
der sie im Ruhezustand lag, verlassen und fiihrt eine erzwungene regulare Prazessionshewegung 
aus, deren Kreiselmoment, mit dem die Scheibe beim Umlaufen auf die Welle wirkt, zu bestimmen | 
ist. Unter der Voraussetzung, daB die Scheibendicke gegeniiber den anderen Scheibenabmessungen | 
vernachlissigt werden kann und die Wellenauslenkung klein ist, liefert die Kreiseltheorie fiir die | 
hier allein wichtigen Falle der synchronen Prazession im Gleich- und Gegenlauf die Kreiselmomente | 


Ke = Oay'(%) @?, Ken = —3 Oay' (xp) @? = — 3 Ke 

(s =synchron, Gl = Gleichlauf, Gn = Gegenlauf). In diesen Gleichungen, die fiir eine in x, befindliche | 

Scheibe gelten, bedeutet QO, das aquatoriale Massentragheitsmoment der Scheibe und y'(x) die 

erste Ableitung der Auslenkungsfunktion. 
Geht man zu einer kontinuierlichen Scheibenverteilung iiber, so wirken auf ein mit Scheiben 

besetztes Wellenelement der Lange dx die Fliehkraft dF = uz 1°(x) y(x)w*dx und das Kreisel- 


moment dK¢,— 2 wz r(x) y'(x) w2dx baw. dKe, = — - pe a r4(x) y’(x) w* dx, wobei w die 


Masse des Scheibenmaterials pro Volumeneinheit ist. Dabei wird angenommen, daf sich die auf 
der Welle sitzenden Scheiben trotz dichter Aneinanderreihung nirgends beriihren mégen. 


2. Die Integrodifferentialgleichung der zweifach kugelig gelagerten Welle. Gegeben sei eine an 
heiden Enden kugelig gelagerte Welle der Lange | mit einer durch r= r(x) in 0<x<I definierten 
kontinuierlichen Scheibenverteilung. Vorausgesetzt wird: 

a) die Biegesteifigkeit der Welle, EJ, ist konstant; 

b) die Masse der Welle ist vernachlassigbar gegeniiber der Scheibenmasse. 

Zur Herleitung der Integrodifferentialgleichung! wird die gemaB Abb.1 in x==t durch eine 


Fliehkraft dF = 27r(t) y(t)w?dt und ein Kreiselmoment dK¢, = 2 a ri(t) y'(t) w? dt belastete 


Welle betrachtet. Aufzusuchen sind die Auslenkungsfunktionen y,(x,t) in 0 <x <t und y,(x,t) in 
t<« <I. In den Lagern ist die Wellenauslenkung gleich Null, d.h. y, und y, miissen den geome- 


[ Feld Feld 2 


trischen Randbedingungen 


¥(0,t) = y(I, t) = 0 (la) 
genigen. Ferner verschwinden dort 
wegen der vorausgesetzten freien Lage- 
rung die Biegemomente, was zum Aus- 
druck kommt durch die dynamischen 
Randbedingungen 


Abb. 1. (Ol Dea YALE) ee (1b) 


Zur Bestimmung von y, und y, ersetzt man voriibergehend dF und dK¢, durch die endlichen 
Gréfen F und K und berechnet die Biegemomente M,(x) und M,(x) fiir die Felder 1 und 2. Nach 
kurzer Rechnung ergibt sich 


M(x) =>) 2+ x, (2a) 

My(x) = (lea) p= =a), (2b) 
Mit £¥ — — FO folgt aus (2a) und (2b) 

OU) = Fl) a (3a) 

eS ait) + ail. (3b) 


') Siehe FuBnote 1 von S. 111. 
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Die Gleichungen (3a) und (3b) erfiillen die Bedingungen (1b). Durch zweifache Integration von 
(3a) und (3b) nach « erhalt man die gesuchten Funktionen y, und y,. Die dabei auftretenden vier 
Integrationskonstanten sind so zu bestimmen, dai (la) erfiillt wird, ferner y,(t,t) = Yo(t,t) und 
y,(t,t) =y;(t,t) ist; die beiden letzten Bedingungen sichern den sprung- und knicklosen Verlauf der 
Auslenkungsfunktionen in x=t. Nach Ausfiihrung dieser Rechnungen erhalt man 


Fx 


Yi(%t) = Faq (4) (2 Ut t? — x?) sen +2*2+3«£2—6xtl), 


K 3 
| éEyi * 
Fiihrt man wieder statt F und K die Elementarbelastungen dF und dK; ein und beriicksichtigt, 


daB t alle Werte zwischen 0 und / annehmen kann, so findet man durch Superposition aller Einzel- 
auslenkungen 


Yo(%>t) =a») (2 lx — x? — 7?) 31+ 3x20 +2 x12?—32)), 


y(«) =p wo | | om (l— x) t (21x — a2 — 2) 7°) y(t) dt 


; (4) 
—- btw | [ aE (?@—3 281+ 3xe2+ 2x2 <3 21) r4(t) y'(t) dt 


0 


1 Z F 
| o4Bg1 aF Axe —6xrhro yal, 


d dij d*y 1 dy uaa? ls 
Mit eben yori | fate, rae 10, a ae ig fae? tay (5) 


erfolgt der Ubergang zu dimensionslosen Veranderlichen; damit geht (4) iiber in 


[ 4 r@E—P— a) (0) (2) dr 


i [ ¢O—NERr— PB) ax 


(& —3& + 3&r? + 2& —3 7?) of(1) n'(x) dr 


Cx | a (& + 2€ +3 2?— 6 &r) oX(t) y(t) i 


Nach Einfiihrung der Kerne K(é,t) und K,(&,7) mit 


K(é,t) =~ (I—7)é Qt — 2 —®&) = Ky (G7) fiir E< tr, 
=4 (18) 1 2&2 #— 1) = Ky(67) fiir E> 13 F 
- ( 
K,(é,1) = 5 ( + 26 + 3€0* — 6 £7) fiir €<r, 
=, (8 $e 43 é cP 4 2 F — 3 72) fir ES 7 
folgt aus (6) : 
n(é) =u / K(E,2) @°(t) n(z) dt — | K(E,7) @X(t) n/t) dr). (8a) 
LO) 0 
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Dies ist eine lineare homogene Integrodifferentialgleichung mit dem symmetrischen Kern K(é,7) 
und dem unsymmetrischen Kern K,(é,7). Der Kern K(é,r) wird auch »»Musterkern® genannt; er 
gibt die Auslenkung an, die an einer Welle von der Biegesteifigkeit 1 und der Lange 1 durch eine 
an der Stelle r(0<1<1) senkrecht zur Ruhelage der Welle wirkenden Kraft der Grofe 1 hervor- | 
gerufen wird. Der erste Term in (8a) stellt den von der Fliehkraftbelastung herrithrenden Anteil 
der Auslenkung dar, dem sich die von den Kreiselmomenten verursachte zusatzliche Auslenkung 
(zweiter Term) iiberlagert. Aus der fiir synchrone Prazession im Gleichlauf abgeleiteten Glei- | 
chung (8a) erhalt man sofort die entsprechende Gleichung fiir synchrone Prazession im Gegenlauf, | 


wenn man wegen Ke, =— 3 Ke, den zweiten Term von (8a) mit — 3 multipliziert: 


1 Leta: | 

(8) =u || K(Ex) 922) nl) dr +3. Ky(Esr) oF(e) 1!) de (8b) 
0 0 

Ohne Kreiselwirkung ist K,(&,7) = 9, also 


1 


= u [K(E,r) 0°(t) (7) dr. (8) | 


Zwischen K(é,r) und K,(é,7) bestehen verschiedene Beziehungen!; hier sei nur eine erwahnt, welche 
die Elimination des Kernes K,(é,7) aus (8a) und (8b) ermoglicht. Es ist 


(1 
oh — (& + 2& + 3&2? — 6&2) | fst | 
EEO 4K, é,1) _[s fiir (9) | 
; ([F@—se + see tze oa) fee ee 


n(é) = | [ Ken) eX) ale) de —] | ALP 0F(2) of a 


0 


und nach Ausfiihrung einer partiellen Integration weiter zu 


— 


1 1 
1 pn hes ih if d ' 
n(g) =u | [KG nee) nle) de — 7 KE 2) OTe) S| [KE F (0M) (0) dr 
} 6 
1 (10a) 
uw} { (Ez) {0%(r) n(x) + 4% (oe) 91) de — 1 K(E,2) ot) 9!) "> 
j ’ 8 1) 4 dt Q / { A ? ls Ui =O ms 
é 
J 
Das Glied — = K(&, 7) o4(t) 4'(t) Ie verschwindet im allgemeinen nicht; ausschlaggebend fir — 
. t™=0 5 : 
die einzusetzenden Zahlenwerte der oberen und 
x-Feld A-Feld h-Feld unteren Grenze ist der Definitionsbereich von (6), 


d.h. Anfangs- und Endpunkte der Intervalle, in 
welchen 0(&) nicht identisch Null ist. : 

Im allgemeinen liegt der Funktion 0() ein line- 
ares Bildungsgesetz zugrunde (Abb.2), das sich mit | 


tgo) =i durch 0(f) =o, +ié (aS éS fP) 
aasdricken 1aBt. 


Gleichung (10a) macht fiir den in Abb.2 dar- 
gestellten Rotor unter Beachtung der Bedeutung 
des Kernes K(&,r) folgende Aussage: | 
Die Auslenkung 7(€) stellt sich ein als Folge einer kontinuierlichen Belastung 


Abb, 2. 


(é) = « |o*E)'€) + 7 Fe (et n'®)] im «SES 


1 Siche FuBnote 1 von S.111. 
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und zweier Krafte (senkrecht zur Ruhelage der Welle gerichtet) 


1 A j i 
Py =gugla)n(a) inE=a, Py = —4 woK(B)n'(B) in E=8. 
In diesem speziellen Fall ist also 
I , T= 1 t= fp 
— = Kx) of) n'@) [= =— | K(E0) ofl) n' [7 ". 


, 


Zwischen P,, Pz und dem Gewicht G der in x <  < £ kontinuierlich verteilten Last — x ul o£) n'()) 
besteht folgende Beziehung : 


B 
1 / / 1 ii ) 1 y 
G =u | (MG) 1G) a = woXB)nl(B) ~ | uoKa)n'(@) =—(Pp + Pa), 
d 
d.h. die Summe der Vertikalkomponenten der Kreiselmomentbelastung ist Null, was ja von vorn- 


herein nicht anders zu erwarten ist. 


Mit (9) folgt aus (8b) analog 


ne) =u] [KE olen — | gleterrte| ae + | KEM otervten 
0 


o{ (LOb) 
Fiir einen Rotor gemaB Abb. 2 gilt jetzt 
a(8) =u (0X8) nl8) —F 3 (0 1'®) ino SEXP, 
P, =— ugha) 7a) inE=a, Py=+uo(A)n(8) inf =f. 
3. Variationsproblem und Differentialgleichung. Ausgehend vom Hamiltonschen Prinzip 
6H = fv —L) dt =0 (t = Zeit) , (11) 


wobei V die potentielle, L die kinetische Energie des betrachteten Systems ist, gelangt man zu 
dem durch eine Differentialgleichung samt Rand- und Ubergangsbedingungen definierten Eigen- 
wertproblem. 
Fiir die im kritischen Zustand rotierende Welle ist im Falle synchroner Prazession im Gleichlauf 
l l 


= ar mw? [ r2(x) y?(x) dx + = ae | o ied 4) ee r4(x) vx) dx , (12) 

0 Ce 

I 
1 a, 
er eee (13) 
é 

Mit (12) und , folgt aus (11) nach Division durch oe 

0) i i BIy(x x)— pnw (Pio y?(x) + ; r4(x) ; Tax) ya) dx dt = 0". (14) 


1 
Da y(x) nebst allen Ableitungen nach x nicht von der Zeit abhangt, ferner der Term — yw 7 w)* > x M(x) 
die Variation nicht beeinfluBt, reduziert sich (14) auf 
l 
1 


9 | [BIy’%e) —m 08 (2) 4) — Ge) yAA))| dx = 0. (15) 


Vermége der Transformationsgleichungen (5) geht (15) ttber in 
1 


| (x Me epee 7 ue" 7/*—ue *nf\ ae = 0. (16a) 


0 
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Fiir synchrone Prazession im Gegenlauf und unberiicksichtigte Kreiselwirkung erhalt man ent- 


sprechend 1 
d | (n'" — + up'n’? —ug? n*) dé = 0, ey 
0 
i 
6 f(y!” — uo"7?) dé = 0. (16) 
0 
Eine Auflésung dieser Variationsprobleme kann iiber die entsprechenden Euler-Lagrangeschen 
Differentialgleichungen 
n\*) pee uo" 7 = 0, (17) 
Mea Ns ee Nath 17 
" z Ulo*n') — uo?n =9, (17a) 
fhe + u(o*n')’ —ug?n = 9 (17b) 


erfolgen. 
Die Weiterfiihrung der Uberlegungen erfolgt fiir den speziellen Fall eines Rotors gemaB Abb. 2, 


fiir welchen 0, =0,=9, @,5£0 gilt. An die Stelle der Gleichung (17a) treten fiir synchrone Pra- 
zession im Gleichlauf die Differentialgleichungen 


yi =0 (J =%M) >| 
ny) — = wei mi)’ — vei m = 0, | 
da in jedem der drei Felder ein anderes Bildungsgesetz fiir @ vorliegt, folglich die Auslenkungs- 


funktion 7 in die drei Komponenten 7,, in 0<& <a, n, maSéLPf und y, inf SEX) zerfallt. 
Die Lésungen 7, und 7,, von (18) haben den bereits bekannten Bedingungen fiir die Lager 


ml0) = (0) = 71) = nut) = 0 (19) 
zu geniigen. Ferner sind 7, und 7, baw. 7, und 7, an der Stelle « baw. / durch je vier Ubergangs- 
bedingungen miteinander verkniipft. Je zwei Ubergangshedingungen fiir ~ und § ergeben sich 
aus der Tatsache, da die Auslenkungsfunktion tiberall sprung- und knicklos verlauft, die Funk- 
tionen 7,, 4, und 7, also die Bedingungen 


NulX) — Nad)» nA(B) a HAP)» | (20) 

no) = n2(a); m(P) =nAB) J 
erfiillen miissen. Die Gleichungen (19) und (20) gelten auch fiir synchrone Prazession im Gegenlauf 
und unberiicksichtigte Kreiselwirkung. Die restlichen Ubergangsbedingungen lassen sich aus den 
entsprechenden Variationsgleichungen gewinnen. Dies erfolgt hier fiir einen Rotor gemaB® Abb.2 
im Falle synchroner Prazession im Gleichlauf. Gleichung (16a) fiihrt zu 


(18) 


(04 
|| nie d& + j (ni 42 +> ue! na —u 0; nt) dé + | Nu we =U (21) 
0 
Nach Ausfiihrung der veces und Division durch 2 folgt heat 
sae B 1 
417 Qo fi 7 " il U7 , uv " 
| me onkde + (7 int + = wot ni Oni —uel mons] dé + | ni onde =0. (22) 
) x B 


Durch partielle mia lassen sich folgende Glieder dieser Gleichung umformen: 


4st 


a 
I ne One dé = yy Ons) — mx’ One|, + Jn neds. 


B B 
fe Piel It uv 1B AAA r < 
J ni oni dé = nf dni\’ — ni’ ona |, + Jn? dm, a8, 
a a 
(23) 


, he Ae uv ae Yh 1 : 
[6 dn dE = rip Oral, — Nhe Nal [mi On dE 


B 
- wjo nh Oni dé =~ w ob nf On|? — uf (ohni)' bmp a. 
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Beriicksichtigt man hierbei die Rand- und Ubergangsbedingungen (19) und (20) und bedenkt 
dafs diese weiterhin in jeder beliebigen Nachbarlage nur dann erfiillt bleiben, wenn die Waristoten 


bn;, On; (J =A, mu) den Bedingungen 
6,0) = dn, (1) = 0, 
010) = dn,(x) = dna), dnAlB) = dn,(B) = dn(8). 
bn dx) = Oya(a) = bn'(x), (8) = dn, (B) = dn'(B) 
geniigen, so folgt damit aus (23) 


[u: On dE = n(x) dn’(ax) — v(x) > dn(o) + fn dn, dé, 
0 

B 

Ju Oni dE = 73(B) + dn'(B) — 3x) > dn'(a) 


$s m P 
— nx (B)- dn(B) + 7"(a)- dn(a) + J n$? dn, dé, 


Be u v / 4th iY 24. 

J cdi A& = —nftB)- b1f(8) + rhe(B)- On(8) + [rf Ong A sy 
f 

oes eer 

uf hn dnt a& = —; wol(B) ni(B) dn(8) 


a 


1 pact iis 
mide 03(a) n3(x) O(a) Sark uf (02.2)' On, dé. 


a 


ast 4st ' 1 
— dn(a) |e) — nf" (0) + | w ob(a) ni(o) 


— dn(8) Pi") —1he'(B) #0118) 418) 25) 


B : 1 


Mit (24) ergibt sich aus (22) nach Zusammenfassung der mit d(x), d7'(a), dy(8) und dn'(f) be- 
[ 1 ry) F : f . = 
i ie u(92. 2)’ — u 9; | dn, dé +- | my’ On, dé = 0. 


hafteten GHieder 6 
dna) [oelo) — nfla)] + 91/(B) [lB) — XB) ; 
eee 
‘ | mon, dé + | 
0 By es ; | 


Beachtet man jetzt, da®B (25) fiir voneinander unabhingige Variationen 67;, dj (j =x, A, mu) erfiillt 
sein inuB, wobei diese den oben genannten Bedingungen geniigen, so erkennt man, da die Inte- 

_granden der drei Integrale fiir sich verschwinden miissen — womit man wieder zu (18) gelangt —. 
und da®B ferner jede der eckigen Klammern gleich Null sein mu. Die letzte Aussage liefert die 
noch fehlenden Ubergangsbedingungen fiir synchrone Prazession im Gleichlauf: 


nila) = nla) nle"(a) = ni (0) — | wolla) ni(0) | 
! 26 

u uv ‘tf 407 J 4 / ( a) 

WlB) =nAB)» Nu (B) =a (6) — 4 & Ox{B) MAP). | 


_ Fur synchrone Prazession im Gegenlauf baw. unberiicksichtigte Kreiselwirkung erhalt man analog 


uM uv 44s 4/7 3 / 
milan) = nila) «9h (0) = iC) +g oAa) nO) | 


a “ 5 ; (26b) 
nil) = nb). CAB) = nie (B) + > wohl) lB) | 
bzw. 
AO). =O) on = Ne (0) = na (%)> | (26) 


ni(B) = nilB)»  2'(B) =n (B)- J 
Aus den Gleichungen (26), (26a) und (26b), die fiir einen Rotor gemi® Abb.2 abgeleitet wurden, 
lassen sich sofort die entsprechenden Gleichungen fiir den allgemeinsten Typ der Mehrfelderwelle 


angeben. 
9 
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Durch die Differentialgleichung (17) baw. (17a) baw. (17b) einschlieBlich der zugehérigen Rand- 
und Ubergangsbedingungen einerseits und die Integralgleichung (8) bzw. die Integrodifferential- 
gleichung (10a) baw. (10b) andererseits wird das gleiche Eigenwertproblem dargestellt. Diese 
Aquivalenz kann man kontrollieren, indem man gemiB (8) bzw. (10a) baw. (10b) fiir die einzelnen 
Felder die Auslenkungsfunktionen aufstellt, diese viermal nach é differentiiert und die Werte der 
in den entsprechenden Rand- und lilereanthediacumeen vorkommenden Funktionen und ihrer 
Ableitungen berechnet. 


4, Eine Verallgemeinerung des Verfahrens von R. Grammel zur angeniherten Berechnung der 
kritischen Drehzahlen unter Beriicksichtigung der Kreiselwirkung. Bei der angenaherten Auflésung 
von Eigenwertproblemen spielen die sog. ,,zulassigen Funktionen“ eine wesentliche Rolle; sie sind 
fiir die hier zu behandelnden Falle wie folgt definiert: 


Eine reelle Funktion v(&) hei®t zulassige Funktion des vorliegenden Eigenwertsproblems, 
wenn v(&) viermal stetig differentiierbar ist und den Randbedingungen v(0) = v(1) =v'’(0) =v’’(1) = 0 
genugt. 


Endlich viele linear unabhingige zulassige Funktionen v,== 0 bilden in der Form 
P(E) = S anve(€) 
kad 


mit zunachst noch unbestimmten Koeffizienten a, einen ,,Naherungsansatz‘‘. 


Ausgangspunkt der Naherungsrechnung (wieder fiir den speziellen Fall eines Rotors gemaf 
Abb.2 fiir synchrone Prazession im Gleichlauf) ist die Variationsgleichung (22), die sich mit (23) 
umformen lat zu 


aX 


[nb dn,d& + font), — ni" Ons, 


0 


i 1 Oxo 2 n i ate) TIL-«& iB it Ete 
# I (m =a! (02m) — vein) Oma d& + 1. 0a|, — 2 ONa|y + | woin dnl. pag 


417 ® il 
=f fil mdm AE + mp Mg — The ON |, = 9- 


Die ae ee Auflésung des Eigenwertproblems vollzieht sich folgendermafen: 


1) Man ersetzt in (27) die Funktionen 7;(j =x, A, 4) durch solche Iterationen, die sich auf Grund 
der Vorschrift 


B 
2 | Ky (é,7) 103 


a“ 


(rot) +4 < (ofle)¢ ier, 


T(E) =— KeslE, B) | 0X8) 9'(B) + Koalé a) + of(a) g'(0 


+ | Rute nfe arp) +44 7, (At) 


nid) = —KoslEsB) OSB) ¢'(B) + Kor (Ea) + ofa) ¢"( _ 
(28) 


a i Ko, (&, 7) eile) Q(t) + 2 ve (03(z) : dt. 
nAE) = — Ko2(E, B) - ex(B) p(B) + Koo(é, o) = o8(a) (a) 


B 
+ | Kos(é2) leit) gt) +> < (eke) 9") 


a 


— 
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und der dazu aquivalenten Vorschrift 


ne (yO, | 
4 ered peed 
a(é) = oAl€)  (€) + | Ge (X(€) @'(E)), (29) 
nw(é) = 0 | 
einschlieBlich der Rand- und Ubergangsbedingungen 
Ti(0) = HL(0) = 7,(1) = (1) = 0. 
(or) == ,(a) » nip) = oe; 7AP)- | 
Hilo) = 7(0), 73(B) = 7,(B). 
—/ —4 u" ( (30) 
Nx(%) = 173(0), vil (8) = 7.(B) - 
bend 4 borat Wl —f fp lot Waid il i 
Na (& (a Ves UP (0) — oi(a) @" (a) 0 "Nu () = an (B) —= @4(B) Y (B) | 
J 
darstellen lassen, wobei y = ) ay %, = a ay, V, ist. 
haa = 
2) Man ersetzt die Variationen 6y;, Onj (j =x, A, w) durch Variationen von der F orm 
ie [osese. | 
AG) =k . 
oy (€) = Pa bare (N=) én lati se ep) 6 = 0,13 (31) 
li sx L ee, | 
lb 
hierbei gilt also 
| oi.(0) = 0F1,(1) = 0, | 
On,(o%) = O(a) = O(a), — 07,(8) = 07,(6) = On(B), (32) 


d7,Aa) = O72(%) = Om (x), O18) = On7,(B) = On (P)- | 
3) Der durch diese Vorschriften entstandene Ausdruck wird gleich Null gesetzt. 
Geht man mit den unter 1) und 2) gemachten Vorschriften in (27) ein, so folgt unter Beachtung 
von (30) und (32) zunachst fiir die Glieder vom Typ 177, 677, i usw. — wobei jetzt u statt wu geschrie- 
ben wird, um hervorzuheben, da es sich um eine ‘pine lepaar handelt — 
ppp Nh PP s/t 1 =r 1); — |p 
"x ON, e emt fog OM |, =i "4 oii, |. = na ’ oy. He Ne 57, |. Ig fle On, |, “ A U0, On; Hie 
= O01] (a) [7:(o) — 112(a)] + 07° (B) [n2(B) — 7,8) 
Ca Sis de 1 = mart o= SPF. ALLA = maw 
— éii(a) li (0) — f(a) + 7H ofa) ms(0)| + o7(6) |7ic(8) —7i"(B) + 7% 0X(8) THA) 


— Ife (0) O7,(0) + Me (0) 07,40) + MAL) OL) = Mh (1) ONj,A1) 


= — 00) |e) 78") + 4 Hobo) riled] + dm(B) |B) — HB) + 5% 088) THB) 
= — 6n(0) |— offa) o'(o) + | # ofa) Hila)| + O78 )— 5 ett8) (P+ 4 HMB) THA) 


Weiter verschwinden in (27) wegen der ersten und dritten Gleichung (29) die Integrale (encund 
; 0 


1 
[.... Man erhalt somit zusammenfassend und gemaB Vorschrift 3) das Ergebnis gleich Null 


setzend 
B 
[[2oe@ +5 4 (GOO) — FH REO HE) —# AO HO] aE 
— dii(a) | — + 00) (a) + | HX) Halm)| + 9718) ae of(B) o'(B) + 4 #eh(B)AAB)| = 0. (33) 
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Mit (31) geht (33) wber in 


woraus wegen der Willkiirlichkeit der da; folgt 


an Ded ee Lit dpe ete ae ae 
[ete oe + ag oH #) — i |e (0B) HB) + CHE) nn(8)|| vl) a8 | a 
— + foi(8) 06) 9 OF +t [oA CHE) THO], =O, = 1,-- 4 7). 
Mit den Umformungen 
— 7% [ le) ge (HO HE) dé =— Fi [vl QE) HOY. + | 4 (HE) HE) MAN 
‘ a B 
— i | Tale) ge loH(6) 08) a8 
ye ee ee 1 : 1 f 
zl vil£) ge (OME) gE) 4& = vil) OnE) @(E — | vi(E) o2(€) y'() dé 


folgt aus (34) 


| v8) gel) #8) — 4 ME) C8) 0 — 4m) | Ge (GE) MMB) + GE) HB] aE 
—i|— ofa) ol8) 8) = 0 @=1,-+4 2). 


ja 
Setzt man fiir 7, p und gy’ die entsprechenden Werte ein, so ergibt sich aus der letzten Gleichung 
n 


nach Vorziehen von |)’ 


(o5(€) vi(@)) + oH(€) v8) | 


(35) 


= . paper 01(B) v;(B) (= Ko (6, B) = 03(B) v z(B) )+ Ko2(P, a) 4 a8 03(a) vy (cx) 


+ | Kol 79} ofl) eale) + Fe (oll) vee) dr] as 


—EEE~ 


1 


+ | Koa(Be1) [eter ete) + ge (ete sk(n)| ar) 
+ oh(a) vf(«) (= Kx(a, B) > @h(B) i(B) + Koala. a) + ofa) ve(a) 


2 i Kon(0s 1) [@i(2) wale) +p ge (0h) vi(e)] dr) lo, 
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mit v,(€) =— Kyl, B) = of(6) v}(B) + Koo(é, a) = o3(a) vj(x) 


+ | Kote) [oily oe) + 5 (ote ass) ar 


+ | KolEs1)/ eX) oe) + geen) ofp) de 
PL =F ella) of(a), Pt; =—-+ off) vif). 


(ein Stern soll auf synchrone Prazession im Gleichlauf, zwei Sterne auf synchrone Prazession im 


Gegenlauf hinweisen) geht (35) iiber in 
n B 
>t | [ (EO 2@ n@ — FeO HO Mw A inl) a 


k= 
— u( P¥;: 0,,(8) + Pr: 7.) Si (haa bee aie (36) 


Hier ist v;,(£) die Wellenauslenkung ina < £ < £, hervorgerufen von der ina < é< f kontinuier- 
lich verteilten Last q;(&) und den in € = « und € = 6 senkrecht zur Ruhelage der Welle angreifenden 
Einzelkraften P;, und P},. Demnach ist 


t =| at) Bald) dé + Ph O(B) + PET ,(a) (37) 


die von der aus qj(), P3; und P}; bestehenden Belastung bei der Auslenkung v,;,(&) geleistete 
doppelte Arbeit. Mit (37) und 


Ai, ef et dé , Kiel Gh eco tale (38) 
folgt aus (36) : ‘ 
| > a (40-5 Kin —it Bh) = 0 hee ae (39) 
Durch Nullsetzen 8 ee dieses linearen homogenen Gleichungssystems 
fiir die a; erhalt man eine Gleichung n-ten Grades fiir u: 


| 4i,— > Ki, — i re Oi best. en) (40a) 


Die Auswertung der Gréen B?,, verursacht schon bei einfachen Beispielen sehr viel Rechenarbeit. 
Diese 1a8t sich, wie R. Grammel! gezeigt hat, unter Verwendung der Biegemomente umgehen. Ist 
Mj = M (PX), P3;, q{) das zur Belastung Px;, P3;,q; gehérige Biegemoment bei reiner Biegung, 
so ist 


] 
+, = | MP Mi dé. (41a) 
0 


Fiir synchrone Prazession im Gegenlauf und unberiicksichtigte Kreiselwirkung gelten die vor- 
stehenden Ausfiihrungen in entsprechender Form. Die sich einstellenden Endergebnisse lauten: 


Ain += Ky, —% Bit | =0, (40b) 
A;, —U Bix = 9, (40) 
1 
oye Mus de, (41 b) 
0 


1 Siehe Fufnote 3 VOTUES aelalelle 


mit 


bzw. 


Fehrle: Kritische Drehzablen gewisser Rotorformen Ingenieur-Archiv 


Ba mM My, dé . (41) 
Mj* = M(P, PRK gf"), af") = ek) (6) — > ae (HE) HE) 
Prt =—> ola) ofa), PEF = oh) 24(8) 


M;=M(qj), —-aj(€) = @2(€) o(€) - 


eee durch 9(6) = Gop +ipé in % = ESB (p—1,...m) (GX Pi Sutr— pia 
0 < a1; bm < 1) gegebenen Scheibenverteilung erhalt man ebeatalle die Gleichungen (40) bis (40b). 
Dabei bedeuten jetzt 


An = D> uae Kee eres 


p=1%p p=1 %p 


; | 
Bi == | MPMe di Ms MPs Ge Ps eee ora ee 
0 


1 
Bi = [MEME dg, Mf* = M(P&j, PRG. «5 PELys PRs ahs --+ aad) 


i = 1 Oop) vi(ap), PE» =—q ob Bp) (Bp) 
aplE) = ob) (2) +7 al) HO) — (Hp SESH), 
Pi = + of(B,) Vi(Bp)> Par aes O(Xp) vi(ap), (p=1,--» m), 
GE) = e8 YO-FR(GOYO) OSESH) G=L--+ 0). 


5. Eine Gegeniiberstellung der Verfahren von W. Ritz und R. Grammel. Nach der Methode von 
W. Ritz’ wird das Integral 


l 
If == [ (a8 + metn't—u orn) dé 
0 


(synchrone Prazession im Gleichlauf) zu einem Extremum gemacht, indem man hierin J) durch 


n * 
> a vKI) (j = 0, 1, 2) ersetzt und die Operation a =) (t= 1... ., n)taustulire: 
=a 


bzw. 


Mit 


da; 
1 
oe n 2 1 n /\2 a n 2 
r= {|S ari) ar A) —itot (7 aw ms) | db 
i aaa 4 ba k=1 
1 n 
6l* : . , 7 nm n 
oa =2| v > a, Vv, + = Uo*d; > ay Vj, — U GO? > a,v, | dé = 0 
6 k=1 k=1 k=1 
nm : 
it GB oife 1 = Lege nae (6) 2 
ay, VU; Vy ae pet 04 U; vy, — Ug" ¥; my tte es (Ret ase (42) 
k=1 


1 


1 ll 
Tong) . / 
Aj, === J 0” VU; UE dé 5 Ki, — { o} v; UE dé 5 Fy, = Ah v; VE dé 
0 0 


0 


wobei in A;, und K;, nur tiber die Bereiche zu integrieren ist, in welchen 9 = 0, folgt aus (42) 


a Siche Fugnote 2 von Ss: Lit, 
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Die Gleichung n-ten Grades fiir u lautet 


i 1 | 

Fai (4a Kix)/=0 Minis aes Leto. Th) (43a) 

die entsprechenden Gleichungen fiir synchrone Prazession im Gegenlauf und unberiicksichtigte 
Kreiselwirkung A . 3 f 

|Pe—# (Aa + 7 Ku) = 0, (43b) 

| Fin — u Ay,||=0. (43) 


Fir die Lésungen der Ritzschen Gleichungen (43) bis (43b) 1aBt sich unter gewissen mathema- 
tischen Voraussetzungen exakt beweisen!, da® diese fiir unberiicksichtigte Kreiselwirkung und 
synchrone Prazession im Gegenlauf obere Schranken, fiir synchrone Prazession im Gleichlauf obere 
Schranken fiir die positiven, untere Schranken fiir die negativen Kigenwerte sind. Fiir die Lésungen 
der nach der Methode von R. Grammel baw. deren Verallgemeinerung gewonnenen Gleichungen (40) 
bis (40b) gelten folgende Aussagen: 

a) Ohne Kreiselwirkung und synchrone Prazession im Gegenlauf gilt: In denjenigen Fallen, 
in welchen sich die Funktionen des Naherungsansatzes nach den Eigenfunktionen des zugehérigen 
Eigenwertproblems entwickeln lassen, ist unter Verwendung dieser Entwicklung mit Hilfe eines 
von H. Kneser? abgeleiteten Satzes beweisbar, daB die n Wurzeln von (40) bzw. (40b) obere Schran- 
ken fiir die n ersten Eigenwerte darstellen und genauer sind als die mit dem gleichen Naherungs- 
ansatz berechneten entsprechenden Ritzschen Naherungen. 

b) Fiir synchrone Prazession im Gleichlauf ist keine allgemein giiltige Aussage gefunden worden. 

Da jedoch in den praktisch interessanten Fallen die den Beweisen zugrundeliegenden Voraus- 
setzungen nicht erfillt sind, ist man auf die bei der Durchrechnung verschiedener Beispiele, deren 
strenge Losungen bekannt sind, gewonnene Erfahrung angewiesen. Im Falle Ritz* hat sich dabei 
gezeigt, daf} man auch hier obere Schranken fiir die positiven Eigenwerte erhalt. Im Falle Grammel 
konnte in verschiedenen Beispielen die unter a) gemachte Aussage bestatigt werden. Fiir synchrone 
Prazession im Gleichlauf wurden aus (40a) Naherungen gewonnen, die z. T. oberhalb, z. T. unterhalb 
der strengen Lésungen liegen, deren Abweichungen von letzteren dem Betrage nach kleiner sind 
als bei den entsprechenden Ritzschen Lésungen. 

6. Zusammenfassung. In der vorliegenden Arbeit erfolgt die Berechnung der kritischen Dreh- 
zahlen von kontinuierlich besetzten Wellen konstanter Biegesteifigkeit mit einer Verallgemeinerung 
der von R. Grammel entwickelten Methode vermége zweier miteinander aquivalenter [terations- 
vorschriften: die eine bedient sich der Integrodifferentialgleichung, die zweite der Differential- 
gleichung samt Rand- und Ubergangsbedingungen. Bei beiden Vorschriften wird die Wellenaus- 
lenkung ermittelt, die zu der mit einer linearen Kombination zulassiger Funktionen gebildeten 
Belastung gehért. Mit diesen Vorschriften und einer weiteren fiir die auftretenden Variationen 
geht man in die aus dem Hamiltonschen Prinzip folgende Variationsgleichung ein. SinngemaBe 
Ubertragung der Grammelschen Methode fihrt zu einem linearen homogenen Gleichungssystem fiir 
die Unbekannten des Naherungsansatzes, dessen verschwindende Koeffizientendeterminante eine 
algebraische Gleichung fiir die Eigenwerte darstellt. 

Die charakteristische Eigenschaft der Methode von R. Grammel, obere Schranken — und 
zwar genauere als der Ritzsche Verfahren — zu liefern, tritt auch bei Bericksichtigung der Kreisel- 
wirkung im Falle synchroner Prazession im Gegenlauf stets auf. Im Falle synchroner Prazession 
im Gleichlauf dagegen kénnen die Naherungslisungen bei dem verallgemeinerten Grammelschen. 
Verfahren z. T. oberhalb, z. T. unterhalb der entsprechenden strengen Lésungen liegen. 

Unter Beachtung des physikalisch-technischen Sachverhaltes lassen sich die hier abgeleiteten 
Naherungsgleichungen sowohl auf mehrfach gelagerte Wellen mit konstanter Biegesteifigkeit und 
kontinuierlicher Scheibenverteilung als auch auf Wellen (zwei- und mehrfach gelagert) mit ver- 
anderlicher Biegesteifigkeit und Einzelscheiben iibertragen. Das sich im zuletzt genannten Fall 
einstellende Endergebnis unterscheidet sich nur unwesentlich von demjenigen der Methode von 


A. Traenkle4. 
Eine Fortsetzung dieser Arbeit bringt die Durchrechnung eines Zahlenbeispiels. 


(Eingegangen am 26. August 1955.) 
Anschrift des Verfassers: Dr. Ludwig Fehrle, Miilheim/Ruhr, Dickswall 2b ILI. 


1 &. Kamke, Math. Z. 45 (1939), S.759; 46 (1940), S. 231; 48 (1942/43), S. 67. 
2 H. Kneser, Math. Z. 45 (1939), S. 590. . 

3 Siehe FuBnote 1 von S. I11. 

4 A. Traenkle, Ing.-Arch. 1 (1930), S. 499. 
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Stiitzmomenteneinflu&felder durchlaufender elastischer Platten 
mit zwei frei drehbar gelagerten Randern 


Von G. Hoeland. 


1. Einleitung. EinfluBfelder elastischer Platten sind bereits auf die verschiedensten Arten er- 
mittelt worden. Wahrend u.a. Bittner! von bekannten Reihenansatzen ausging, bediente sich 
Pucher® der Singularitatenmethode. Alle diese Arbeiten behandeln jedoch im wesentlichen nur 
Kinfeldplatten. Lediglich Bittner untersuchte naherungsweise den Kinflub der Nachbarfelder. In 
der vorliegenden Arbeit? soll nun ein Weg aufgezeigt werden, um die EinfluBfelder fiir die Momente 
iiber den Zwischenstiitzen durchlaufender elastischer Platten mit zwei frei drehbar gelagerten 
Randern zu ermitteln. Auf die Berechnung der Feldmomente durchlaufender Platten wird nicht 
eingegangen, da bei diesen der Einflu®B der Durchlaufwirkung mit sehr guter Genauigkeit einfach 
bestimmt werden kann. 

Ausgehend von der Singularitét des EinfluBfeldes fiir das Stiitzmoment am Querrand eines 
Plattenhalbstreifens wird die Singularitat des EinfluBfeldes fiir das Moment iiber der Zwischen- 
stiitze zweier aneinandergrenzender Halbstreifen gewonnen. Diese Singularitat bildet dann den 
Ausgangspunkt fiir die Ermittlung der EinfluBfelder symmetrischer und unsymmetrischer Platten- 
paare mit starrer oder elastischer Zwischenstiitze. SchlieBlich wird noch ein Weg angegeben, um auch 
fiir den Fall engliegender elastischer Zwischenstiitzen das StiitzmomenteneinfluBfeld zu ermitteln. 


141 . 


ra 


Abb. la. Es ist a=ay, aes b= by =f. 


2. Das EinfluSfeld fiir das Moment iiher der starren Zwischenstiitze zweier aneinandergrenzender 
Halbstreifen. Wir betrachten zuerst das Stiitzmomentencinfluffeld fiir einen am Querrand ein- 
gespannten Plattenhalbstreifen. Es laBt sich durch die Funktion 


iy = aD e e “*sinausina x (2,1) 
ley ora 

darstellen, wobei w die EinfluBordinaten, x und y die Koordinaten, u die Lage des Aufpunktes (der 

Punkt, fiir den das EinfluBfeld bestimmt werden soll) in x-Richtung, | die Stiitzweite der Platte und 


' E. Bitmer, Momententafeln und EinfluBflachen fiir kreuzweise bewehrte Eisenbetonplatten. Wien 1938. 

” A. Pucher, Die MomenteneinfluBfelder rechteckiger Platten. Deutscher AusschuB fiir Stahlbeton, 
Heft 90. — A. Pucher, Ing.-Arch. 12 (1941), S. 77. — A. Pucher, Ing.-Arch. 14 (1943), S. 246. 

* Gekirzte Fassung einer von der Fakultat fiir Bauwesen der Technischen Hochschule Hannover geneb- 
migten Dissertation, Referent: Prof. Dr.-Ing. habil. A. Pfliiger, Korreferent: Prof. Dr.-Ing. W. Zerna. 
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a =nz/l eine von n abhangige Konstante bedeuten. Das EinfluBfeld wurde zuerst von Pucher! 
angegeben. Abb. 1a ist eine skizzenhafte Darstellung der Umgebung des Aufpunktes als Héhen- 
schichtlinienplan. Dabei sei der Halbstreifen an der x-Achse eingespannt. 

Nun denken wir uns das EinfluBfeld fiir das Stiitzmoment iiber der starren Zwischenstiitze 
zweier aneinandergrenzender Halbstreifen. Die Platte soll dabei kontinuierlich durchlaufen. 
Abb. 1b moge wieder die Umgebung des 
Aufpunktes darstellen. An der Zwischen- * 80 -tach 
stiitze («-Achse) sei die Platte frei dreh- / 
bar gelagert. / ras eave ‘ 

Betrachten wir eine Last im Punkte 4 / | < 

25st. SN : a 
' Pak \ 
iN 


~ 


auf dem Halbstreifen, so muf sie im Auf- PA / Ge 
punkt des Einfluffeldes das gleiche Mo- i i ee 
ment hervorrufen wie zweiLasten gleicher j 
Grée in den entsprechenden Punkten B 
und B’ der Durchlaufplatte im Aufpunkt | ff i 
dieser Durchlaufplatte. AuSerdem mu | / / \ \ 
eine Last im Punkte B das gleiche Mo- | / / 
ment im Aufpunkt heryorrufen wie eine Hs f Sie 
andere Last P gleicher Art und GréBe | / \ \ 
im Punkte B’ des gleichen EinfluBfeldes I| \\ 
der Durchlaufplatte (vgl. Abb. 1b). 

Das EinfluBfeld fiir das Moment iiber j / j ae \ 
der starren Zwischenstiitze zweier an- 
einandergrenzender Halbstreifen hat also me 
die gleiche Form, aber nur die halben 
Ordinaten wie das Stiitzmomentenein- i “Gd 
fluBfeld des am Querrand eingespannten 
Halbstreifens. Es ]a8t sich durch die 25 
Funktion 


WwW = Wy 


~ 
—— > + -e-%Ysinausinax (2,2) 
WN, 2.5% N 


darstellen und ist beiderseits der Zwi- 
schenstiitze gleich. In Abb, 2 ist dieses 
EinfluBfeld als Héhenschichtlinienplan 
angegeben. i 
Die gleichen Uberlegungen, die hier 
fiir das EinfluBfeld der aneinandergren- 
zenden Plattenhalbstreifen angefihrt 
sind, lassen sich analog fiir jedes andere 
symmetrische Plattenpaar und ebenso 
fiir die klassische wie fiir jede andere 
Singularitat des Stiitzmomenteneinflub- 
feldes anstellen. Es ergibt sich daraus 
allgemein, das die EinfluBfelder fiir die 
Momente iiber starren Zwischenstiitzen 
symmetrischer durchlaufender Platten- 
paare aus den Stiitzmomenteneinflufi-_. 
feldern der entsprechenden eingespann- Abb. 2, my-Einflubfeld fiir die Mitte iiber der starren Zwischenstiitze eines 
ten Einfeldplatten gewonnen werden Rens 
kénnen. In vielen praktisch vorkommenden Fallen wird man, wenn die Seitenverhaltnisse zweier 
aneinandergrenzender Platten nicht zu sehr voneinander abweichen, mit guter Naherung auch fiir 
unsymmetrische Platten entsprechend vorgehen kénnen. 


yY 


ZwischensturZze 


1s 


3. Das EinfluBfeld fiir das Moment iiber der Zwischenstiitze zweier aneinandergrenzender Recht- 
eckplatten. In Abschnitt 2 wurde das EinfluBfeld fiir das Stittzmoment zweier aneinandergrenzender 


iA, Pucher, EinfluBfelder elastischer Platten, Tafel 9. Wien 1951. 
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Plattenhalbstreifen hergeleitet. Die sehr einfache Singularitat dieses EinfluBfeldes erfillt bereits 
die Randbedingungen an den Langsrandern, wenn die Platte an diesen frei drehbar gelagert und 
starr gestiitzt ist, sowie die Kontinuitatsbedingungen 


3 02 a a, 
0 Winks 0 Wyechts OW inks OW, echts 
= - 9 = — 

Orvy2 eRe) B a 
oy links oy recbts OV inks CV pechts 


an der Zwischenstiitze. Damit ist die Ermittlung der EinfluBfelder fiir das Moment an der Zwischen- 
stiitze zweier aneinandergrenzender Reckteckplatten auf die Erfillung der zusatzlichen Rand- 
bedingungen an den der Zwischenstiitze gegeniiberliegenden Randern zuriickgefiihrt, wenn man das 
EinfluBfeld (2,2) als Singularitat benutzt. Die Erfiillung dieser zusdtzlichen Randbedingungen 
erfolgt durch biharmonische Funktionen. Dabei ist es auch hier zweckmaBig, solche Funktionen 
zu wihlen, die fiir sich bereits méglichst viele Randbedingungen erfiillen, also Funktionen von der 
Form 


w= > Sign) 1B OX (3,1) 
n=1,2.. 


Funktionen dieser Art sind: 


w= » (AnGofay+ Bpzy Ginay + C, Sina y + Din yCojay)sinax. (3,2) 
n=1,2.. 
Dabei bezeichnen i das Plattenfeld und n das Reihenglied. 
Fir die Randquerkraft ergibt sich mit N als Plattensteifigkeit unter Vernachlassigung der 
Querkontraktion 


oder 

N & [An(—a) Sina y + Br (a? Sina y —a@ y Cofj ay) 
“2. (3,3) 

+ Cj, (— a) Coj ay + Din (@ Gof ay —a’ y Sina y)] sina x. 

Die Randbedingungen an den Gegenrandern kénnen beliebig vorgegeben sein. Haben wir 

jeweils einen eingespannten Rand im Abstande y, = b, und y, = by, so lauten sie 


Th fur cy, = 0 und “y7—b). (3,4), (3,5) 
Ow a 
aye =) fT 5b, (3,6) 
fir die eine Rechteckplatte mit w, = wy, + w, und 
W, = 0 fir ya=-0 “und styp==be, (3,7), (3,8) 
Ow, ~ 
ay» = 0 fiir Mol = by (3,9) 


fiir die zweite Rechteckplatte mit w. = wo, + w,. Dabei bezeichnen w, und w, die EinfluBordinaten 
in den Feldern 1 und 2, wo, und wo, die Singularteile [entsprechend (2,2)] sowie w, und w, die Regular- 
teile [entsprechend (3,2)] in diesen Feldern. 

An der Zwischenstiitze lauten die Kontinuitatsbedingungen 


@ on, 

c= eG it Vy = Vo = 0 tp == (h—= my (3,10) 
V1 OY» zs s 1 2° ’ 

Pw, Cou ane ( muy 

at oy a Sa 0| xy =1— xp. (3,11) 


Die Randbedingungen und die Kontinuitatshedingungen werden fiir jedes Reihenglied n einzeln 
erfillt. Die Singularitat selbst erfiillt die Bedingung (3,10) nicht fiir jedes Reihenglied einzeln son- 
dern nur fiir die ganze Summe. Bei der Ermittlung der Konstanten geniigt es aber, bei dieser Be- 
dingung die Regularteile alleine zu beriicksichtigen. 

Die Randbedingungen (3,4) und (3,7) fiihren auf 4;, = 0, da die Singularitat fiir sich bereits 


die Randbedingung erfiillt und die wbrigen Funktionen (die Glieder mit B;,, C;, und D;,) eben- 


falls fiir Y1 = Y¥2 = 0 schon Null sind. Damit ist fiir jeden Wert n ein System von 6 Gleichungen 
mit 6 Unbekannten zu lésen. 
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Fur den Sonderfall u — 41 fallen die Glieder mit geradem n heraus und es sind nurmehr die 


Konstanten fiir n = 1, 3, 5,... zu berechnen. Fiir den Fall b, =41 und b, =1 wurde das Ein- 
fluBfeld berechnet. Es ergaben sich dabei folgende Werte fiir die Konstanten (w= 41): 
y= Bo — ©0,2364)1 , y 
Ce == 02310, i Sl aac Ice cam Des 
fire = Neat, | 
Gy = 7 0.0297), 
Ds, = — 0,2614/I, 
Pigs a. =! 05585) 1058/1"; 
Cys = 0, (02-1073 > 
Dis eA L021ss 10-9; 
Cos = — 0,558- 10-3, 
De. = — 0,022- 10-6/1. 


Das entsprechende Einflu®feld ist in Abb. 3 
dargestellt. 

Grundsatzlich der gleiche Weg ist ein- 
zuschlagen, wenn das EinfluBfeld fiir eine 
Zweifeldplatte zu bestimmen ist, bei der ein 
oder beide Gegenrader nicht eingespannt 
sondern frei drehbar gelagert sind. Es sind 
dann lediglich die entsprechenden Rand- 

bedingungen zu erfiillen. 


4. EinfluBfeld fiir das Moment iiber einer 
einzelnen elastischen Zwischenstiitze. Haben 
wir einen Plattenstreifen mit einer einzelnen 
elastischen Querstiitze und wollen wir das 
‘EinfluBfeld fir das Moment iiber dieser 
Querstiitze ermitteln, so gehen wir am 
besten wieder vom EjinfluBfeld der an- = 
einandergrenzender Halbstreifen (Abb. 2)  *””Liicr"Garchlaufenden Platte. Seitenverhiltnis Til und 1:03. 
aus. Fir dieses EinfluBfeld gilt (2,2). 

Die beiden Seiten iiben auf die zuerst starr angenommene Zwischenstiitze eine Belastung aus, 
die gleich der Randquerkraft der-beiden Felder ist. Da jetzt die Stiitze elastisch sein soll, erleidet 
sie unter diesen Randquerkraften q, eine Verformung, der die Platte wiederum anzugleichen ist. Die 


-Randquerkraft 


1 Pw q Ow ae 
y= P =— Noe + 2a are (vy = 0) 
-wirkt auf beiden Seiten gleichgerichtet, ist also zu addieren. 
An der Zwischenstiitze ergeben sich die folgenden Rand- und Kontinuitatsbedingungen: 


WQu.Tr. = Wy = Wy fir y,=y.=0, (4,1), (4,2) 
Cw, Ow, a 
——$ = — = = fur Mi = Yo as 0 5 (4,3) 
OY, OY» . 
Pw, Fw, a = a 4.4: 
Ba Os fOr yp V3 0 (4,4) 


Dazu kommen noch die Randbedingungen an den der Zwischenstiitze gegentiberliegenden Randern, 
wenn solche vorhanden sind. Die Stiitzungen an diesen Randern kénnen elastisch oder starr sein. 

Betrachten wir den Sonderfall, daB das EinfluBfeld fiir die Mitte tiber der Querstiitze zweier 
aneinandergrenzender Halbstreifen gesucht ist (im anderen Fall ware analog zu Abschnitt 2 und 3 
zu verfahren), so wird A, =— C, und B, =— D,, (der Index i kann hier entfallen) und daraus 


wes Dy (A, e—* + D, y e—*%) sina x. (4,5) 


n=1, 3... 
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Fiir den elastischen Quertriiger ist p = 2 (qyo ++ gy). Dabei stellen q,, den Anteil aus der Singularitat 
w, und q, den Anteil aus dem Regularteil w dar. Durch Einsetzen der entsprechenden Ableitungen 
erhalt man fiir y = 0 


p=—2N >) (— + sin au —— A, a = DD; a) sin a x (4,6) 


n=1,3.. 


1 und aus der Differentialgleichung 


der Biegelinie des Quertragers 


ao 


sap 2 
g Wou.Tr. : —2N » Y a’ eS, 
Ox as FE Jou. ar: ; l 
n=153 6. 


+ A,a? + D, «| sinax. (4,7) 


Durch Integration ergibt sich aus dieser 
Differentialgleichung die Biegelinie des 
Quertragers, wenn die beiden Langsrander 
der Platte (und damit auch die Enden des — 
Quertragers) frei drehbar gelagert sind, zu 


== 2'Nin tm ait Pop. 
nf Tra) = S——_——_— 
WQu.Tr. E Jou. Ty. 2 , e a7 ina 


1 I Ns 
| A, he) sin a x 


und mit a =nz/l 


_—2Nl YN /(—1. 
WQu.Tr. = E Jqu.tr. = 7 sin au 


1 
= A, =~ + D 


) sin ax. (4,8) 


n n2 7 
Setzen wir zur Abkiirzung 
2 N1/ EJqucrr. SS '5) 


so erhalten wir aus (4,8) 


ae Be 
We Our ae ee sin au + = 
w= SiS. 
Dane 
=i Dot) sin ax. 
Die Randbedingung (4,1) liefert 
oo ae | ° An De \ ; 
2 (S5 BUM ON aia ! sinax 
aa S A, siihax — 0% (4,9) 
n=1,3.. 


Abb, 4. my-EinfluBfeld fiir die Mitte tiber der elastischen Zwischenstiitze —— ° : ‘ iz n . 
if eines Plattenstreifens, : Dabei 1st die Bedingung (4,2) selbstandig 
miterfiillt, da wir vom symmetrischen Plat- — 
tenfeld ausgingen. Fiir den gleichen Sonderfall ergibt sich aus (4,3), wenn wir beriicksichtigen, dab 
die Singularitat diese Kontinuitatsbedingung bereits fiir sich erfiillt, 


Pay (— Ana + D,) sin €%:== 0 ez = fir a= } 


und damit die einfache Bedingungsgleichung 
Ala = Dr. (4,10) 
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Setzen wir nun in (4,9) ein, so erhalten wir 


~ 


oO 
Y ‘ y : J 
——, sinau — 2 A, — }sinax — Auer Hees 
2; & a re x A, sinax = () 


SS a 


nm 


Fie us sin au in a 
i nn (a. + 2/n a) (na + 2y)na (4,11) 
sowie mit (4,10) 
ite y sin au 
a (na+2y)1 (4,12) 


als Bedingungsgleichungen fiir A, und D,. In Abb.4 ist ein solches EinfluBfeld fiir y =1 dargestellt. 
Der Grenzfall eines vollkommen weichen Quertragers (E Jgu.rr, =) liefert y = oo. Damit 
ergibt sich fiir w entsprechend w = w,-+ w [Gleichung (2,2) und (4,5)] 


—— »3 & * sinau-+ 4, + Day) e— ¢Y sin ax (4,13) 
n=1,3..° 
mit Ens | AL a snau — = sin au 
ae ‘  na(na+2y) 2n7 
und nD. = y sin au he sin au 
"yr 00 “ (nz -+ 24) 2] 
inn 1 = 1 re lies Mee we ; 
ia Ix (Lay) - e sinausinax. (4,14) 
n=1,3 


Die Gleichung (4,14) stimmt, abgesehen vom zweiten Vorzeichen, mit der von Pucher! an- 
gegebenen Lésung fiir das EinfluBfeld m, fiir die Mitte des Plattenstreifens iiberein. Wir erhalten 
also im Falle des Grenziiberganges E Jo,.7,, > 0 das EinfluBfeld des Plattenstreifens. 


Bererch3 


Abb. 5. m bezeichnet den Punkt, fiir den das Einflufifeld ermittelt werden soll. 


5. Einfluffeld fiir das Moment iiber einer von vielen engliegenden elastischen Zwischenstiitzen. 
In diesem Abschnitt soll ein Weg aufgezeigt werden, um angenahert das EinfluBfeld fiir das Moment 
tiber einer elastischen Zwischenstiitze zu ermitteln, wenn in gewissen Abstanden weitere elastische 
Zwischenstiitzen vorhanden sind (orthotrope Platte mit einer Schar Trager). Die Abb.5 mége ein 
solches System im Grundrif darstellen. Das EinfluBfeld soll fiir den Punkt m ermittelt werden. 

Zur Berechnung des EinfluBfeldes werden die ersten beiden Plattenfelder (Bereiche ] und 2 
in Abb.5) beiderseits des Quertragers, auf dem der Aufpunkt m liegt, als isotrope Platte betrachtet. 
Fiir sie gilt also die Differentialgleichung AAw=p/N (mit p=0). Der weitere Bereich (Bereich 3 
nach Abb.5) wird als Kontinuum betrachtet. Fiir diesen Bereich mége die Hubersche? Differential- 
gleichung 


tw Aw Ow 


(auch hier ist p=0) Giiltigkeit haben. Es bedeuten darin B, und By die Biegesteifigkeiten in den 
beiden Koordinatenrichtungen und H die Drillsteifigkeit (im Sonderfall y—0) der orthotropen 
Platte. 


) 


1) A. Pucher, Ing.-Arch. 14 (1943) S. 259. 
2) Vel. K. Girkmann, Flachentragwerke, S. 292. Wien 1954. 
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Die einzelnen Plattenfelder sind an den Randern x =0 und x=/ starr gestiitzt und frei drehbar 
gelagert. An den Zwischenstiitzen sind sie elastisch gestiitzt und in die Nachbarfelder elastisch 
eingespannt. Am ersten Plattenfeld wirkt die Singularitat (2,2). Sie ruft an den Querrandern Ver- 
formungen hervor, denen die Nachbarfelder anzupassen sind. Dabei miissen die folgenden Rand- 
und Kontinuitatsbedingungen erfillt werden: Uber dem Quertrager 1 (y, = 0) 


(5,2) 


(5,3) 


W Platte 1 = Wau. Trak 


OW pili Cw Pl 
wobei die Kontinuitatsbedingung 5 = 5 : 
Yi a 


sprechenden Ansatzes gleich miterfillt ist; iiber dem Quertrager 2 (y, =6, vy, = 0) 


fiir y,;—=0 infolge des der Symmetrie ent- 


Wy = U2 ; 5,4) 
Wz = WQu.Tr.2 > (5,5) 
Pee OO 
OY We ROV ee (5,6) 
Ow, Cw, Saat 
N, ay = Maye (meistiist NG Ney) (5,7) 
tiber dem Quertrager 3 (vy. = b, v3 = 0) 
Wy = W3, (5,8) 
Wz = WQu.Tr.3 > (5,9) 
aE cae 510 
Oy. Os” (>, 
OW, Cws 
2 oy? ney dy? g (5,11) 


Zur Erfiillung dieser Rand- und Kontinuitaétsbedingungen werden wieder entsprechende Funk- 
tionen und deren Ableitungen benétigt. w, besteht wieder aus einem Singularteil und einem — 
Regularteil, w. und w, dagegen bestehen nur aus Regularteilen. Singular- und Regularteile der 
Felder J und 2 haben die gleiche Form wie im Abschnitt 3. Entsprechend kénnen sie auch durch 
die gleichen Funktionen dargestellt werden. Die Erfiillung der Huberschen' Differentialgleichung (5,1) 
der orthotropen Platte im Bereich 3 dagegen erfolgt durch Funktionen, wie sie bereits von Corne- 
lius? benutzt wurden: 


we= 2 (A;3, cosyays + Bs, singays)e°*%ssina x 


n=1,2.. 


Nh eshia coe een is ee 
mit « = a) B,/B, und ot Vel type) NaN 


y 


Fir die Randquerkraft qy ergibt sich bei der orthotropen Platte? 


aw Boye 
dy = B 4:6 229. RB.) 
ts aroye ( tee) Ox? Oy 
und mit y, =v, =0 
Aw Ow 
qy = Be uv 2 0 we : 
y of oy? Ox oy 


as wir die entsprechenden Ableitungen in die Gleichung fiir q, ein, so wird fiir den Sonderfall ’ 
PS | 


h3 = = [43, (— Ou — 2 Oy?) + B3, (up —2P)] a? | BS B, sin ax. 


‘ K. Girkmann, Flachentragwerke, S. 292. Wien 1954. 
2 W. Cornelius, Der Stahlbau 21 (1952), S. 21, 43 und 60. 
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Durch Einsetzen der Ableitungen in die Randbedingungen (5,2) bis (5,11) erhalten wir die Bedin- 
gungsgleichungen zur Ermittlung der Konstanten 4;,, Bin, C;, und D;,. Das EinfluBfeld setzt 
sich dann im Bereich 7 aus dem Singular- und dem Regularteil zusammen wihrend es in den 
Bereichen 2 und 3 durch die Regularteile alleine dargestellt wird. 


6. Vergleich mit dem Naherungsverfahren der DIN 1075. In der DIN 1075 (Fassung vom August 
1951) wird gefordert, Fahrbahnplatten nach der Plattentheorie zu berechnen. Fir die durch- 
laufenden Platten werden dabei Naherungsverfahren angegeben. Nun ist es fiir die Praxis von 
Interesse zu erfahren, wie genau diese Naherungsverfahren sind. Es wurden deshalb die Stiitz- 
momenteneinflufBfelder fiir symmetrische Plattenpaare verschiedener Seitenverhaltnisse L,/ly aus- 
‘gewertet und mit den Werten verglichen, die sich nach den Verfahren der DIN 1075 (ohne Faktor a 
nach Tabelle 3 der DIN 1075) ergeben. Die nach dem Naherungsverfahren ermittelten Stiitz- 
momente Mz, waren bis zu 75°, gréRer als die nach dem exakten Verfahren ermittelten Werte Ms, 
der Durchlaufplatte. 


@ 
18 
16 eee LEA me 
AS Ly MZ. 
y Yy 
14 - : fe 
WWE, 
| Ze He oy Ve 
| EE 
| Zi, 
42 
| 
40 U 
: 30 50 70,0 75,0 200™ 
ME ; : F : 
Abb. 6. @ = Mo,’ wobei My = Einspannmoment der Einfeldplatte, 
t 


Ss 
Ms¢ = Stiitzmoment der Zweifeldplatte; Briickenklasse 60. 


Die Abb.6 gibt einen Uberblick ttber die GréBe der Abweichungen in Abhangigkeit von der 
Stiitzweite der Platten. Alle ermittelten Werte fir die Briickenklasse 60 lagen innerhalb des 
schraffierten Bereiches der Abb.6. Die Darstellung zeigt, daB die Abweichungen bei Stiitzweiten 
von 4 bis 10 m besonders groB sind. Es dirfte deshalb zweckmafig sein, die in der Vorschrift 
angefihrten Naherungsverfahren zu erganzen und die Vorschrift entsprechend abzuandern!. 


7. SchluBbemerkung. Es wurden Verfahren gezeigt, um StitzmomenteneinfluBfelder durch- 
laufender elastischer Platten mit zwei frei drehbar gelagerten Randern zu berechnen. Ausgehend 
vom EinfluBfeld des Plattenhalbstreifens wurden die Gleichungen zur Ermittlung der EinfluB- 
felder iiber einer starren Zwischenstiitze aufgestellt und ein solches EinfluBfeld berechnet. Bei der 
Ermittlung dieser EinfluBfelder spielt es fiir den Rechenaufwand grundsatzlich keine Rolle, ob 
die Platten an den Gegenrandern frei drehbar gelagert, elastisch oder fest eingespannt, elastisch 
oder starr gestiitzt sind oder auch freie Gegenrander haben. Es sind nur jeweils die entsprechenden 
Randbedingungen zu erfiillen. 

Fur den Sonderfall der zur starren Zwischenstiitze symmetrischen Plattenfelder ergibt sich 
eine vereinfachte Lésung derart, da} lediglich die Ordinaten der EinfluBfelder der entsprechenden 
eingespannten Platten zu halbieren und beiderseits der Zwischenstiitze aufzutragen sind. Dieses 
Verfahren kann jedoch nicht angewandt werden, wenn die Platte mit einer verdrehungssteifen 
Zwischenstiitze fest verbunden ist. 

Weiter wurden Lésungen fiir die Ermittlung von Einfluffeldern fiir die Stiitazmomente tber 
elastischen Zwischenstiitzen angegeben. Dabei wurden sowohl einzelne als auch engliegende Zwischen- 
stiitzen beriicksichtigt. 


1 Naheres hieriiber ist in der Dissertation des Verfassers zu finden. 
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die einzelnen Plattenfelder durch Befriedigung der Kontinuitétsbedingungen an den Zwischen- 
stiitzen zu verkoppeln, wahrend bei den Stiitzmomenten die Singularitaten mit in die Kontinuitats : 
bedingungen eingehen. 

Zum Schlu8 wurden Vergleiche mit den bisherigen Naherungsverfahren gefiihrt. Dabei wurde 
gezeigt, das sie im allgemeinen, besonders bei groBen Einzellasten, zu ungiinstige Werte liefern, 
Es wurde deshalb auf einen Abinderungsvorschlag fiir die DIN 1075 hingewiesen, durch den die 


in dieser Arbeit gewonnenen Erkenntnisse Beriicksichtigung finden kénnen. 


(Eingegangen am 13. September 1955) 
Anschrift des Verfassers: Dr.-Ing. G. Hoeland, Dortmund, Liebigstr. 4 ~ 
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Acht- und neunstellige Tabellen zu den ellip- 
; tischen Funktionen. Dargestellt mittels des Jacohischen Parameters q. 


Eight und Nine Place Tables of Elliptical Fune- 
tions. Based on Jacobi’s Parameter q. | : - 


Von Dr.-Ing, M. Schuler, Professor emeritus an der Universitat Gittingen, und Dr. phil. habil. H. Gebe- 
lein, Dozent, Bamberg. Mit einem englischen Text von LauritzS.Larsen, B. S. Mit 11 Ab- 
bildungen. XXIV, 296 Seiten 4°, 1955. Ganzleinen DM 58,— 


Inhaltsverzeichnis (deutsch und englisch): Vorwort. Inhaltsverzeichnis. Einfihrung. Ta- 
belle I: Funktionen G (q*, z) laufend nach z (mit Angabe der zugehérigen Werte q und @). Tabelle Tl: 
Funktionen G (q#, =) laufend nach q*. Zugehérige Werte q und O mishmiiongectellt. Tabelle IIT: Funk- 
tionen H (q°, 2) laufend nach z (mit Angabe der zugehérigen Werte g und @), Tabelle IV: Funktionen 
H (q°, z) laufend nach q°, Zugehérige Werte q und O.zusammengestellt. Tabelle V: Jacobische elliptische ‘ 
Funktionen laufend nach z (mit Angabe der zugehérigen Werte fiir O, —Ig cos O = —Ig.k’, K und K/E). 
Tabelle VI: Jacobische elliptische Funktionen laufend nach q. Zugehirige Werte fiir —lg k’ und 0. Zu- 
gehérige Werte fiir K und K/£. Tabelle VII: Tafeln fiir die Umrechnung zwischen dem Legendreschen 
Modul © und dem Jacobischen Parameter g. 


Fiinfstellige Tabellen zu den elliptischen Funk- 


e 
tionen. Darsgetellt mittels des Jacobischen Parameters q. 


Five Place Tables of Elliptical Funktions. 


Based on Jacobi’s Parameter q. 


Von Dr.-Ing. M. Schuler, Professor emeritus an der Universitat Gottingen, und Dr. phil. habil, H. Gebe- 
lein, Dozent, Bamberg. Mit einem englischen Text von Lauritz S. Larsen, B.S. Mit 11 Ab- 
bildungen. XI, 114 Seiten Gr.-8°. 1955. Ganzleinen DM 29,60 


Inhalts “ erzeichnis (deutsch und englisch): Vorwort. Inhaltsverzeichnis. Finfithrung. Ta- 
belle I: Jacobische elliptische Funktionen laufend nach z mit Angabe der zugehorigen Werte ©. Tabelle IT: 
Jacobische elliptische Funktionen laufend nach q. Beigefiigt: Werte fiir © und —Ig cos ©. Werte fiir K 
und K/E. Tabelle III: Funktionen G(q z) und H (2) laufend nach z (mit Angabe der zugehérigen 
_Werte O). Tabelle IV: Funktionen 6 (q,2) und H (q, z) laufend nach q. Tabelle V: Tafel fiir die Umrech- 


nung zwischen dem Lrgendreschen Modul © und dem Jacobischen Parameter q. Tabelle VI: Tafeln der 


f Koeffizienten fiir die Interpolation nach Everett. i 
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Berechnung und Gestaltung yon 


G ifed 
mS 
Von Dr.-Ing. E. F. Gobel, Karlsruhe. (Konstruktionsbiicher. Hrsg. Prof. Dr.-Ing. K. Kollmann, 
Karlsruhe. th Band.) ; . 
a 
Zweite, neubearbeitete und erweiterte Auflage. Mit 110 Abb. VI, 86 Seiten Gr.-8°. 1955. DM 9— 


Inhaltsiibersicht:; Vorwort. Benutzte Formelzeichen. — Einfiihrung: Gewinnung und Verar- 
beitung von Gummi. Herstellung von Gummifedern. Verwendungsméglichkeiten. — Eigenschaften: 
Elastische Verformungsfahigkeit. Darnptaug von mechanischen Schwingungen. Dampfung von Kérper- 
schall-Schwingungen. Schubmodul und Elastizitatsmodul. Shore-Harte. Festigkeit. Einflu8 von Warme, 
Kalte und angreifenden Mitteln. — Berechnung statisch beanspruchter Gummifedern: Berechnungsgrund- 
| lagen. Schubbeanspruchte Gummifedern. Druckbeanspruchte Gummifedern. Kombinierte Schub- und 


Druckbeanspruchung. Sonderfalle. — Berechnung dynamisch heanspruchter Gummifedern: Grundlagen 


der Schwingungsisolierung. Dynamische Federkonstante und Eigenschwingungszabl. Berechnung einer 


gummigefederten Freischwinger-Siebmaschine. Berechnung eines Vibrationsrostes mit gummigefederter a: 


Zusatzmasse. — Konstruktionsbeispiele : Konstruktionsrichtlinien. Allgemeiner Maschinenban. Werkzeug- 


maschinenbau. StraBenfahrzeuge. Schienenfahrzeuge. Cummikupplungen. CGeratebau, Fertigungs- 


technik. — Priifung yon Gummi und Gummifedern. — Schrifttum. — Quellennachweis. 


Seit dem Erscheinen der ersten Auflage hat die Verwendung von Gummi in der Form von Gummifedern 


betrachtlich zugenommen. Es wurden weitgehende theoretische Erkenntnisse gewonnen und viele prak- — 


tische Erfahrungen gesammelt. Der Begriff Gummifeder hat sich analog dem Begriff Metallfeder all- 
gemein eingefihrt. 


{ 


In der vorliegenden zweiten Auflage ist das bis heute Erreichte dargestellt worden mit besonderer Beriick- 


sichtigung der Belange des Konstrukteurs. Die Berechnungsgrundlagen wurden durch einige wichtige 


Falle erweitert und auf schwingungstechnische Aufgaben ausgedehnt. Beispiele aus der Praxis, die zahlen- 
, maBig vollstandig durchgerechnet sind, zeigen die Anwendung der Gleichungen. Der Abechuits Kon- 


struktionsbeispiele konnte erheblich vergréBert werden. Neu hingugekommen sind die Konstruktions- 


richtlinien und Ausfiihrungen iiber die Priifung von Gummifedern bei statischer und dynamischer Bean- 


spruchung. Die wichtigste, in den letzten 20 Jahren erschienene Literatur wurde in einem Verzeichnis 
zusammengestellt, ebenso die von der einschlagigen Industrie und von wissenschaftlichen Instituten in 


dankenswerter Weise zur Verfiigung gestellten Bild- und Tabellenunterlagen. 
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